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Chapter  1 
Introduction 


Chapter  1  contains  a  statement  of  the  waveguide  mode  converter  problem. 
This  problem  is  that  treated  in  [1-3]  generalized  to  the  case  where  an  arbi¬ 
trary  number  or  waveguide  modes  can  propagate  in  the  circular  waveguide 
and  where  the  excitation  is  an  arbitrary  mixture  of  the  forward  traveling 
waves  of  these  modes.^  The  general  method  of  solution  is  the  same  as  that 
described  in  Section  1.2  of  (3).  A  computer  program  was  written  to  imple¬ 
ment  this  solution.  This  program  was  obtained  by  modifying  the  program 
that  is  described  and  listed  in  (3).  ChiH>ter  2,  which  is  based  on  Chapter  2  of 
[3],  contains  instructions  for  using  the  computer  program.  These  instructions 
should  enable  one  to  use  the  program.  In  other  words,  th^  should  enable 
one  to  modify  the  input  data  to  suit  his  needs,  to  run  the  program  with  the 
modified  input  data,  and  to  interpret  the  results. 

The  computer  program  consists  of  a  main  program,  the  subroutines 
MODES,  BESIN.  BES,  INTERPOL.  PHI.  and  DON,  the  function  subpro¬ 
gram  FXY,  and  the  subroutines  DECOMP  and  SOLVE  in  that  order.  The 
subprograms  MODES.  BESIN.  BES.  I.NTERPOL,  PHI,  DON,  FXY,  DE¬ 
COMP,  and  SOLVE  are  exactly  the  same  as  those  described  and  listed  in 
Chapters  4  to  11  of  [3].  The  mam  program  of  the  present  report  was  ob¬ 
tained  by  modifying  the  main  program  in  [3|.  The  formulas  used  in  the  main 
program  to  compute  the  excitatioo  vector,  the  normalized  coefficients  of  the 
propagating  modes  in  each  rectangular  waveguide,  the  normalized  coefficients 
of  the  propagating  modes  in  the  circular  waveguide,  and  the  two  normalized 

Un  [1-3],  only  tlw  even  and  even  sad  odd  TEn  modes  can  propagate  in  the 
circular  waveguide,  and  the  esciMiwo  m  the  forward  traveling  wave  of  the  even  TMoi 
mode.  The  amplitude  of  thw  waw  •  unity 
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<»thogoiial  components  of  the  dectric  field  in  eadi  apertuxe  are  presented  in 
Chapters  3  to  6,  respectively.  The  differences  betneen  the  main  program  of 
the  present  report  and  the  main  program  in  [3]  are  described  in  Chapter  7. 
The  entire  main  program  is  listed  in  Cluster  7. 

He  who  wants  merely  to  use  the  computer  program  is  advised  to  read 
only  Chapters  1  and  2  of  the  present  report  and  to  ccmsult  Ch^ters  1  and 
2  of  [3]  only  when  referred  there.  The  user  will  have  to  go  into  Chapters  3 
to  7  of  the  present  report,  and/or  Chapters  3  to  11  of  [3]  only  if  he  wants  to 
modify  the  computer  program  or  if  unforeseen  difficulty  arises  while  running 
the  program.  If  such  difficulty  does  arise,  the  basic  theory  in  [1]  and  [2]  could 
be  helpful. 

1.1  Statement  of  the  Problem 

There  is,  as  shown  in  Fig.  1,  a  circular  waveguide  that  is  doaed  at  one  end. 
Two  symmetrically  placed  apertures  in  the  lateral  wall  of  this  waveguide  are 
backed  by  rectangular  waveguides  of  identical  dimensions  but  terminated 
with  loads  that  may  be  different.  The  interiors  of  the  left-hand  rectangu¬ 
lar  waveguide,  the  right-hand  rectangular  waveguide,  and  the  dxcular  wave¬ 
guide  are  called  regions  1,  2,  and  3,  respectively.  Homogeneous  space  of 
permeability  41  and  permittivity  e  exists  in  all  of  these  repons.  The  circular 
waveguide  is  of  radius  a  and  is  terminated  by  a  perfectly  conducting  wall 
at  2  s  L3.  An  arbitrary  number  of  modes  can  propagate  in  the  circular 
waveguide.  Both  rectangular  waveguides  run  parallel  to  the  x-axis.  Both 
have  the  same  cross  section  (— |  <  y  <  — f  <  r  <  f)  where  c  <b  and  b 

is  such  that  only  the  TEio  dominant  mode  can  propagate  in  each  rectangu¬ 
lar  waveguide.  As  seen  from  the  circular  wav^ide,  the  aperture  Ai  whidi 
feeds  the  left-hand  rectangular  waveguide  in  Fig.  1  is  the  surface  for  which 
(p  s  a,  X  —  ^0  <  ^  <  T  +  — I  <  2  <  f )  where  p  and  ^  are  the  cylindrical 
coordinates  related  to  x  and  y  by 
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As  seen  from  the  drcuUr  waveguide,  the  ^>ertttie  Aj  which  feeds  the  ri^t- 
hand  rectangular  waveguide  is  the  surface  for  whidi  (p  =  a,  — ^  4  ^ 
dot  — I  <  2  ^  §)•  However,  as  seen  by  the  rectangular  waveguides,  the 
curved  apertures  Ai  and  A3  are  ^proximated  by  plane  iq>ertures. 

The  excitation  is  most  generally  a  mixture  of  the  r-traveling  (traveling 
in  the  x~  direction)  waves  of  all  the  propagating  mocks  of  the  circular  wave* 
guide.  This  mixture  is  specified  by  giving  the  time-average  power  of  each 
propagating  mode  and  its  phase  at  r  s  I13.  By  assun^tion,  the  z-traveling 
waves  at  z  3=  —c/2  are  exactly  this  mixture.  Concerning  the  fields  in  the 
part  of  region  3  for  which  z  >  —c/2  and  in  all  of  repcms  1  and  2,  it  does 
not  matter  how  this  mixture  of  z-traveling  waves  is  produced.  We  let  it 
be  produced  by  a  sheet  of  electric  current  located  at  z  <  —c/2  in  the 
circular  waveguide.  For  simplicity,  we  assume  that  the  circular  waveguide  is 
terminated  at  z  <<  —c/2  by  a  matched  load,  i.e.,  any  —z-traveling  wave  in 
the  region  for  which  z  <  —c/2  is  never  reflected. 

The  voltage  to  current  ratio  of  the  TEio  mode  at  z  s  in  region  1  is 
taken  to  be  Zx.  All  other  rectangular  waveguide  modes  are  evanescent.  The 
voltage  to  current  ratios  of  the  evanescent  modes  at  x  ss  — do  not  come  into 
play  because  Lx  is  taken  to  be  so  large  that  any  evanescent  wave  emanating 
from  the  termination  at  x  s  — £|  will  have  negligibly  small  amplitude  upon 
arrival  at  aperture  Ai.  The  voltage  to  current  ratio  of  the  TEio  mode  at 
X  ss  L3  in  region  2  is  taken  to  be  Z3.  Here,  Lj  is  taken  to  be  so  large  that 
any  evanescent  wave  emanating  from  the  termination  at  x  =  L2  will  have 
negligibly  small  amplitude  upon  arrival  at  aperture  A3. 

The  problem  is  to  find  out  bow  much  time-average  power  of  each  prop¬ 
agating  mode  is  reflected  in  the  circular  waveguide,  how  much  time-average 
power  is  transmitted  into  the  rectangular  waveguides,  and  what  are  the  mag¬ 
nitudes  of  the  ^  and  z-components  of  the  electric  field  in  the  apertures  Aj 
and  A3. 


1.2  Method  of  Solution 

The  method  of  solution  is  that  described  in  Section  1.2  of  [3].  In  this  method, 
the  apertures  are  closed  with  infinitely  thin  perfect  conductors  and,  as  shown 
in  Fig.  2,  a  magnetic  current  is  placed  on  each  side  of  each  of  these  conduc¬ 
tors. 
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Fig.  2.  Top  and  side  views  of  the  situation  equivalent  to  that  of  Fig.  1. 
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Chapter  2 


Instructions  for  Using  the 
Computer  Program 


Written  in  FORTRAN,  the  computer  program  is  available  on  diskette.  On 
disl«tte,  the  computer  program,  which  consists  of  a  main  program  and  some 
subprograms,  is  stored  in  the  file  DEC.92.  Exemplary  input  data  are  stored 
in  the  file  DEC92.DAT.  These  files  are  named  DEC.92  and  DEC92.DAT  be¬ 
cause  December  1992  is  the  date  of  the  present  report,  the  report  in  which 
the  computer  program  is  described  and  listed.  The  main  program,  the  sub¬ 
routines  MODES,  BESIN,  BES,  INTERPOL,  PHI,  and  DON,  the  function 
subprogram  FXY,  and  the  subroutines  DECOMP  and  SOLVE  are  stored  in 
order  in  the  file  DEC.92 

There  are  two  modules  of  input  data  in  the  file  DEC92.DAT.  The  first 
module  of  input  data  is  preceded  by  the  two  comment  statements 

C  THIS  IS  THE  FIRST  MODULE  OF  INPUT  DATA. 

C  PUT  IT  IN  THE  FILE  IN. DAT. 

and  followed  by  the  comment  statement 

C  THIS  IS  THE  SECOND  NODULE  OF  INPUT  DATA. 

The  second  module  of  input  data  is  preceded  by  the  two  comment  statements 

C  THIS  IS  THE  SECOND  NODULE  OF  INPUT  DATA. 

C  PUT  IT  IN  THE  FILE  BESIN.DAT. 
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The  last  line  of  the  second  module  of  input  data  is  the  last  line  in  the  file 
DEC92.DAT. 

To  use  the  computer  program,  first  follow  the  instructions  given  in  the 
comment  statements  in  the  preceding  paragraph,  i.e.,  create  input  data  files 
named  IN.DAT  and  BESIN.DAT,  put  the  first  module  of  input  data  in 
IN.DAT,  and  put  the  second  module  of  input  data  in  BESINJDAT.  The 
data  in  INJ)AT  are  read  by  statements  in  the  main  program.  The  data 
in  BESIN.DAT  are  read  by  statements  in  the  subroutine  BESIN.  Next, 
create  output  data  files  named  OUT.DAT  and  BESOUT.DAT.  Then,  pve 
the  command  or  commands  necessary  to  run  the  computer  pit^am  that 
resides  in  DEC.92.  Running  the  program  causes  output  data  to  be  writ¬ 
ten  in  the  output  data  files  OUT.DAT  and  BESOUT.DAT.  The  data  in 
OUT.DAT  are  written  by  statements  in  the  main  program.  The  data  in 
BESOUT.DAT  are  written  by  statenwnts  in  the  subroutine  BESIN. 

The  contents  of  the  first  module  of  input  data  that  the  user  was  instructed 
to  take  from  the  file  DEC92.DAT  and  put  in  the  file  IN.DAT  are  merety 
exemplary  data.  One  can  modify  the  input  data  in  the  file  IN.DAT  to  suit 
his  needs.^  The  input  data  in  the  file  IN.DAT  are  described  in  Section  2.1.1. 
The  contents  of  the  second  module  of  input  data  that  the  user  was  instructed 
to  take  from  the  file  DEC92.DAT  and  put  in  the  file  BESIN.DAT  are,  as 
stated  in  Section  2.1.2.  not  subject  to  change.  If  the  program  runs  without 
difficulty,  only  the  final  output  data  need  to  be  interpreted;  intermediate 
output  data  can  be  ignored.  The  final  output  data  are  described  in  Section 
2.2.1. 

The  first  module  of  input  data  in  the  file  DEC92.DAT  is  listed  in  Section 
2.1.3.  The  second  module  of  input  data  in  the  file  DEC92.DAT  is  not  listed 
in  the  present  report:  this  module  is  exactly  the  same  as  the  data  listed 
under  the  heading  “Listing  of  the  second  module  of  the  sample  input  data” 
in  Section  2.1.3  of  |3].  The  complete  main  program  is  listed  in  Chapter 
7.  However,  the  subprograms  MODES.  BESIN,  BES,  INTERPOL,  PHI, 
DON,  FXY,  DECO.MP,  and  SOL\'E  are  not  listed  in  the  present  report 
because,  except  for  replacement  of  PH1(100).PH2(100),PH3(100),PH4(100) 
by  PHl(150),PH2(150).PH3(I50).PH4(150)  in  the  first  common  statement 
in  the  subroutine  PHI,  they  are  the  same  as  in  [3].  The  output  data  that 


^Modification  of  the  input  data  in  IN.DAT  may  requin  an  increase  in  the  storage  area 
allocated  to  some  arrays.  Minimum  allocations  of  arrays  are  given  in  Section  2.3. 


are  written  in  the  file  OUT.DAT  when  the  computer  program  is  run  with 
all  the  input  data  exactly  as  taken  from  the  file  DEC^.DAT  are  listed  in 
Section  2.2.2.  The  output  data  that  are  written  in  the  file  BESOUT.DAT 
are  not  listed  in  the  present  report  because  they  are  exactly  the  same  as  the 
output  data  described  under  the  heading  ‘‘Sample  output  written  in  the  file 
BESOUT  by  statements  in  the  subroutine  BESIN*  in  Section  2.2.2  of  [3].^ 


2.1  The  Input  Data 

There  are  two  modules  of  input  data:  the  input  data  in  the  file  IN.DAT  and 
the  input  data  in  the  file  B^IN.DAT. 


2.1.1  The  Input  Data  in  the  File  IN.DAT 

The  input  data  in  the  file  IN.DAT  are  read  from  there  by  means  of  two  groups 
of  statements. 


The  first  group  of  statements 

The  first  group  of  statements  that  read  input  data  from  the  file  IN.DAT  is 

REAO(20.10)  B.C.L1,L2.L3,BKM,XM,ZL1,ZL2 
10  FORMAT (4D14. 7) 

REA0(20.144}  KAM.BKA0.DBKA,KE3M,IfPHI,NZ 
144  F0RMAT(I4.2D14.7.3I4) 

READ(20, 146} (KE3(I) .Z-1 ,K£3N) 

146  F0RMAT(15I4) 

The  above  statements  occur  early  in  the  main  program.  The  data  read  in 
by  them  are  the  szime  as  the  data  in  Section  2.1.1  of  [3].  The  following 
description  of  these  data  is  based  on  the  description  in  Section  2.1.1  of  [3]. 

^There  is  no  file  named  BESOUT  in  [3].  The  file  BESOUT.DAT  is  meant.  There  are 
three  instances  in  Section  2.2.2  of  [3]  where  BESOUT  is  mistakenly  written  instead  of 
BESOUT.DAT. 


8 


In  the  first  of  the  nbove  three  tesd  sUtements, 


B  =  i 

a 

(2.1) 

C*  - 

a 

(2.2) 

Ll*^ 

a 

(2.3) 

II 

(2.4) 

4|o 

II 

a 

(2.5) 

The  fields  in  the  rectangular  waveguides  are  expanded  as  linear  onnlM* 
nations  of  all  rectangular  waveguide  modes  whose  cutc^  wavenuodbers  do 
not  exceed  BKM/6.  The  cutoff  wavenumber  ol  both  the  TM^  and  TE^ 
rectangular  waveguide  modes  is  ib^  given  by  (see  eq.  (A.8)  of  [1]) 


Hence,  the  values  of  p  and  q  are  limited  by 

<  BKM.  (2.7) 

The  field  in  the  circular  waveguide  is  expanded  as  a  linear  combination 
of  a  finite  number  of  circular  waveguide  modes.  This  number  is  controlled 
by  the  input  variable  XM.  The  even  and  odd  TMr*  and  the  even  and  odd 
TEr«  circular  waveguide  modes  are  taken  for  all  nonnegative  integer  values 
of  r  and  s  such  that  at  least  one  of  the  values  and  x’„  does  not  exceed 
XM.  Actually,  the  odd  TMr«  mode  and  the  odd  TEr.  mode  exist  only  for 
r  >  1.  Here,  is  the  root  of  the  Bessel  function  J,.,  and  x'„  is  the  s*^ 
root  of  Jf  where  J'  is  the  derivative  of  Jr  with  respect  to  its  argument.  If 
both  of  the  values  xoi  and  Xqi  exceed  XM,  then  execution  is  aborted  because 
XM  is  deemed  to  be  too  small.  The  cutoff  wavenumber  IcS*  of  the  even  and 
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odd  TM^  circular  waveguide  modes  and  the  cutoff  wavenumber  of  the 
even  and  odd  TEr*  circular  waveguide  modes  are  given  by 


uTti  _ 

”  “  a 

(2.8) 

II 

(2.9) 

Whereas  the  previously  described  first  seven  variables  in  the  first  read 
statement  are  real,  the  last  two  variables  in  this  read  statement  are  complex. 
These  complex  variables  are  ZLl  and  ZL2  given  by 

ZLl » 

(2.10) 

ZL2  -  z,y^ 

(2.11) 

where  Zi  and  are  the  impedance  loads  at  x  s  — £|  and  Li  vn  Fig.  2 

and  is  the  characteristic  admittance  of  the  TE^o  rectangular  waveguide 

mode.  Equations  (2.28)  and  (2.29)  of  [3]  are 

(2.12) 

(2.13) 

where  q  is  the  intrinsic  impedance  ^ven  by 

• 

(2.14) 

The  first  seven  variables  in  the  first  read  statement  are  obviously  dimen¬ 
sionless.  The  variables  ZLl  of  (2.13)  and  ZL2  of  (2.13)  are  dimensionless 
because  Zi/17  and  Z2/17  are  dimensionless.  Therefore,  all  the  variables  in  the 
first  read  statement  are  dimensiooless. 

The  variables  KAM,  BKAO.  and  DBKA  in  the  second  read  statement  are 
such  that  the  waveguide  mode  con^wrter  problem  is  solved  KAM  times,  once 
for  each  of  the  following  values  of  ka: 

ka  s  BKAO  +  (KA  -  1 )  •  DBKA,  KA  =  1, 2, 3,  • ,  KAM.  (2.15) 
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Here,  k  it  the  propagation  constant  and  a  is  the  radius  of  the  circular  wave¬ 
guide. 

The  variables  K£3M,  NPHI,  and  NZ  in  the  second  read  statement  and 
{K£3(I),  I  s  1, 2,  •  •  • ,  K£3M}  in  the  third  read  statement  ccmtrol  the  calcu¬ 
lation  and  the  writing  out  of  the  magnitudes  of  the  ^  and  r-components  of 
the  normalized  electric  field  in  the  apertures  Ai  and  A^.  These  magnitudes 
are 


i£i**>i 

I4"'l 

(2.16) 

i4"’l 

(2.17) 

where  is  the  root  mean  square  value  of  over 

the  cross  section  of  the  circular  waveguide  at  r  s  0.  Here,  is 

the  electric  field  of  the  z-traveling  waves  that  would  exist  in  the  circular 
waveguide  if  the  apertures  were  closed.  In  (2.16),  is  the  ^-component 
of  the  electric  field  at 

(*,!)=  +  +  =  2,.  -.NPHlI.o)  (2.18) 

in  aperture  i4i,  and  is  the  r-o>mponent  of  the  electric  field  at 

(*.»)=  +  (2.19) 

in  aperture  ^i.  In  (2.17),  is  the  ^component  of  the  electric  field  at 

(*.*)  •=({*.(-!  +  2yp~,|.,) .  J  =  1,2.  -  ,NPHl}  ,o)  (2.20) 

in  aperture  A2,  and  is  the  r-component  of  the  electric  field  at 

(♦.»)=  (0,{|(-l+2j^),J  =  l, 2, -.Nz})  (2.21) 

in  aperture  A2. 
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The  quantities  {K£3(I),  Is  1,2,  •  •  • ,  KE3M}  are  podtiue  integers  such 
that 


KE3(1)  <  KE3(2)  <  KE3(3)  •  •  •  <  KE3(KE3M  - 1)  <  KE3(KE3M)  (2.22) 

where  K£3(K£3M)  is  equal  to  either  KAM  or  KAM  + 1.  The  aperture  field 
magnitudes  (2.16)  and  (2.17)  are  calculated  and  printed  out  only  fiw  the 
values  of  Ira  of  (2.15)  for  whi^ 


KA  =  {KE3(1),  KE3(2),  KE3(3),  •  • ,  KE3(KAE)} 

where 

f  KE3M-1,  KE3(KE3M)*KAM  +  1 
KAE  s  \  _ 

I  KE3M,  KE3(KE3M)  «  KAM. 


(2.23) 

(2.24) 


The  second  group  of  statements 

The  second  group  of  statements  that  read  input  data  frcnn  the  file  INJ)AT 
18 

DO  48  KA«1.KAH 
IIEA0(20.109)  Rim.ItTEIl 
109  F0mUT(l0I4) 

liEAO(20 . 109)  (STNRCR)  .Ml  .Rnn) 

REA0(20 . 109)  (STERCR)  .Ml  .RTCR) 

DO  156  Ml  .Rim 
sHAX-srmCR) 

IF(SKAX.BI.O)  GO  TO  156 
Rl-R-1 

READ(20 . 164)  (BTMECS  .R)  .S-1  .SNAX) 

164  F0RMAT(4E14.7) 

IF(Rl.EQ.O)  60  TO 

REA0(20 .  IM)  (I1N0(S  .R)  .Ml.SMa) 

156  coirniiUE 

DO  167  R«l.RnR 
SMAX-STER(R) 

IF(SMAX.EQ.O)  CO  TO  167 
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Rl-R-1 

REAO(20.164)(BTEE(S.R)  .S-l.SMAX) 

IF(Rl.EQ.O)  60  TO  167 
READ(20.164)  (BTEO(S.R)  .S-l,SiUZ) 

167  covniniE 

48  CONTINUE 

The  index  KA  of  DO  loop  48  in  the  ebove  group  ci  FORTRAN  tUtements 
obtains  ka  of  (2.15).  The  first  statement  inside  DO  loop  48  in  the  above 
group  of  statements  occurs  at  about  540  lines  into  the  main  program. 

The  quantities  BTME(s,r  +  1),  BTMO(s,r  +  1)«  BTE^St**  +  !)• 
BTEO(s,  r+ 1)  indicate  th^  the  excitation  of  the  circular  waveguide  indndes 
the  z-traveling  waves  whose  dectric  fidds  are 

\  /  V  y/2^  / 

g-. 

where  £2?****  the  dectric  fidds  of  the  s- 

traveling  even  TM»,  wave,  the  ;*travding  odd  TMr«  wave,  the  s*traveling 
even  TEr«  wave,  and  the  ;>travrling  odd  TEr«  wave,  respectivdy.  These 
dectric  fidds  are  given  by  cqs.  (B.l).  (B.26),  (B.35),  and  (B.55)  of  [1].  FW- 
thermore,  0^  and  arc.  fcspectivdy.  the  wavenumber  and  the  char¬ 

acteristic  impedance  of  the  TM»«  mode.  Moreover,  0^  and  are,  re¬ 
spectivdy,  the  wavenumber  and  the  diaracteristic  admittance  of  the  T£r« 
mode. 

The  range  of  values  of  r  and  a  in  (2.25)  and  (2.26)  should  cover  all  the 
propagating  modes  in  the  circular  waveguide.  However,  he  who  does  not 
know  which  modes  propagate  in  the  circular  waveguide  will  not  know  what 
the  maximum  values  of  r  and  s  are  in  (2.25)  and  (2.26).  Hence,  he  will 
not  know  what  the  values  of  RTMR.  RTER.  {STMR(R),  R  s  1,  RTMR} 
and  {STER(R),  R  «  1.  RTER)  should  be.  Concerning  the  quantities  in  the 
previous  sentence,  the  Ibllowaig  action  is  taken  in  the  computer  program.  If 
one  or  more  of  these  quantities  is  too  small  so  that  the  the  range  of  values  of 
r  and  s  in  (2.25)  and  (2.26)  does  not  cover  all  the  propagating  modes,  then 
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the  coeffidents  of  the  miewng  propagating  inodes  are  set  equal  to  sero.  If 
one  or  more  of  the  quantitki  mentioned  in  the  third  firom  tl^  last  sentence 
is  too  large  so  that  nonpropagating  modes  are  included  in  (2.25)  and  (2.26), 
then  these  nonpn^Migating  modes  are  igmued. 

2.1.2  The  Input  Data  in  the  File  BES1N.DAT 

The  input  data  in  the  file  BESIN.DAT  are  read  from  there  by  of 

statements  in  the  subroutine  BESIN.  One  who  mereh  uses  the  oofiq>uter 
program  does  not  have  to  concern  cmeself  with  these  input  data  because  he 
will  never  have  to  change  any  of  them. 

2.1.3  The  Exemplary  Input  Data 

The  exemplary  input  data  consist  the  two  modules  of  input  data  in  the 
file  DEC92.DAT.  The  user  was  instructed  to  put  the  first  module  in  the  file 
IN.DAT  and  to  put  the  second  module  in  the  file  BESIN.DAT. 

The  first  module  of  input  data  is  for  the  structure  ci  Fig.  2  with 


i.u 

a 

(2.27) 

-  =  0.5 
a 

(2.28) 

~  =  1.838694 
a 

(2.29) 

MO 

(2.30) 

3^ 

It 

(2.31) 

ka  s  5.6. 

(2.32) 

The  value  of  Ls/a  in  (2.29)  was,  for  want  of  anything  better,  taken  fr<»n 
eq.  (2.49)  of  [3].  In  (2.30)  and  (2.31),  is  the  diaracteristic  admittance  of 
the  TEio  rectangular  waveguide  mode  so  that  Z\  and  Z2  are  matdied  loads. 
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Consequently,  the  electromagnetic  fidd  in  the  drcular  waveguide  will  not 
depend  <m  other  L\Ot  provided  that,  as  assumed  in  Ch^ter  1  [1],  L\ 

and  Ira  are  large  enough  so  that  any  evanescent  wave  emanating  frmn  either 
the  termination  at  x  s  — Li  or  that  at  x  »  £a  will  have  neglipbly 
amplitude  up<m  arrival  at  the  pertinent  iq)ertnre  in  the  drcular  waveguide. 

The  value  of  ha  in  (2.32)  was  obtained  by  setting 

KAM  »  1  (2.33) 


BKAO  «  5.6 

(2.34) 

DBKA  s  0.0. 

(2.35) 

The  values  of  the  variables  BKM,  XM,  NPHI,  NZ,  KE3M,  and  KE3(1)  in 
the  exemplary  input  data  are 

BKM  »  15.0 

(2.36) 

XM  «  40.0 

(2.37) 

NPHI »  81 

(2J8) 

NZ  =  21 

(2.39) 

KE3M  =  1 

(2.40) 

KE3(1)  =  1. 

(2.41) 

Because  of  the  values  of  KE3M  and  KE3(1)  given  by  (2.40)  and  (2.41),  the 
magnitudes  of  the  ^  and  2>components  of  the  normalized  electric  field  in  the 
apertures  Ai  and  A3  evaluated  at  ka  s  5.6  will  appear  in  the  output  data. 
For  ka  —  5.6, 

RrMR  =  2 

(2.42) 

RTER  =  3 

(2.43) 

STMR(l)  :=  2 
BTME(1,1)  =  1.0  +  ;0.0 
BTME(2,1)  =  0.9  +  i0.0  , 

'  (2.44) 

STMR(2)  =  1 
BTME(l,2)s:0.8+i0.1 
BTMO(l,2)  *0.7+j0.15  J 

1  (2«) 
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STER(1)  «  1 
BTEE(1,1)  *  0.6  +  i0.2 

STER(2)  *  2 
BTEE(1,2)  aO.l+iO.8 
BTEE(2,2)  a  0.15+i0.7 
BTE0(1,2)  a  0.2+i0.6 
BTEO(2,2)a0.5+i0.3  . 

STER(3)  a  1 
BTEE(l,3)a0.4  +  i0.4 
BTE0(1,3)  a  0.3+i0.5 


(2.46) 

« 


(2.47) 


(2.48) 


The  v«lues  of  the  B’s  in  (2.44)-(2.48)  were  arbitrarily  cKbaen. 

Only  the  first  nodule  of  input  datn  in  the  file  DEC92.D AT  is  listed  bdosr. 
The  second  module  of  input  daU  in  the  file  DEC92.DAT  is  exactly  the  same 
as  listed  in  Section  2.1.3  of  [3]. 


Listing  of  the  first  module  of  input  data  in  the  file  DEC02.DAT 


O.llOOOOOIKOl  0.5000000P»00  0.40000000^2  0.40000000^02 
0.1838604PH)1  0.i500000D«02  0.40000000«02  0.2000000PH>1 

o.ooooooon»oo  o.soooooop»oi  o.oooooooo^ 

1  0.560000<X)-»01  0.00000000-»00  2  02  22 

2 

2  3 

2  2 
2  2  2 

0.2000000B«02  O.OOOOOOOE'MM  O.OOOOOOOE^OO  O.0OOOO0OE*OO 
O.SOOOOOOE-»00  0.2000000E'»00 
0.7000000E-»00  0.2SOOOOOE-»00 
0.«000000E«00  0.2000000E^00 

0.2000000E«00  O.OOOOOOOE^OO  0.2E00000E*00  0.7000000E«00 
0.2000000E<HX>  O.OOOOOOOE^OO  O.SOOOO00E«OO  0.3000000E400 
0.4000000E^00  0.4000000E<»00 
0.3000000E<»00  0.6000000E«00 
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2.2  The  Output  Data 

The  output  d«U  consist  of  the  input  <Uta,  aotae  intermediate  output  data, 
and  the  final  output  data.  The  input  data  were  described  in  Sections  2.1.1 
and  2.1.2.  The  final  output  data  are  described  in  Section  2.2.1.  The  inter¬ 
mediate  output  data  are  included  in  the  exemplary  output  data  listed  and 
discussed  in  Section  2.2.2. 


2.2.1  Description  of  the  Final  Output  Data 

The  final  output  data  consist  of  the  normalized  coeflSdents  of  the  TEio  waves 
in  the  rectangular  waveguides,  the  normalized  coefficients  of  all  the  waves 
of  all  the  propagating  modes  in  the  circular  waveguide,  the  magnitudes  of 
the  4  end  x  components  of  the  normalized  electric  field  in  the  apertures, 
the  electrical  length  ka,  the  time-average  incident  power,  the  time-average 
transmitted  power,  and  the  time-aver^se  reflected  power. 

The  electric  field  of  the  TEjo  waves  in  the  left-hand  rectangular  wav^uide 
in  Fig.  2  is  written  as 


(ClOUT)£5=-e^'^*»*  +  (ClIN)£^e-'ft*** 


(2.49) 


Here,  is  the  total  electric  field  in  the  waveguide  and  the  subscript  “prop** 
denotes  the  part  that  consists  of  all  the  waves  of  the  forward  and  backward 
traveling  TEio  modes  that  propagate  in  the  waveguide.  In  (2.49), 
the  electric  field  of  the  — z-traveling  TEjo  mode,  and  u  the  electric 
field  of  the  r-traveling  TEio  mode.  These  fields  are  given  by  eqs.  (A.  15) 
and  (A.14)  of  [1].  The  quantities  0xo  and  in  (2.49)  are,  respectively, 
the  wavenumber  and  the  characteristic  admittance  of  the  TEio  rectangular 
waveguide  modes.  Also, 


a  sin 


(2.50) 


The  quantities  CIOUT  and  CUN  in  (2.49)  are  the  normalized  coefficients  of 
the  TEio  waves  in  the  left-hand  rectangular  waveguide.  The  electric  field  of 
the  TEio  waves  in  the  right-hand  rectangular  waveguide  in  Fig.  2  is  [£^^^]pt«p 


written  M 


The  quantities  C20UT  and  C2IN  in  (2.51)  are  the  nocmalised  coefficients  of 
the  TEio  waves  in  the  right-hand  rectangular  waveguide. 

The  coefficients  ClOUT,  ClIN,  C20UT,  and  C2IN  in  (2.49)  and  (2.51) 
are  variables  in  the  main  program.  They  axe,  for  ha  s  BKA,  written  out  by 
means  of  the  statements 

U1111E(21,69}  CIOUT.CIIM 

69  F0RIUTC*C10trr-\2E14.7.\  Cini*\2E14.7) 
imiTE(21,70)  C20VT.C2IN 

70  F0RH4T(»C200T«*.2E14.7.*.  C2II-\2E14.7) 

Computer  program  variables  ClOUTS,  ClINS,  C20UTS,  C2INS,  and  PT 
are  defined  by 

C100TS-C100T*aW  JO  (ClOUT) 

ClIliS»ClIN*CONJG  (ClIN) 

C20UTS-C20UT*C0ir  JO  (C20UT} 

C2I1IS-C2IN*C0NJG  (C2IH) 

PT-C10UTS-C1IMS+C20UTS-C2INS 

and  printed  out  according  to 

VRITE(21,82)  ClOUTS, ClDrS.C20UTS,C2IMS,Fr 
82  FORMATC* ClOUTS- ’,E14. 7, \  C1IIS-* ,E14.7, » ,  C20UTS-*. 
1E14.7/»C2IHS-'  ,E14.7,  * ,  l»T-»,E14.7) 

The  previous  three  groups  of  FORTRAN  statements  are  inside  DO  loop  48 
whose  index  KA  obtains  ka  of  (2.15).  The  quantities  ClOUTS,  ClINS, 
CLOUTS,  and  C2INS  are,  respectively,  the  time-average  powers  of  the  out¬ 
going  TEto  wave  in  the  left-hand  rectangular  waveguide,  the  incoming  TEio 
wave  in  the  left-hand  rectangular  waveguide,  the  outgoing  TEio  wave  in  the 
right-hand  rectangular  waveguide,  and  the  incoming  TEio  wave  in  the  left- 
hand  rectangular  waveguide.  Here,  ’‘outgoing”  means  going  away  firmn  the 
pertinent  aperture,  and  “incoming”  means  coming  in  toward  the  pertinent 


(C20UT)£^e^^  -t-  (C2IN)£g“e->^ 


(2.51) 


aperture.  The  quantity  PT  is  the  tune-average  power  Pt  transmitted  into 
the  rectangular  waveguides.  Here,  Pt  is  given  by 

+  (2.52) 

where  pf^^  and  P**'  will  be  pven  by  (4.7)  and  (4.8),  respectively. 

The  electric  field  of  all  the  waves  o{  all  the  propagating  inodes  in  the 
circular  waveguide  is  written  as 

t^”U= 

ex: 

rmO  tml 

+EE 

rmlasl 

+EE 

raO«ail 

+EE 

r«l  «sl 

Here,  is  the  total  electric  field  in  the  circular  waveguide,  and  the  sub¬ 
script  ‘‘prop”  denotes  the  part  that  consists  of  all  the  waves  of  all  the  prop¬ 
agating  modes  in  the  circular  waveguide.  In  (2.53),  is  the  electric 

field  of  the  z-traveling  even  TMr«  circular  waveguide  mode,  is  that 

of  the  —z-traveling  even  TM»,  circular  waveguide  mode,  is  that  of 

the  z-traveling  odd  TM„  circular  waveguide  mode,  is  that  of  the  — z- 

traveling  odd  TMr*  circular  waveguide  mode,  is  that  of  the  z-traveling 
even  TEr,  circular  waveguide  mode.  £^"  is  that  of  the  -z-traveling  even 
T£r«  circular  waveguide  mode.  is  that  of  the  z-traveling  odd  T£r« 

circular  waveguide  mode,  and  is  that  of  the  -z-traveling  odd  T£r« 

circular  waveguide  mode.  These  fields  are  given  by  eqs.  (B.l),  (B.2),  (B.26), 
(B.27),  (B.35),  (B.36),  (B.55).  and  (B.56)  of  [1].  The  quantities  0^  and 
in  (2.53)  are,  respectivrly.  the  wavenumber  and  the  characteristic 


(BTME(s,r  -h  -  (CTME)£y*“e“^^^ 


(BTMO(s,r  -H  -  (CTMO)£y*~e“;^^ 

(BTEE(s,r  +  l))£l^e^^”^  +  (CTEE)£^“e-'^^ 

■V'"— ■■  IW— i— — — 1— i^— A—— 111  III  W—iW— Jaij  lffww.1  i» 

(BTEO(s.r  l))£j;^e>^”fa  -f  (CTEO)£^e“><>y^ 

/vTEw 

V 


(2.53) 
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impedance  of  the  TM,«  circular  waveguide  modes.  FHirthermore,  0^  and 
are,  respectively,  the  wavenumber  and  the  characteristic  admittance 
of  the  TEra  circular  waveguide  modes.  The  double  summations  with  re> 
spect  to  r  and  s  in  (2.53)  are  over  all  the  pnqxagating  modes  of  the  circular 
waveguide.  The  quantities  BTME(s,r+l),  CTME,  BTMO(s,r+l),  CTMO, 
BTEE(s,r+l),  CTEE,  BTE0(8,r+l),  and  CTEO  in  (2.53)  are  the  normal- 
ized  coefficients  of  all  the  waves  of  all  the  propagating  modes  in  the  circular 
waveguide.  EUtch  of  these  coefficients  depends  cm  r  and  a.  The  B*s  were  read 
in;  CTME,  CTMO,  CTEE,  and  CTEO  will  be  given  by  (7.149),  (7.155), 
(7.189),  and  (7.195),  respectively. 

The  normalised  coefficients  CTME,  CTMO,  CTEE,  and  CTEO  are  vari¬ 
ables  in  the  main  program.  They  and  their  magnitudes  squared 

CTNES-CTME*CON JG (CTME) 

CTM0S-CTM0*C0N JG (CTMO) 

CTEES-CTEE*COir  JG  (CTEE) 

CTE0S-CTE0*C0MJG (CTEO) 

are  written  out  by  the  means  of  the  statements 

DO  172  R«1,RTM 
SMAX*STM(R) 

IF(SMAX.Eq.O)  GO  TO  172 
Rl-R-1 

DO  173  S«1,SMAX 
WRITE(21,160)  R,S 
160  F0RMAT(’R-M4,',  S-M4) 

VRITE(21.162)  CTME.CTMES 
162  FORMAT (' CTME- ',2E14. 7,’,  CTMES-' ,E14.7) 

IF(Rl.EQ.O)  GO  TO  173 
WRXTE(21.158)  CTMO.CTMOS 
158  F0RMAT('CTM0-',2E14.7,',  CTMOS-’ ,E14.7) 

173  CONTINUE 
172  CONTINUE 

DO  180  R-l.RTE 
SMAX-STE(R) 

IF(SMAX.EQ.O)  GO  TO  180 
Rl-R-1 
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DO  181  S-l.SHAZ 
HRnE(21.160)  R.S 
VRnE(21.16S)  CTEE.CTEES 
165  F0RBUT('CrEE-*.2E14.7.».  CTEES— .E14.7) 
IF(Rl.EQ.O)  GO  TO  181 
imnE(21,90)  CTEO.CTEOS 
90  F0RIttT(’CTE0«*,2E14.7,‘.  CTEDS-* ,E14.7) 

181  comiruE 

180  CONTIHUE 


The  above  group  of  FORTRAN  statements  is  inside  DO  loop  48  whose  index 
KA  obtains  ka  of  (2.15).  The  quantities  CTMES,  CTMOS,  CTEES,  and 
CTEOS  are  the  time-average  powers  of  the  waves  of  the  — s-traveling  even 
and  odd  TM,,  and  even  and  odd  TEr,  propagating  modes  of  the  circular 
waveguide.  For  example,  CTMES  is  the  time-average  power  of  the  wave  of 
the  — x-traveling  even  TM,«  mode  in  the  circular  waveguide. 

Computer  program  variables  E3A1PS(J),  E3A1ZS(J),  E3A2PS(J),  and 
E3A2ZS(J)  are  defined  by 


E3A1PS(J)  = 

E3A1ZS(J)  = 

E3A2PS(J)  = 

E3A2ZS(J)  = 


^1) 


i  z*0 

|£M«1 

’  ^  X 

/  *-  (  *  + 

l.  =  o 

f  e  =  0 

(2.54) 

(2.55) 

(2.56) 

(2.57) 


where  |£^***(,I'^,Q)|ni»  ihc  root  mean  square  value  of  the  magnitude  of 
the  transverse  electric  field  of  the  r-traveling  waves  that  would  exist  over  the 
cross  section  of  the  circular  waveguide  at  r  =  0  if  the  apertures  were  closed. 


In  (2.54)  and  (2.55),  and  ate,  lespectivelijr,  the  magnitudes  of 

the  ^  and  z-components  of  the  electric  field  at  the  indicated  values  of  ^ 
and  z  in  the  left-hand  aperture  in  F^.  2.  In  (2.56)  and  (2.57),  and 

are,  respectivdy,  the  magnitudes  of  the  ^  and  z-cmnpmients  oi  the 
electric  field  at  the  indicated  values  of  ^  and  z  in  the  right-hand  ^>erture  in 
Fig.  2.  The  right-hand  aides  of  (2.54)-(2.57)  will  be  the  magnitudes  of  the 
right-hand  sides  of  (6.6),  (6.8),  (6.7),  and  (6.9),  respectively. 

The  dimensioned  variables  ^AlPS,  E3A1ZS,  E3A2PS,  and  E3A2ZS  in 
the  main  program  are  written  out  the  statements 

VRnE(21 . 140)  (E3A1PS(  J) ,  >1  ,HPBI) 

140  F0RMAT(*E3A1PSV(5E14.7» 

VRnE(21 .141)  (E3A1ZS(J) .  J-l.MZ) 

141  F0RMAT(*E3A1ZS*/(SE14.7)) 

VRITE(21 . 142)  (E3A2PS(  J) .  >1  .VPHI) 

142  F0RMAT(*E3A2PSV(SE14.7)) 

WRITE(21 .143) (E3A2ZS(J) . J-1 .HZ) 

143  F0RMAT(»E3A2ZS»/(SE14.7)) 

The  above  group  of  statements  is  inside  DO  loop  48  whose  index  KA  obtains 
ka  of  (2.15).  However,  this  group  of  statements  is  executed  only  for  those 
values  of  KA  on  the  right-hand  side  of  (2.23). 

The  computer  program  variables  BKAPLT(KA),  PINC(KA),  PTRAN 
(KA),  and  PREFL(KA)  are  set  equiJ  to  ka,  the  time-average  power  of  the 
z-traveling  waves  that  would  exist  in  the  circular  waveguide  if  the  aper- 
tures  were  shorted,  the  time-average  power  transmitted  into  the  rectangular 
waveguides,  and  the  time-average  power  reflected  in  the  circular  waveguide, 
respectivdy.  Here,  ka  is  given  in  terms  of  KA  by  (2.15),  and  all  these  time- 
average  powers  are  evaluated  at  this  value  of  ka.  The  computer  program 
variables  cited  in  the  first  sentence  of  this  paragraph  are  vmtten  out  by  the 
statements 

VRITE(21 .149)  (BKAPLTCl)  .I-l.KAM) 

149  F0RNATC'BICAPLT'/(SE14.7)) 

URITE(21 .189) CPINC(I) .1*1 .KAM) 

189  F0RMAT('PIKC»/(SE14.7)) 

«RITE(21 .92)  CPTIUN(I)  .1-1  .KAM) 

92  F0RMAT(*PTRAirV(5El4.7)) 
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IIRnE(21.139)  (FREFL(I)  .I-l.KAN) 

139  F0RIUT(»FREFLV(5E14.7)) 

AU  three  of  the  variables  PINC(I),  PTRAN(I),  and  PREFL(I)  are  calculated 
independently  so  that  satisfaction  of  the  equation  of  conservation  time- 
average  power 

PINC(I)  *  PTRAN(I)  +  PREFL(I)  (2.58) 

gives  some  assurance  that  the  values  of  these  three  variables  ate  accurate. 


2.2.2  The  Exemplary  Output  Data 

The  output  data  that  resulted  when  the  computer  program  was  run  with  the 
exemplary  input  data  of  Section  2.1.3  are  presented  in  here  in  Section  2.2.2. 


Listing  of  the  output  data  written  in  the  file  OUT.DAT 

B .  C  .U  .L2  .U  .BXR.XH  .ZLl  .ZU 

O.llOOOOW^l  O.SOOOOOOO^OO  0.400000(»403  0.4000000IH02 
0.1S386B4D401  O.IBOOOOOO^OS  0.4000000IK03  0.100000(»401 
O.OOOOOOOD400  O.IOOOOOOIKOI  O.OOOOOOOIHOO 
KAN>  1,  BKAO>  0.6S00000IK01.  OBKA*  O.OOOOOO0D4OO 
KE3N«  1.  IPBX«  81.  IZ«  31 
KE3 

1 

IIIAZ«  3 

m 

6  S  3 

KTM«  6.  BTE«  11.  Kl>  16 

BXB«  0.6160000E401 

TUC 

0.1388117E*01  O.OOOOOOOE^OO  0.8773033S«00  0.0000000B4O0 
0.630301SE-»00  0.0000Q00E-»00  0.4TE6368E«00  O.OOOOOOOE4-00 
0.4834891E<i>00  0.00000Q0E-»00  O.OO0OO00E<H)O-O.86O1761B-i-0O 
-0.30098S4E-»00  0.0000000E<»0O-0.11S7968B«01  0.0000000B’K>0 
•0.17780e3E01  0.000O00OB«OO-O.S088OO7B<»00  O.OOOOOOOB^OO 
-0.7304006B<K)0  0.0OO0000B'»0O-0.1139867E«01  O.OOOOOOOB-i'OO 
-0.iei3379B4Ol  0.000000OB-K>O-O.3103489E-»01  O.OOOOOOOB-^00 
~0.70098S9E*01  0.0000000E«00-0.307356eB4Ol  O.OOOOOOOB-fOO 
0.1388117E-»01  O.OOOOOOOE-»00  0.8773033E«00  O.OOOOOOOE-i-00 
0.6303015E-»00  0.0000000E«00  0.47863S8E<»00  O.OOOOOOOE^OO 
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0.4t24Wll400  0.0000Q00I400  O.OOOOOOOI^OO-O.aSOlTSlB^OO 
•0.20098S4I400  O.OOOOOOOI^OO-O.llSTMM^Ol  O.OOOOOOM^OO 
-0.17TtQ«»«01  0.0000000l-»00>0.i0aa007l«00  0.0000000I400 
-0.730400SI400  O.OOOOOOOI-»00-0.113Mi7|«OS  0.0000000l«00 
-0.1613SrM<H)l  0.0000000l-»00>0.210a48w«01  O.OOOOOOM^OO 

-o.2oonBi»«oi  o.oooooooi^ooHi.aoTassn^oi  o.oooooooi^oo 
Km>  2.  im-  3 
sm 

2  1 

snt 

12  1 

■me  i.s) 

0.1000000D401  O.OOOOOOOD^  O.0OOOOOOO«OO  O.OOOOOOOfrHM 
BTHB(  2.S) 

O.3OOOO0MKO0  O.IOOOOOOIKOO 

■me  2,s) 

0.7000000D400  0.1M00000400 
tTElC  I.S) 

0.60000000400  0.2000000l>«00 
BTBCC  2.S) 

0.1000000D400  0.8000000D400  0.1SO0Oa»«OO  0.70000000400 
tTE0(  2.S) 

0.200000Q04QO  0.00000000400  O.tOOOOOQD^OO  o.sooooocsy^oo 
BTEE(  3.S) 

0.4000000D400  0.40000000*00 
BTEOC  3.S) 

0.30000000*00  0.10000000*00 
KTN«  3.  ITS-  1 
STM 

2  11 
STB 

12  111 
T1 

0 . 35 18446E*00-0 .34f37«St*00-0 . OOSOOlOt-Ol-O . M93077E-02 

0 . 4838421E-01*0 .4TMl3TE>«t  •« .  StOOtTlB-Ot-O .  S33926SE-02 

0.424907SE-0l>0.3330TiSB*0O-«.t42e0i0t*0O  0. 1330089E*01 

0.1630048E-01  0.l90299flE»00  0  99018211-01  0.41096898*00 

0.11890188-01  0.7SillllB-«l-«.S9T9t88««0  0.49830838*00 
-0 . 33343118-01-0 .4703018-01 •«  94498SK-01-0 . 19886138-01 

0 . 48466118-01-0.11128818-01-8 .37833778-01-0 .37108128-01 
-0.36368998*00  0.37681678-01  0.17836888-01-0.38014868-01 

0.75642698-01  0.63437068*80  8.78466968-01-0.17820458*00 
-0.40013198-02  0.13880108*00  8. 33986788-01-0.64304678-01 
-0.1800006E*00-0.3U9993S-01  8.87446398*00  0.33780308*00 
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0.23S2414l«0(H>.49e74MBH)l  0.13916001400  0.19336091400 
0.1064943I400H). 10063031-01  0. 17606416400-0 .4616SaOB400 
-0.10000016400  0.37016476400  0.20040746-01-0.11343336400 
-0.16620626-01  0.10673036400  0.26034646-01-0.01766426-01 
-0.63374146400  0.11000766400  0.14100766400-0.14604306400 
T 

0 . 16133666400-0 .60622066-01-0 . 30400366-01-0 . 60300016-03 
0.70300066-01-0.66463006-01-0.976336a-01  0.37040306-01 
-0.10637066400-0.10060476400-0.03200776400-0.20001666400 
-0.32072306400-0.67301066-01  0.600670n-01-0. 17431406400 
-0.30207066-01-0.10323016-01  0.26227166400-0.63270626400 
-0.96171646-01  0.10300606400  0.41333046-01-0.30163666-01 
-0.23337666-01  0.00430176-02  0.14104046-01-0.11330106-01 
0.11606066400-0.26000306-01-0.44100136-01  0.36060006-01 
0.61036106-01  0.30419416400  0.23041416-01-0.30643136-01 
0.66003626-02  0.16207036400  0.10000746-01-0.71246346-01 
-0.26130616400  0.22136306-02-0.14074006400  0.36092006400 
-0.40003466400  0.26120406400-0.60716636-01-0.67667066-01 
-0.41675036-01  0.23003606-01-0.61063206-01  0.61060366400 
0.6067063B-01-0. 33003096400-0. 26003616-01  0.16400326400 
0 . 7602720E-02-0 . 07 166026-01-0 . 76606176-02  0 . 40330076-01 
0 . 16460466400-0.60063276-01-0 . 62970806-01  0 . 42607776-01 
C100T«-0.76246976400-0. 37006606400.  CltB»  0.00000006400  0.00000006400 
C20OT— 0.13807736400  0.33464566400.  C3ZI>  0.00000006400  0.00000006400 
CIOOTS*  0.6686664E«00,  CIZIS*  0.00000006400.  C300TS>  0.13137466400 
C2ZIS«  0.00000006400.  PT"  0.70004016400 
%•  1.  Sa  1 

rniE«  0.26611256400-0.10600176-01.  mCES«  0.70930296-01 
CTMEa-O. 73388766400-0. 10600176-01.  CTKESa  0.638701«400 
E«  1,  S*  2 

rni6*  0.12723686400-0.00110446-03.  XTII6S«  0.16386476-01 
CTMEa-O. 77276336400-0. 98110446-03.  CTII6S*  0.60725026400 
6«  2.  5-  1 

ZTH6«  0.12984506400  0.60663476-01.  XTIIXS>  0.31712586-01 
CTH6>-0. 67015506400-0. 30337636-01.  CTICESa  0.46002826400 
ZTMO>  0.40367746-01-0.63906696-03.  X71I0S>  0.166024<»-02 
Cni0a-0.66963336«00-0. 16630076400.  CTOOSa  0.45964436400 
Rs  3.  S>  1 

ZTII6>  0.43987316400  0.76150756-03  .  671I6S>  0.10354396400 
Cni6«  0.43087216400  0.75150756-03.  C71IXS>  0.10354396400 
ZTMOa-O. 24833886-01  0.7602&ii6-0t.  ZTlIOSa  0.63664646-03 
0760—0.34832886-01  0.7503w7»i;  -0l.  CTMOS>  0.63664646-02 
!•  1,  S«  1 

6766—0.78697206-01  0.20176186400.  X7E6S«  0.91312386-01 
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CTII--0.67M972I400  0.91TSlt3IH)l.  CTBS* 
I*  2.  S«  1 

nil*H>.M69223l-01  0.30932881-01.  HIIS* 
CTII^.10S9932l400H).789397ei*00.  CRB* 

nn*  0.1700T1R+00  o.ioooorii^oo.  xrqs* 

CRO»-0.290383<«-01-0.4679980I«00.  CRB* 
I*  2.  3*  2 

XRI»-0.649«7i98-01  0.230132SI-01.  XRB* 
CRI*H).214M73l400-0.imi>R^.  CRB* 
XR0>  0.39093411-01  0.34912141-01.  XROS* 
CRO— 0.4910399l«00-0.29S19rtl«00.  CROS* 
3*  3.  3*  1 

XRI*-0.7980280B-01  0.2899089t«00.  XRIS* 
CRB*-0.479S029B400-0.1113133B«00.  CRB* 
XRO— 0.1447939l<»00  0.8917S9R-01.  XROS* 
CR0*-0.4447089S40O-0.4138341l«00.  CROS* 
1*  4.  S*  1 

XRB*  0.3949763X400-0.097837R-01.  XRB* 
CRB-  0.29497S3B<»00-0.9979S7R-01.  CRBS* 
XRO— 0.3324401B400-0.164S972B-01.  XROS* 
CRO— 0.2S24401B«00-0.1S43«72B-01.  CROS* 
B*  S.  3*  1 

XRB— 0.241943SB-01  0.8780342B-01.  XRBS- 
CRB— 0.241943SB-01  0.6760343E-01.  CRBS* 
XRO— 0.3080200B-»00  0.3961476B-01.  XROS* 
CRO— 0.30S0299B400  0.3961479B-01.  CROS* 
?XI*  0.5036000E-»01  PIN*  0.7959231B400  PB* 
PTQTAL*  0.5935001E«01 
PR*  0.7899401E4>00.  PBIU*  0.7659330E«00 
PBX2 

O.OOOOOOOE^OO 


0.1993496E-i>00 

0.3929991E*00 

0.5890489E-i>00 

0.78639e2E400 

0.981747R400 

0.11780B7E-f01 

0.137444R^01 

0.1670799Et01 

0.1797149E-K)1 

0.1993498E'»01 

0.2159e4SE«01 

0.23S81B4E*01 

0.2S52544E«01 


0.3926991E-01 

0.23891968«00 

0.4319990E^00 

0.6383185E-»00 

0.824e680E'»00 

0.1031018E<»01 

0.121736R-»01 

0.1413717E'»01 

0.1610066E<»01 

0.1809416E-K)1 

0.2003798E401 

0.21991188*01 

0.239S494E«01 

0.2S91814E«01 


0.7883982E-01 

0.3748883E<M>0 

0.47123898400 

0.6975884E<H)0 

0.8939380E400 

0.106028R<H>1 

0.138683R4'01 

0.14S298R«01 

0.1840338E«01 

0.1846980E401 

0.20420368401 

0.22383BX401 

0.24347348401 

0.29310848401 


0.4990484B400 

0.1009991B-01 

0.9303103B400 

0.499B71SI-01 

0.21999018400 

0.498973B-02 

0.S0S92SB400 

0.299071B-02 

0.2937106B400 

0.99192731-01 

0.23944S7B400 

0.20391SB-01 

0.39909200400 

0.7942011B-01 

0.7943011B-01 

0.S42997B-01 

0.94299780-01 

0.S14093B-02 
0. 61409298-02 
0.438099B-01 
0.43699080-01 
0.6146091B401 


0.117909R400 

0.31416038400 

0.6106088B400 

0.70086838400 

0.90320798400 

0.109966R401 

0.129600R401 

0.14B226R401 

0.19899088401 

0.18849608401 

0.20813068401 

0.22778668401 

0.247400tt401 

0.267036tt401 


0.16707998400 

0.36342928400 

0.649778R400 

0.74912838400 

0.94247788400 

0.113982R401 

0.133617R401 

0.16316298401 

0.17278798401 

0.19242268401 

0.21206768401 

0.23199268401 

0.26132748401 

0.27009248401 


0 
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e.a74MMi«oi 

o.aMutMfOi 

0.9l41ini««l 

z 

0.000000W«40 

o.7tssMai«oe 

0.UT0TMi«01 

0.3SM1MI«01 

0.S14iiMI«01 

I9A1M 

0.1440Slil*01 

0.1S7TMn«01 

o.iaisi7ii«oi 

0.10161831401 

0.86314131400 

0.7646am400 

0.740SS10I400 

0.77086141400 

0.81806668400 

0.63141388400 

0.76880618400 

0.71676768400 

0.64678738400 

0.86613738400 

0.86887048400 

0.88314708400 

0.68038388400 

6381ZS 

0.13738648401 

0.10388618401 

0.03608378400 

0.13308818401 

0.1836633E401 

E3i3PS 

0.66038748*01 

0.84680128*01 

0.61144368*01 

0.73011478*01 

0.180892U400 

0.31606088400 

0.44637068400 

0.61631808400 

0.51366798400 

0.46670348400 


0.37881638401 

0.a66ttl38401 


0.18707668400 

0.64347788400 

0.17278768401 

0.38133748401 


0.14377188401 

0.13814828401 

0.11786778401 

0.67680748400 

0.82801028400 

0.78636068400 

0.78313968400 

0.78783868400 

0.8164669(400 

0.81784038400 

0.77400868400 

0.70830838400 

0.63736078400 

0.88831048400 

0.86191608400 

0.88206418400 


0.13600368401 

0.98334368400 

0.98397798400 

0.13116068401 


0.98487938*01 

0.79061918*01 

0.46343078*01 

0.90333188*01 

0.30733208400 

0.34680608400 

0.46813028400 

0.83181068400 

0.80894948400 

0.48641668400 


0.38374331401 

0.30337838401 


0.31418638400 

0.10668878401 

0.18646888401 

0.a67038tt401 


0.1tt666W401 

0.13218318401 

0.11386468401 

0.6440tfl«400 

0.80703188400 

0.780600n400 

0.78824788400 

0.76837118400 

0.83346448400 

0.81161668400 

0.76146388400 

0.69066668400 

0.63863348400 

0.88134018400 

0.88891608400 

0.88137088400 


0.12234808401 

0.64163188400 

0.90614828400 

0.13668868401 


0.66733888*01 

0.73396848*01 

0.48293158*01 

0.11004208400 

0.33468198400 

0.37367388400 

0.48338868400 

0.83368688400 

0.46670748400 

0.44881718400 


0.36867088401 

0.30630838401 


0.47133888400 

0.13866378401 

0.30430388401 

0.38274338401 


0.14173338401 

0.13888838401 

0.1068709(401 

0.61068068400 

0.76023718400 

0.74814388400 

0.76488388400 

0.80377168400 

0.83383778400 

0.804626M400 

0.74806668400 

0.67668688400 

0.61488968400 

0.87817818400 

0.86680888400 

0.56074838400 


0.11688608401 

0.61630068400 

0.10864728401 

0.14689868401* 


0.93861838*01 

0.66193088*01 

0.49519218*01 

0.13197068400 

0.36365938400 

0.36030878400 

0.49664568400 

0.52390608400 

0.48661188400 

0.44360618400 


0.39089738401 

0.31023338401 


0.63831888400 

0.14137178401 

0.21961188401 

0.29648138401 


0.13966888401 

0.13829138401 

0.10888408401 

0.88080608400 

0.77879318400 

0.74777468400 

0.77180068400 

0.80643388400 

0.82380778400 

0.76877338400 

0.73413638400 

0.66263178400 

0.60501968400 

0.56697358400 

0.55460938400 

0.55044978400 


0.11037488401 

0.91486488400 

0.11328108401 

0.15109948401 


0.89915438*01 

0.57889748*01 

0.59063838*01 

0.18868388400 

0.29061588400 

0.43333388400 

0.50803208400 

0.51948188400 

0.47596088400 

0.44148878400 
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0.443141»«00  0.4«tllT7l«00  0.4M3S4TB400  0.4TtrS4tl«00  0.4ttl7Stl«00 
0.i27l062l«00  0.SS2UtSI«00  0.i0441fn«40  0.IMimt»00  O.TOCTtSOt^OO 
o.7t7iasii«oo  0.t3S0Siii400  o.to4oi»i«oo  o.tmMii^oo  o.ioMsm^oi 
0.U3H8««01  0.ia»3MI«01  0.ia07St2l«01  0.1W14att«01  0.1473MSB401 
o.isiatmtoi  o.taaoiaM^oi  o.iroiT43t*4i  o.iTfrMti«oi  o.nawiii*oi 
0.1tTtttll«01  O.ltai0S7l«Ol  0.1Miani««l  0.iWI>W*01 

O.MOIMU^OI 

0.Mtia4n«00  0.t73C0MI«00  0.t3»7tfl«00  e.77«Ti7l«00  0.nTt044i*00 
0.«7a4SSSI«00  O.MOaam^OO  0.44M41SI^  0.4S«ST4M«00  O.Mia0MB»OO 
0.eS3tatll*00  0.SfM2SMI«O0  0.fatti0SI«O0  0.4414Tiai««0  0.M04tm*00 
0.4SWM1I400  0.S3704tll«00  0.0474171400  0.7003141400  0.7030231400 
0.80403661400 
BUFLT 

0.66000001401 

PXfC 

0.80360001401 

mil 

0.70004011400 

PUFL 

0.61460611401 

Diocuoiion  of  the  above  output  data 

First,  the  input  data  read  in  by  the  first  group  of  statements  listed  in  Section 
2.1.1  are  written  out.  Then,  the  variable  NMAX  and  the  array  MM  are 
written  out  after  being  calculated  in  the  subroutine  MODES  (see  Chi^iter  4 
of  [3]).  Next,  the  nundber  KTM  of  waveguide  modes  used  in  eadi  rectangular 
waveguide,  the  number  KTE  of  TE  waveguide  noodes  used  in  eadi  rectangular 
waveguide,  K1  defined  by 


K1  =  KTM  +  KTE 

and  BKB  defined  by 


(2.59) 


BKB  s  kb  (2.60) 

are  written  out.  The  elements  >  s  1,2,  •  •  • ,  A’2},  which  will  appear 

in  (4.12)  and  (4.17),  are  written  under  the  heading  YREC.  Here,  is  the 
it*^  element  of  Y'  +  Y^  where  F*  and  Y^  are  the  rectangular  waveguide 
admittance  matrices  given  by  eqs.  (2.25)  and  (2.27)  of  [1]. 
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Next,  the  iaput  del*  reed  in  by  the  secoml  group  ol  etetementt  hated 
in  Section  2.1.1  are  written  out.  Then,  RTM  of  (7.29),  RTE  ci  (7.47), 
{STM(R),  R  »  1,2, ••  • , RTM)  of  (7.30)  and  {STE(R),  R  »  1,2,  • , RTE} 
^  (7.48)  axe  written  out.  The  values  td  RTM,  RTE,  STM,  and  STE  are,  re- 
spectiudy,  those  RTMR,  RTER,  STMR,  and  STER  for  which  the  ranges 
cf  values  of  r  and  s  in  (2.2S)  and  (2.26)  eo>uer  all  the  propagating  but  no 
evanescent  modes  in  the  circular  waveguide. 

The  normalised  elements  of  the  excitation  vector  given  by 

(3.13)  are  written  out  under  the  heading  TI.  The  dements  {iT^}  are,  as  can 
be  gleaned  £n»n  Section  7.4.4,  ordered  as 

{7?™,  <  .  1, 2.  ■  •  • ,  KTM),  {i?™,  i  - 1,2,  •  •  • ,  KTE) 

{/?™,  i  .  1,2,  -,KTM),{7?“,  i  -  1,2,  -,KTE).  (2.61) 

The  dements  cS  V/y/rj  where  V  is  the  column  vector  that  satisfies  the  matrix 
equation  (7.11|  are  printed  under  the  heading  V.  The  dements  of  are 
ordered  as  the  elements  of  the  excHatum  vector  are  ordered: 

{ i  *  1, 2, .  •  • ,  KTM),  i  *  1, 2, ... , KTE) 

t  =  1,2, . . .  ,KTM},  i »  1,2, . . .  ,KTE}.  (2.62) 

The  written  values  of  CIOUT,  CllN,  C20UT,  C2IN,  ClOUTS,  ClINS, 
C20UTS,  C2INS,  and  PT  are  described  in  Section  2.2.1. 

The  16  lines  following  the  line  where  PT  is  written  contain  the  normalized 
coefficients  of  the  TM  propagating  inodes  in  the  drcular  waveguide.  As  indi* 
cated  in  Section  2.2.1,  CTME  and  CTMO  are,  respectively,  the  normalized 
coefficients  of  the  reflected  waves  of  the  even  and  odd  TMr.i^  pn^>agating 
inodes  that  exist  in  the  circular  waveguide  of  Fig.  2.  Also,  CTMES  and 
CTMOS  are  the  squares  of  the  magnitudes  of  CTME  and  CTMO,  respec¬ 
tively.  Other  quantities  written  in  the  16  lines  cited  in  the  first  sentence  of 
this  paragraph  are  XTME  of  (7.146),  XTMO  of  (7.152),  XTMES  of  (7.150), 
and  XTMOS  of  (7.156). 

Now,  XTME  and  XTMO  are,  respectivdy,  the  normalized  coefficients 
<d  the  waves  of  the  even  and  odd  TMr-i^  propagating  modes  radiated  by 
the  combination  of  the  magnetic  currents  — and  in  the  circular 

waveguide  of  Fig.  2.  In  other  words,  XTME  is  the  contribution  to  CTME 
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due  to  the  magnetic  cunents  in  the  circular  waveguide,  and  XTMO  is  what 
the  magnetic  currents  contribute  to  CTMO.  Also,  XTMES  and  XTMOS  are 
the  squares  of  the  magnitudes  of  XTM£  and  XTMO,  respectivdy. 

The  next  44  lines  contain  the  normalised  coefficients  ci  the^propagating 
modes  in  the  circular  waveguide.  As  indicated  in  Section  2.2.1,  CTEE  and 
CTEO  are,  respectively,  the  normalised  coefficients  of  the  reflected  waves  of 
the  even  and  odd  TEr.i  j  propagating  modes  that  exist  in  the  circular  wave* 
guide  of  Fig.  2.  Also,  CTEES  and  CTEOS  are  the  squares  of  the  magnitudes 
of  CTEE  and  CTEO,  respectively.  Other  quantities  written  in  the  44  lines 
dted  in  the  first  sentence  of  this  paragraph  are  XTEE  oi  (7.186),  XTEO  ci 
(7.122),  XTEES  of  (7.190)  and  XTEOS  of  (7.196). 

Now  XTEE  and  XTEO  are,  respectively,  the  normalised  coefficients  ci 
the  waves  oi  the  even  and  odd  TEr.i4  propagating  modes  radiated  by  the 
combination  of  the  magnetic  currents  — and  — in  the  circular  wave* 
guide  of  Fig.  2.  In  other  words,  XTEE  is  the  contribution  to  CTEE  due  to 
the  magnetic  currents  in  the  circular  waveguide,  and  XTEO  is  what  the 
magnetic  currents  contribute  to  CTEO.  Also,  XTEES  and  XTEOS  are  the 
squares  of  the  magnitudes  of  XTEE  and  XTEO,  respectively. 

The  next  three  lines  show  the  values  (d  PIN,  PRM,  PR,  PTOTAL,  PTA, 
and  PRMA.  Here,  PIN  is  the  time*average  power  Piae  of  (5.44)  of  the  s> 
traveling  (incident)  waves  in  the  circular  waveguide,  PRM  is  the  time-average 
power  Prm  of  (5.34)  radiated  by  the  magnetic  currents  in  the  circular  wave¬ 
guide,  and  PR  is  the  time-average  power  Pr  of  (5.45)  reflected  in  the  circular 
waveguide.  Also,  PTOTAL  is  the  sum  of  Pr  and  the  time-average  power  Pt 
transmitted  into  the  rectangular  waveguides.  Here,  Pt  is  ^ven  by  (2.52)  in 
which  Pf^^  and  are  given  by  (4.7)  and  (4.8).  The  variable  PTA  is  Pu 
where  Pu  is  an  alternative  expression  for  the  time-average  power  transmitted 
into  the  rectangular  waveguides: 

+  (2.63) 

where  P^^  and  are  given  by  (4.12)  and  (4.17),  respectively.  The  variable 
PRMA  is  Prm«  where  Prm»m  is  the  alternative  expression  (5.43)  for  the  time- 
average  power  radiated  by  the  magnetic  currents  in  the  circular  waveguide. 

The  computer  program  variable  PT  described  in  Section  2.2.1  and  PIN, 
PRM,  PTOTAL,  PTA.  and  PR.MA  cited  in  the  first  sentence  of  the  previous 
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paragraph  should  satisfy 


PT*PTA 

(164) 

PIN  *  PTOTAL 

(2.65) 

PRM  »  PRMA. 

(2.66) 

The  numerical  values  of  these  o(»nputer  program  variables  in  the  eacemplary 
output  data  satisfy  the  above  equations  within  roundoff  error. 

The  J*^  of  the  NPHI  entries  written  under  the  PHI2  is  given 

by 

“NPHi-r 


The  of  the  NZ*^  entries  written  under  the  Z  is  given  by 


,(i)  _  (J  - 

^  ”  NZ-1  * 


(2.68) 


The  quantities  J  —  1, 2,  •  •  • ,  NPHI}  are  useful  abscissas  for  plotting  the 
ordinates  written  under  the  headings  E3A1PS  and  E3A2PS.  The  quantities 
J  ss  1, 2,  •  •  • ,  NZ)  are  useful  abscissas  for  plotting  the  cwdinates  written 
under  the  headings  E3A1ZS  and  E3A2ZS. 

The  of  the  NPHI  entries  written  under  the  heading  E3A1PS  is 
£3A1PS(J)  of  (2.54),  and  the  of  the  NZ  entries  written  under  the  heading 
E3A1ZS  is  E3A1ZS(J)  of  (2.55).  The  of  the  NPHI  entries  written  under 
the  heading  E3A2PS  is  E3A2PS(J)  of  (2.56),  and  the  of  the  NZ  entries 
written  under  the  heading  E3A2ZS  is  ^A2ZS(J)  of  (2.57). 

The  I^  of  the  KAM  entries  written  under  the  heading  BKAPLT  is  ka 
given  by  the  right-band  side  of  (2.15)  with  KA  replaced  by  I: 


ibas:BKAO-f  (I-1)*DBKA  (2.69) 

The  P**  of  the  KAM  entries  written  under  the  heading  PINC  is  the  value  of 
PIN  that  was  calculated  with  ka  given  by  (2.69).  The  I***  of  the  KAM  entries 
written  under  the  heading  PTRAN  is  the  value  of  PT  that  was  calculated 
with  ka  given  by  (2.69).  The  P**  of  the  KAM  entries  written  under  the 
heading  PREFL  is  the  value  of  PR  that  was  calculated  with  ka  given  by 
(2.69). 
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Output  data  written  in  the  file  BESOUT.DAT  by  itatementa  in 
the  mibroutine  BESIN 

These  output  data  are  the  same  as  those  described  in  Section  2.2.2  of  [3]. 

2.3  Minimum  Allocations 


The  fninimum  storage  space  that  must  be  allocated  to  some  arrays  in  the 
cfxnputer  program  depends  cm  the  values  of  the  input  variables  B,  C,  BKM, 
KAM,  KE3M,  NPHI,  NZ,  RTMR,  RTER,  (STMR(R),R»1,RTMR),  and 
(STER(R),  R  s  1,  RTER).  There  is  also  dependence  mi  the  Tnaximnm  of  the 
values  of  ka  in  (2.15).  Fiurthennore,  dependence  on  XM  arises  when  XM  is 
extraordinarily  large  (see  the  last  two  paragr^ihs  of  this  section). 

Minimum  allocations  in  the  main  program  are  given  by 

HM(liMAX)  «BMV(KTE)  .BN1I2(KTE)  .PRi(3*PNAX)  ,PH2(3*PHAX) . 

PR3(3«PMAX)  .PH4(3«PMAX)  .03(I1MAX)  .G4(1IMAX}  .OTM(IIMAX)  .DTE(VKAX) » 
miPCNPHI)  ,E3AlZ(irZ)  .E3A2PCNPRI)  ,E3A2Z(1IZ)  .y(K2*K2)  ,TI(K2) » 
V(K2)  ,Y11EC(K2)  .GTHCMMAX)  .GTECHMAX)  ,TNP(1IMAX)  .HOKIIMAX)  , 
TEPCMMAX)  .TEMCNMAX)  .OQTMdlMAX)  ,0QTE(NMAX)  ,PHI1(XPHI) , 

PHI2(NPHI)  .Z(irZ)  .PTRAN(KAM)  .PREFL(KAM)  ,BKAPLT(KAM)  .SINPCPMAX)  . 
SINQCPHAX)  ,E3A1PS(IIPHI)  .E3A2PS(1IPHI)  .E3AlZS(IfZ)  ,E3A2ZS(1IZ) , 
IPS(K2)  ,KE3(KE3M)  ,STM(RTM)  .STECKTE)  ,STMR(RTMR)  ,STER(RTER) , 
ITME(KTE*irrME)  ,ITMO(KrE*irmO)  .ITEEP(ICrE*IITEE)  ,nEEZ(iaE*irrEE) , 
ITEOP(KTE*KTEO)  ,ITE0Z(KTE*HTE0)  ,PINC(KAM)  ,BTM(S™m,RTM)  , 
BT£(STEH,BTE)  .HTNE(NAX(STMN.STMR)  ,MAX(RT1f,RTMR)), 

BTHO (MAX (STMM , STHR)  .MAX  (RIR .RTMR)  ) , 

BTEE (MAX (STEM, STER)  .MAX(RTE,RTER))  . 

BTEO (MAX (STEM. STER)  . MAX (RIE. RTER)) 


where  “MAX”  denotes  the  maximum  value  of  the  quantities  that  follow  in 
parentheses.  In  the  above  list  of  minimum  allocations, 


NMAX 


=  i+{ 


the  maximum  value  of  n  such  that 
(nx)  (B/C)  <  BKM 


) 


and 


KTE  =  -1  +  < 


the  number  of  combinations  of  nonnegative 
int^ers  m  and  n  such  that 

^(mx)’  +  ((nx)B/C^  <  BKM 


(2.70) 


(2.71) 
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FVirtlieniiofe, 


A  V  1  .  f  majritnum  value  <d  p  such  Uiat  1  .. 

PMAX  =  1  +  |  ,,<BKM  / 


K2  =  2*(KTM  +  KTE) 


where  KTE  is  given  by  (2.71)  and 


KTM  = 


the  number  of  combinations  of  positive 
integers  m  and  n  such  that 

Jimry  +  ((n»)(B/C))’  <  BKM 


Moreover, 


MTAiTF  —  /  number  of  positive  roots 

“  \  {Jry  r  ss  0, 1,2,  •  •  •}  such  that  x„ 


NTMO 


,of  \ 

<ka  f 


{the  number  of  positive  roots  Xr»  of  1 
{Jr,  r  s  1,2,*«>}  such  that  Zr«  <  ka  j 


MTPn  —  /  number  of  positive  roots  i',  Pf  1 
”  1  ^  *  1,2,***}  such  that  z^,  <  ho  J 

where  ka  is  the  maximum  value  of  the  right-hand  side  of  (2.15): 


.  /  BKAO, 

\  BKA0+(KA.M-1)* 


DBKA 
DBKA,  DBKA 


<0  I 

>0.  / 


(2.73) 


(2.74) 


(2.75) 


(2.76) 


v-TPir  /  number  of  positive  roots  z^,  of  \  . 

NTEE  *  I  ( j;,  ^  s  0, 1, 2, . . .}  such  that  z;.  <  ho  / 


(2.78) 


(2.79) 


In  (2.75)  and  (2.76),  Jr  is  the  Bessel  function  of  the  first  kind  of  order  r,  and 
its  positive  roots  {zr«,  s  as  1,2,  *  *  ■}  ^  ordered  such  that 


Zr]  <  ZrJ  <  Zr3  ' 


(2.80) 
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In  (2.77)  and  (2.78),  ia  the  derivative  of  with  respect  to  its  argument 


and  the  positive  roots  s  »  1, 2,  •  •  •}  of  are  ordered  such  that 

*rl  <*U*”  •  (2.81) 

Other  variables  appearing  in  the  mminwim  allocations  presented  in  the 
previous  paragraph  are 

RTM  3s  1  +  (the  fnaximum  value  of  r  in  (2.75))  (2.82) 

STMM  s  the  maximum  value  of  s  in  (2.75)  (2.83) 

BTE  s  1  +  (the  maximum  value  of  r  in  (2.77))  (2.84) 

STEM  s  the  maximum  value  of  s  in  (2.77) .  (2.85) 


If  there  are  no  roots  in  (2.75),  then  RIM  s  1  and  STMM  s  0.  There  is 
at  least  one  root  in  (2.77);  the  root  Xjj  has  to  be  allowed.  If  it  were  not 
allowed,  then  there  would  be  no  propagating  mode  with  which  to  excite  the 
circular  waveguide. 

Now,  NMAX  of  (2.70),  KTE  of  (2.71),  PMAX  of  (2.72),  K2  of  (2.73), 
KTM  of  (2.74),  RTM  of  (2.82),  and  ^E  of  (2.84)  are  variables  in  the  main 
program,  but  NTME  of  (2.75),  NTMO  of  (2.76),  NTEE  of  (2.77),  NTEO 
of  (2.78),  STMM  of  (2.83),  and  STEM  of  (2.85)  are  not.  All  13  variables 
mentioned  in  the  previous  sentence  are  meant  to  be  int^er  variables.  As 
defined  by  (2.71),  KTE  is  the  number  of  TE  modes  retained  in  the  expansion 
of  the  field  in  each  rectangular  waveguide.  Furthermore,  KTM  of  (2.74)  is 
the  number  of  TM  modes  retained  in  the  expansion  of  the  field  in  each 
rectangular  waveguide.  The  quantities  NTME  of  (2.75),^NTMO  of  (2.76), 
NTEE  of  (2.77),  and  NTEO  of  (2.78)  are,  respectively,  the  numbers  of  even 
TM,  odd  TM,  even  TE,  and  odd  TE  modes  that  propagate  in  the  circular 
waveguide. 

One  who  is  not  inclined  to  calculate  the  precise  values  of  NMAX,  KTE, 
PMAX,  K2,  KTM,  NTME,  NTMO,  NTEE,  NTEO,  RTM,  STMM,  RTE, 
and  STEM  from  their  definitions  is  advised  to  use  educated  guesses  for  these 
variables  in  order  to  specify  allocations.  If  one  can  get  the  program  to  run 
to  completion  without  any  error  messages,  then  he  can  use  the  output  data 
to  easily  obtain  the  precise  values  of  the  variables  mentioned  in  the  previous 
sentence.  The  actual  values  of  NMAX,  KTE,  KTM,  RTM,  and  RTE  appear 
in  these  output  data.  The  first  integer  number  written  under  the  heading  MM 
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is  PMAX.  The  value  of  K2  is  twice  that  of  the  printed  value  of  Kl.  M<mover, 
NTME  is  the  sum  of  the  integer  numbers  written  under  the  heading  STM, 
NTMO  is  this  sum  minus  the  first  integer  number  written  under  the  heading 
STM,  NTEE  is  the  sum  of  the  integer  numbers  written  under  the  heading 
STE,  and  NTEO  is  this  sum  minus  the  first  integer  number  written  under 
the  heading  STE.  Moreover,  STMM  is  the  largest  of  the  integer  numbers 
written  under  the  heading  STM,  and  STEM  is  the  largest  of  the  integer 
numbers  written  under  the  heading  STE. 

Minimum  allocations  in  the  subroutines  are  given  by 

MODES, 

PHI 


DON, 

DECOMP,  and 
SOLVE. 

A  blank  or  labeled  common  block  that  is  used  in  two  or  more  program  seg¬ 
ments  must  be  defined  exactly  the  same  in  each  of  these  program  segments.^ 
Therefore,  any  dimensioned  variable  in  a  common  block  that  is  used  in  two 
or  more  program  segments  must  have  the  same  allocation  in  each  of  these 
program  segments.  That  is  why  each  of  the  arrays  PHI,  PH2,  PH3,  and  PH4 
in  the  subroutine  PHI  was  allocated  3*PMAX  entries  instead  of  2*PMAX 
entries. 

If  XM  is  so  large  that 


/  MM(NMAX),BMN(KTE),  \ 

1  BMN2(KTE)  /  “ 

/  PH1(3*PMAX),PH2(3*PMAX),  1  . 

\  PH3(3  ♦  PMAX),PH4(3  •  PMAX)  |  ““^routme 

D3(NMAX),G4(NMAX),D(NMAX),  ) 

CP(NMAX),CM(NMAX),  \  in  the  subroutine 
DQ(NMAX),G(NMAX)  J 

{  UL(K2  *  K2),SCL(K2),IPS(K2)  }  in  the  subroutine 

{  UL(K2  •  K2),B(K2),X(K2),IPS(K2)  }  in  the  subroutine 


<  XM  (2.86) 

where  is  the  200*^  root  of  then  the  “200”  in  the  fifth  statemmt  in 
DO  loop  19  of  the  main  program  must  be  replaced  by  I  where  I  is  an  integer 

^Hck,  a  program  aegment  is  cither  the  main  program  of  one  of  the  subprogram. 
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at  least  so  large  that^ 


x;j>XM.  (2^) 

AlaOy  the  maxifniini  value  of  the  index  S  of  DO  loop  25  in  the  eubroutine 
BESIN  nnut  be  increased  from  200  to  I.  Fbrthermoee,  the  maxinnim  value 
d  the  indices  S  d  DO  loops  14  and  15  in  the  subroutine  BES  must  be 
increased  from  200  to  I.  Accompanying  minimam  allocatioiis  axe  pven  by 


{XJ(I)XJP(I)} 

{A(I)^P(I)} 

A(I)^P(I),  1 
XJ(I)^P(I)  / 

{A(I)^P(I)} 

{X(I)} 


in  the  m^  program, 
in  tlM  subroutine  BESIN, 

in  the  subroutine  BES, 

in  the  subroutine  INTERPOL,  and 
in  the  subroutine  DON. 


If  XM  is  so  large  that 


*isaa  ^  XM 


(2.88) 


where  z^o,!  is  the  first  root  of  Jig,,  then  the  upper  limit  on  the  index  of  DO 
loop  19  in  the  main  program  must  be  increased  from  500  to  J  where  J  is  an 
integer  at  least  so  large  that 


zjj  >  XM. 


(2.89) 


^See  the  last  paragraph  d  Secuoe  34  of  {31  Note  that  the  first  three  I’s  in  the  cnrly 
brackets  in  the  line  alter  eg  (3  73)  of  JS]  are  rriooioos.  Any  integer  variable  other  than  I 
is  meant  here. 
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Chapter  3 

The  Excitation  Vector 


The  waveguide  mode  converter  ia  excited  by  a  mixture  of  the  z*traveling 
waves  of  the  propagating  modes  in  the  circular  waveguide.  The  elements  of 
the  excitation  vector  are  obtained  by  substituting  the  magnetic  field  of  these 
waves  and  their  reflections  about  the  sh<»t  at  z  s  £3  into  eq.  (2.24)^  of  [1]. 
This  magnetic  field  is  that  whidi  would  exist  in  the  circular  waveguide  if 
both  apertures  were  shorted. 


3.1  The  Electromagnetic  Field  of  the 
Mixture  of  2-TVaveling  Waves  in  the 
Circular  Waveguide 

The  electric  and  magnetic  fields  of  the  mixture  of  z-traveling  waves  in  the 
circular  waveguide  are  and  given  by 


gmj^ 


^An  equation  that  appeal*  in  a  lefercace  i*  ated  by  placing  "eq.”  before  the  equation 
number.  An  equation  that  appear*  in  the  prcwnt  report  is  cit^  by  writing  only  its 
number  in  parentheses  Bowevcr.  an  equMton  number  at  the  beginning  of  a  sentence  is 
always  preceded  by  *EquatMo*. 
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(3.1) 


*S£(^^)'^' 


(3.2) 


where  the  B's  *r»  joastaats,  i»  givea  by  eq.  (B.l)  of  [1], 

is  givea  by  eq.  (B.26)  of  [1],  is  given  by  eq. 

(B.35)  of  [11,  and  {£S^,lC^)  is  given  by  eq.  (B.55)  of  (1).  The  hrst 
sununation  with  respect  to  the  indices  r  and  s  on  the  ri^t-hand  side  of  each 
of  (3.1)  and  (3.2)  is  taken  over  all  the  propagating  TM  modes  of  the  drcular 
waveguide.  The  second  sununation  with  req>ect  to  the  indices  r  and  s  on  the 
right-hand  side  of  each  of  (3.1)  and  (3.2)  is  taken  over  all  the  propagating 
TE  inodes  of  the  circular  waveguide. 

The  quantity  that  multiply  on  the  right-hand  side  of  (3.2)  re¬ 
duces  to  the  magnetic  field  when  z  ^  L^,  This  magnetic 

field  is  that  of  a  TM  wave  whose  time-average  power  is  imity.  The  quantity 
that  BJ^^  multiplies  on  the  right-hand  side  of  (3.1)  reduces  to  the  electric 
field  I\fy^^  when  z  =  £3.  This  electric  field  is  that  of  a  TE  wave 
whose  time-avera^  power  is  unity.  Here, 

2i^,  and  arc  given,  respectively,  by  eqs.  (B.22),  (B.23),  (B.33), 
(B.34),  (B.51),  (B.52),  (B.62),  and  (B.63)  of  [1).  Moreover,  and 
are  given,  respectively,  by  eqs.  (B.25)  and  (B.54)  of  [1]  where 


(3.3) 

(3.4) 

In  (3.3)  and  (3.4),  k  is  the  wavenumber  and 

(3.5) 

tTE  _  *r* 

(3.6) 
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where  x„  is  the  root  of  the  Bessel  fnnctkm  Jr  and  is  the  s*^  root  aS 
J^,.  Here,  is  the  derivative  of  J,  with  respect  to  its  argument. 


3.2  The  Electromagnetic  Field  That  Would 
Exist  in  the  Circular  Waveguide  if 
Both  Apertures  Were  Shorted 

If  the  ^>ertures  Ax  and  in  Fig.  1  were  closed  with  perfect  conductMS,  then 
the  electr(»nagnetic  field  circular  wave¬ 

guide  would  be  the  r-traveling  waves  cm  the  ri^t-hand  sides  of  (3.1)  and 
(3.2)  plus  their  reflections  about  z  »  X3.  Adding  these  reflections  to  the 
right-hand  sides  of  (3.1)  and  (3.2),  we  obtain 


+  E  L  (3.T) 

+2: 2: 

ra  fm,^  \  yJY^^  /  '  ' 


where  (£?^,i£^-).  «id  (g^, 

HJ^)  are  given,  rcspwlively,  by  eqs.  (B.2),  (B.27),  (B.36),  and  (B.56) 
of  [1].  Obviously,  and  £‘®>(.Z^,fl)  of  (3.7)  and  (3.8)  are  propa- 

gating  mode  contributions  to  the  total  electromagnetic  field  in  the  cylindrical 
waveguide.  In  Chapter  4,  (3.7)  will  be  used  in  obtaining  expressions  for  the 
coefficients  of  the  propagating  modes  of  the  total  electromagnetic  field  in  the 
cylindrical  waveguide. 
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3.3  The  Excitation  Vector  in  Terms  of  the 
B’s  in  (3.1) 

Tlie  excitation  vector  is  the  column  vector  on  the  right>haad  side  of  eq.  (2.22) 
of  [1]  where  the  dement  of  is  pven  by  eq.  (2.24)  of  [1]: 

JT*  -  -  /jf^  jd  •  u  (3.9) 

where  t  is  related  to  m  and  n.  The  relatkmship  between  t  and  m  and  n  is  the 
same  as  in  eq.  (2.24)  of  [1].  Substituting  eqs.  (B.l),  (B.2),  (B.26),  (B.27), 
(B.35),  (B.36),  (B^),  ud  (B.56)  of  (1)  lot 

a?**.  Md  ajf^  in  (3.8), « obuin 

-EE  +  E  E  (3.10) 

ft*  SmtA  ft* 


where 


7^ 


A?''eo.(/J™(I,-t)) 


(3.11) 


=  -s^{y',?“ii7'co.(d”(L,-t)) 


-  *))}• 

Substitution  of  (3.10)  into  (3.9)  gives,  upon  multiplication  by  —  j\/q, 


(3.12) 


(3.13) 


where 


a=  1,2 

^  tM^TE 

,  /  =  «,o. 
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Here,  if  »  ^/i/e  is  the  intzinsic  impedance.  In  (3.14),  >>  lC{ 

anted  at  p  as  a.  Suitable  expressions  for  ACw  must  be  found 

before  the  integral  on  the  right-hand  side  (rf  (3.14)  can  be  evaluated. 

3.4  Evaluation  of  Expression  (3.14) 

3.4.1  The  M's  in  (3.14) 

FVom  eq.  (2.13)  of  (1], 


(3.15) 

(3.16) 


where  e^„  and  ef^„  are,  respectively,  the  p-  and  s-cmnpcments  of  Here, 
0  is  either  TM  or  TE,  and  fi™  *und  are  given  by  eqs.  (A.IO)  and  (A.23) 
of  [1].  In  (3.15)  and  (3.16), 


where 


y*+  *  (X  -  ^)*.  +  i 

(3.1T) 

-1- 1 

(3.18) 

r'*’  =  z  -t- 1 

(3.19) 

asin^o 
*•-  ♦.  • 

(3.20) 

With  and  given  by  eqs.  (A.IO)  and  (A.23)  of  |1],  (3.15)  and  (3.16) 
expand  to 


41 


MS . 

-('7r)(T)“(=^)*{=?))  i“) 

-(Tr)©”(=?^)-(=?)>- 


3.4.2  The  KH  in  (3.14) 

Substitution  of  eqs.  (B.23)  and  (B.34)  of  [1]  into  (3.11)  gives 

rj,(e*.)m(r») 


^  ~  kVf 

+l4'^(*?‘/»)cos(r4)}  cos(/9j‘(l3  -  *)) 


(3.25) 


"~2r“  (  2  \  r  rJ.(fcyp)cos(r4) 

t„/9™  [aJr^tiXr.))  kj^p 


~1U*^(*™^)»«(»’4)}  cos(;9™(l3  -  *)) . 
In  obtaining  (3.25)  and  (3.26),  we  used 


(3.26) 


ZTM«o  _ 
r$  — 


(3.27) 


where  k  s  ttiy/pi  is  the  wavenumber  and  17  s  yjp/t  is  the  intrinsic  impedance. 
Equation  (3.27)  results  when  7^  of  (3.3)  is  substituted  into  eq.  (B.25)  of 


« 


42 


[ll-  Sabttitatioo  o£  eqs.  (B.41),  (B.52),  (B.59),  ud  (B.63)  of  [1]  into  (3.12) 

glVM 


^  (ife) 

tin(rd)]  cotifi^iLi  - 1}) 

(^)  «»(»’^)  *in(^“(Is  -  *))}  (3.28) 

^  (zfe) 


V# 

lTB' 


(^)  «»03?(^«  “*))}•  (3.29) 

In  obtaining  (3.28)  and  (3.29),  we  used 

rE 

(3.30) 


yTE«e  _  ^ 

'•  ”  ib7  • 


Equation  (3.30)  results  when  7^  of  (3.4)  is  substituted  into  eq.  (B.54)  of  (1). 
Because 


M^r,)  =  0,  (3.31) 

substitution  of  x  =  x„  into  the  identity  (see  formula  9.1.27  of  (4)) 

Jr^lix)  =  -j;(x)  +  ^Jr{x)  (3.32) 

gives 


^r+l(x„)*-j:(x„). 


(3.33) 
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When  ^  »  a,  it  follows  from  (3.5) 

(3.34) 

and  consequently, 

^(Crt  -  0 

(3.35) 

j;(*?‘<»)  -  j;(»-)- 

(3.36) 

Setting p^am  (3.25)  «nd  (3.26) and naiiig (3.33),  (3.35),  and  (3.36),  we 
obtain 


\^]^- .))  (3.37) 

=  (s-») 

Recall  that 


(3.39) 

When  p  s  a*  it  follows  from  (3.6)  that 

kJ^P  *  (3.40) 

and  consequently, 

Mf^p)  =  J.(*;.)  (3.41) 

r,{kl^p)  -  0.  (3.42) 

Setting  p  =:  a  in  (3.28)  and  (3.29)  and  using  (3.6)  and  (3.40)-‘(3.42),  we 
obtain 


(^a)  ”  *))} 


sin(r^)cos(^“(l3-*)) 


(3.43) 


44 


[^u= fe) 

+“•  (^^)  8m(^™(I.s  -*))}.  (3.44) 


Wanting  to  express  (3.37),  (3.38),  (3.43),  and  (3.44)  in  terms  of  and 
X*  rather  than  4  fint  sdve  (3.17)  for  ^  in  terms  of  y^*: 


Defining 


(3.45) 


(3.46) 


and  recalling  (3.19),  we  have 

L,-z~14-z*.  (3.47) 

Substitution  of  (3.43)  tnd  (3.47)  into  (3.37),  (3.38),  (3.43),  and  (3.44)  giues 

(!)  -  »’*)) 

(3.48) 

(!)  «i>  (^(f  -  »’*))  eo.(^?*(iJ  -  '*)) 

(3.49) 

(t*)  (i)  ™  (i't  - »-)) 

•coeifiJ^iLt  -  2*))  +  M,  «*•  (~(|  “  y*"^))  «“(^J^(^3  -  ^■^))} 

(3.50) 

{^)  “*  te)  “-(ii^’  - 
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4 


•eo.W™(£J  -  *+))  +  «.  sin  (i(i  -  ,■+))  siii(^™(ij  -  s+))}. 

(3.51) 

ExpresrioDB  (3.48)-(3.51)  are  valid  only  in  A\.  To  obtain  exprankns  valid 
in  Alt  we  first  solve  (3.18)  for  ^  in  terms  oi 

*  =  (3^) 

Sttbstituticm  of  (3.52)  and  (3.47)  into  (3.37),  (3.38),  (3.43),  and  (3.44)  gives 

(!)  “•  (!:<»”  -  ‘0 

(»•«) 

(¥)  {“♦  fe)  iry  - 1)) 

.co.(/9«(Lj  -  »♦))  +  u.  (^)  CO.  -  {))  «a(^"(ij  -  .+))} 

(3.55) 


[a-*] 


3gg 


co.(/J:!«(4*  -  .-*))  +  «.  «n  (^(»>*  - 1))  .in(^?(^  -  .+))}. 

(3.56) 


Expressions  (3.53)-(3.56)  are  valid  in  Ai. 


3.4.3  Substitution  of  the  M’s  and  £’s  of  Sections 
3.4.1  and  3.4.2  into  (3.14) 

In  (3.14),  ds  s  athdz.  Substituting  (3.21),  (3.22),  and  (3.48)  into  (3.14) 
and  changing  the  variables  of  integration  from  4  end  z  to  and  z*,  we 
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obtain 


where 


in  whidi 


z/™-™*  -  [-lYniSf 

(3.57) 

(3.58) 

iHHr = I®**-*?* 

(3.59) 

yTM 

(3.60) 

y«e  *= 

(3.61) 

,TM  _ 
*as  *" 

!  jT  “•  (^)  “• 

(3.62) 

In  obtaining  (3.57),  we  replaced  2  by  y/tm^  in  (3.21).  Such  replacement  is 
allowed  because  the  right-hand  side  of  (3.21)  vanishes  when  dther  n  or  m 
is  zero.  The  coefficient  is  so  named  because  it  is  prt^rtional  to  the 
quantity  that  the  fiven  TM^  modal  part  of  the  incident  magnetic  field  in 
the  circular  waveguide  contributes  to  — and  of  (3.13).^ 

As  defined  by  (3.61),  y^e  is  an  integral  with  respect  to  of  the  product  of 
a  aine  and  a  cosine.  As  defined  by  (3.62),  is  an  integral  with  respect 
to  c***  of  a  cosine  and  a  cosine  whose  argument  is  proportional  to  the  TM 
propagation  constant  0^.  Substitution  of  (3.23),  (3.24),  and  (3.53)  into 
(3.14)  gives 

7.*™  tm.  ^  ^;TM.  (3  g3) 

j  ,3 

^Rather  than  indicating  vectocs,  the  underlines  in  this  sentence  bi^ili^t  the  super¬ 
scripts  and  subscripu  of 
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where  1^**  is  given  by  (3.59)  in  which  y^c  i>  given  by  the  right-hsnd  side  of 
(3.61)  with  replaced  by  y^*^.  Since  y**^  and  yl**  are  dnnm^r  variables  oi 
integration,  it  follows  that  the  y^  implicit  in  (3.63)  and  (3.64)  is  given  by 
the  unaltered  right-hand  side  of  (3.61).  Therrfore,  the  1^^  in  (3.63)  and 
(3.64)  is  predsriy  that  given  by  (3.59)  in  whidi  7^,  y^,  and  are  pven 
by  (3.60),  (3.61),  and  (3.62),  respectivdy.  In  obtaining  (3.63),  we  replaced  2 
by  y/mm  in  (3.23).  Sudi  replacement  is  allowed  because  the  right-hand  side 
of  (3.23)  vanishes  when  either  n  or  m  is  zero. 

Substituting  (3.21),  (3.22),  and  (3.49)  into  (3.14),  we  obtain 

^  (3  65) 


where 


rTMo  _  r™.. 

2^.  /*  .  ft 


(3.65) 

(3.66) 

(3.67) 


In  (3.67),  and  zj^  are  given  by  (3.60)  and  (3.62).  In  obtaining  (3.65) 
and  (3.^),  we  replaced  2  by  Cr  in  (3.49).  Such  replacement  is  allowed  because 
the  right-hand  side  of  (3.49)  vanishes  when  r  s  0.  The  coefficient  is  so 
named  because  it  is  proportional  to  the  quantity  that  the  odd  modal 
part  of  the  incident  magnetic  field  in  the  drcular  waveguide  contributes  to 
and  (3.13).  As  defined  by  (3.68),  y^  is  an  integral 

with  respect  to  c(  the  product  of  a  lio^  vid  a  sine.  Substitution  of  (3.23), 
(3.24),  and  (3.54)  into  (3.14)  gives 

yamTM.  ^  (3  69) 


rZTC.TM* 


(2f  ^  7™* 
5  /  • 


(3.70) 


In  obtaining  (3.69)  and  (3.70),  we  replaced  2  by  Cr  in  (3.54).  Such  replace¬ 
ment  is  allowed  because  the  right-hand  side  of  (3.54)  vanishes  when  r  =  0. 
Substituting  (3.21),  (3.22),  and  (3.50)  into  (3.14),  we  obtain 


riTM.TC« 


T\n.n* 


(3.71) 

(3.72) 
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where 

jTM  ^  r”ry.,»?  (3.73) 

“  X  jT  “•  (^^)  “•  (ir'*'"  ■  *0  <*•”> 

z“  =  ?j(‘co.(=^jco.(^”(«-j+))dz+  (3.77) 

z™  =  ?  nn  (2^  j  no  (^™(IJ  -  z+))  <U* .  (3.78) 


In  (3.73),  Vm  is  given  by  (3.68).  The  quantity  is  proportional  to  the 
quantity  that  the  ^component  of  the  £ven  TE^  modal  part  of  the  incident 
magnetic  field  contributes  to  — of  (3.13).  The  quantity  is 
proportional  to  the  quantity  that  the  ^‘component  of  the  even  TE^,.  modal 
part  of  the  incident  magnetic  field  amtributes  to  of  (3.13).  Sub* 

(3.79) 

(3.80) 

(3.81) 

(3.82) 

(3.83) 

(3.84) 


stituting  (3.23),  (3.24),  and  (3.55)  into  (3.14),  we  obtain 

jTTM.TE*  ^  j jr^l^lTM.TEe 

^  ^_jp|^lTE.TEe 

Substitution  of  (3.21),  (3.22).  and  (3.51)  into  (3.14)  giv« 

I  jT^  +  „/TE~| 


where 


/.ijr  -  r™r,.z? 


/TEm 
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where 


In  (3.83)  and  (3.84),  T^,  y*,  xj®,  and  •*«  given  by  (3.75),  (3.61),  (3.77), 
and  (3.re),  respectively.  In  (Staining  (3.81)  and  (3.82),  we  rq>laoed  2  by  Cr 
in  (3.51).  Such  replacement  is  allowed  because  the  ti^t*hand  side  of  (3.51) 
vanishes  when  r  ae  0.  The  quantity  i*  proportMoal  to  the  quantity 
that  the  j^-component  of  the  fidd  I£^  modal  part  of  the  incident  magnetic 
field  contributes  to  — of  (3.13).  The  quantity  is  proportional 
to  the  quantity  that  the  ^-component  of  the  odd  XE^  modal  part  of  the 
incident  magnetic  field  contributes  to  of  (3.13).  Substituting 

(3.23),  (3.24),  and  (3.56)  into  (3.14),  we  obtain 


j3TM,TE«  _ 
^2TE.TEo  ^ 


(3.86) 

(3.87) 


In  obtaining  (3.86)  and  (3.87),  we  replaced  2  by  e^  in  (3.56).  Such  replace* 
ment  is  allowed  because  the  right-hand  side  of  (3.56)  vanishes  when  r  k  0. 


3.4.4  Evaluation  of  the  y’s  Defined  by  (3.61),  (3.68), 
(3.76),  and  (3.85)  of  Section  3.4.3 


Since 

cos  “  2))  =  co8(r0p)co8  +  sin(r^o)rin  (3.88) 

where,  in  agreement  with  eq.  (E.9)  of  [1], 


4o  = 


(3.89) 


(3.61)  and  (3.76)  become 

y«e  =  cos(r^«)  +  sin(r^p)  (3.90) 

Vec  =  «»(r^o)  +  ^2’  sin(r^p)  (3.91) 
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where  end  ere,  respectively,  and  of 

eqs.  (E.14)-(E.17)  of  [1]  with  p  replaced  by  m.  Since 


where  is  given  by  (3.89),  (3.68)  and  (3.85)  become 

y„  *  coe(r^.)  -  sin(r^,)  (3.93) 

Vc  =  coe(r^.)  -  sin(r^,) .  (3.94) 

Equations  (3.90),  (3.91),  (3.93),  and  (3.94)  can  be  rewritten  as 

»«  =  #2  (3.95) 

»«  »  *1  (3.96) 

»-  -  #2  (3.97) 

»«  =  *2  (3.98) 


where  ♦St  ^2t  “8  ♦“  are  given,  respectively,  by  eqs.  (6.28)-<6.31)  of 
[2]  with  p  replaced  by  m. 

3.4.5  Evaluation  of  the  r’s  Defined  by  (3.62),  (3.77), 
and  (3.78)  of  Section  3.4.3 

Equations  (3.62)  and  (3.77)  are  written  compactly  as 

4  =  I  jTcos  c09(St,{Lt  -  Z*))dz*,  6  =  TM,TE  (3.99) 

which  expands  to  (see  formula  401.06  of  (5]) 

i  jC  I  +  “•("‘■c*  +  )}  (3.100) 


where 


Pafuming  the  integratioa  in  (3.100),  we  obUin 

,  «ii(n‘*c  -  +  «■(«.«)  ,  dB(n‘-e  +  ^.U)  - 


“  n*+c  n^c 

(3.103) 

Concerned  ebout  roundoff  error  when  |n^|  is  small,  we  use  fcwmules  401.01 
and  401.02  of  [5]  and  the  identity 


cos X  s  1  —  2sin* 

to  recast  (3.103)  as 

,  sin(n**c)  cos(^,Xj‘)  -f  2  sin^  (»^)  sin(j9*  X?) 


(3.104) 


^  sin(n*-c)cos(^,£^)  -  2sin^  (aft)  sin(^,Xj') 

n^c  '  *  ' 

which  can  be  rewritten  as 

4  =  <51  (3.106) 

where  <51  is  the  tight>hand  side  of  eq.  (3.76)  of  [2]  with  q  replaced  by  n. 
Applying  formula  (401.07)  of  [5]  to  (3.78),  we  obtain 

e?  =  i  jC  {cos(n'^.-+  -  )  -  coe(o’^e+  +  iz* .  (3.107) 

The  evaluation  of  (3.107)  is  similar  to  that  of  (3.100):  . 

^  nn{n^c) )  +  2 sin*  (®^)  8in(^f^ ) 


s^n'^-c)  cosily !$•)  -  2 sin*  (°^)  sin{0^Lt) 


flTE-c 


(3.108) 


Equation  (3.108)  can  be  rewritten  as 

-JE  ^  _^^4)  (3  109) 

where  is  the  right-hand  side  of  eq.  (3.108)  of  [2]  with  q  replaced  by  n. 
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3.4.6  Utilization  of  the  Newly  Evaluated  y's  and  zH 

Substituting  (3.95)  nnd  (3.106)  into  (3.59),  we  obtain 

jTM.  ^  (3  110) 

where  7^  is  given  by  (3.60).  Substituticm  of  (3.97)  and  (3.106)  into  (3.67) 
gives 

JTM.  ^  ^3  Hj) 

Substitution  of  (3.97)  and  (3.106)  into  (3.73),  substitution  of  (3.96)  and 
(3.109)  into  (3.74),  substitution  of  (3.95)  and  (3.106)  into  (3.83),  and  substi¬ 
tution  of  (3.98)  and  (3.109)  into  (3.84)  yield 


(3.112) 

II 

a. 

(<•] 

ViSSFaJ 

1 

(3.113) 

(3.114) 

i7fr  =  (^) 

1 

(3.115) 

where  is  given  by  (3.75).  In  terms  of  the  I's  of  (3.110)-(3.115),  the  final 
expressions  for  of  (3.14)  are  (3.57),  (3.58),  (3.63)-(3.66),  (3.69)-(3.72), 
(3.79)-(3.82),  (3.86),  and  (3.87). 

3.5  Comparison  of  the  Elements  of  the 
Excitation  Vector  with  Those  of  [2] 

The  elements  rather  than  the  elements  of  the  excitation 

vector  itself  are  calculated  and  stored  in  the  array  TI  in  the  computer  pro¬ 
gram  of  the  present  report.  If  =  1  and  all  the  other  B's  in  (3.13)  are 
zero,  then  (3.13)  reduces  to 

= I  ^  =  tm,TE.  (3.116) 
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On  the  other  hand,  the  dements  {— and  {— 
given  by  eqs.  (4.8)  and  (4.9)  oi  [2]  are  calmlated  and  stoted  in  the  array  T1 
in  the  computer  program  d  [2].  HencefiMthjthe  ekmenta  mentioned  in  the 
previous  sentence  will  be  calle^— and  {— and 
written  condsely  as  {— where  the  that  ^ipears  in  a/9  is  either 
TM  or  T£.  In  tins  section,  we  relate  the  dements  {— of  (3.116)  to 
the  dements  {— 

When  *  1  and  all  the  other  B*»  are  sero,  the  exdtation  in  the 
present  report  is,  accmding  to  (3.2),  the  r-travding  wave  whose  magnetic 
fidd  is  JZ<*+>(^,11)  given  by 

(3.117) 

V^w** 

On  the  other  hand,  the  exdtation  in  [2]  is  the  ^-traveling  wave  whose  mag* 
netic  fidd  is  given  by  eq.  (5.2)  of  [1].^  Henc^  the  exdtatkm  cX.  the 

present  report  is  that  of  [2]  multiplied  by  •  Since  the  de¬ 

ments  of  the  exdtation  vector  are  proportional  to  the  exdtatimi,  it  fdlows 
that 


(3.118) 


Mtiltiplying  both  sides  of  (3.118)  by  ^jy/rj  and  using  (3.27),  we  obtain 


(3.119) 


In  Appendix  A,  we  verify  that  as  given  by  (3.116)  satisfies  (3.119) 

with  — given  by  the  right-hand  sides  of  eqs.  (4.8)  and  (4.9)  of  (2j. 


^The  excitation  in  [2]  is  the  same  as  that  in  [1]. 
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Chapter  4 

The  Propagating  Modes  in  the 
Rectangular  Waveguides 


4.1  The  Normalized  Coefficients  of  the 

Propagating  Modes  in  the  Rectangular 
Waveguides 


Dividing  both  sides  of  eqs.  (5.30)  and  (5.31)  of  [2]  by  we 

obtain  expressions  for  the  electric  fields  and  in  the  rectangular 
waveguides: 


evanescent  waves 


(4.1) 


evanescent  waves 


(4.2) 


whne  the  C’s  are  those  of  {2|  divided  by  The  evanescent 

waves  in  (4.1)  are  those  of  all  the  evanescent  noodes  in  the  left-hand  rectan¬ 
gular  waveguide  in  Fig.  2.  The  evanescent  waves  in  (4.2)  are  those  of  all  the 
evanescent  modes  in  the  right-hand  rectangular  waveguide  in  Fig.  2.  In  view 
of  (3.27)  and  eqs.  (5.32),  (5.33),  (5.36),  and  (5.37)  of  (2],t  the  (7’s  in  (4.1) 


^The  kft-hsnd  tide  of  eq.  (5.37)  of  [2]  should  be  C^~. 
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•ad  (4.2)  ue  given  by 


(.IE]  I  (Z.KP  +  1)*^ 

yi  k )  i2(^i??s:oS5n+7S(S5n) 
LIE]  ( _ 

(V  t  J  \2(Z,l^c<«(A.x.)  +  i«B(i9„.,)) 


(«) 

(4.4) 


LIE]  ( _ _ \  fisZ) 

(y  i  J  1 2(Z,ia*  co.(A*«,)  +  ;•  oniA^,))  f  \^  ) 

(  [^]  f  (Z,y^  - 1).-<*»"  )  (V!P\ 

\y  t  )  \2(Z,llPco<Ao».)  + jii»(A.*a))/ \,,fi  ) 


(«) 

(4.6) 


The  C*«  given  by  (4.3)-*(4.6)  we  the  nonnaliaed  coefficients  ti  the  propagat* 
ing  modes  in  the  rectsngulw  waveguides. 


4.2  The  Time- Average  Power  of  the 

Propagating  Modes  in  the  Rectangular 
Waveguides 

The  time-aversM  power  transmitted  into  the  left-hand  rectangulw  wav^uide 
<rf  Fig.  2  is  given  by 

ic;.^r-ic{j**r-  (4-T) 

Here,  the  subscript  stands  for  ’Hransmitted”.  The  time-avwage  power 
(4.7)  is  also  the  time-average  power  radiated  the  magnetic  current 
in  the  left-hand  rectangulw  waveguide  of  Fig.  2. 

The  time-average  power  tansmitted  into  the  right-hand  rectangulw  wave¬ 
guide  of  Fig.  2  is  given  by 

-  lcs?®-f .  (4.8) 

The  time  average  power  (4.8)  is  also  the  time-average  power  radiated  by  the 
magnetic  current  in  the  right-hand  rectangulw  waveguide  Fig.  2. 
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4.3  Alternative  Expressions  for  the 

Time- Average  Power  Radiated  by  the 
Magnetic  Currents  in  the  Rectangular 
Waveguides 


An  altenutive  expression  for  the  time-average  power  radiated  by  the  mM- 
netic  current  in  the  left-hand  rectangular  waveguide  of  Fig.  2  is 
given  by 

flL"  -  -Re  («.9) 


where  “Re”  denotes  the  real  part,  denotes  the  complex  conjugate,  and 
is  given  by  eq.  (2.11)  of  [1]  whidi  is  expressed  in  ^breviated  notation 


(4.10) 

In  (4.9),  the  subscript  “a”  stands  for  “alternative”.  Substitution  of  (4.10) 
into  (4.9)  gives 

fi‘'-Re{5:yriv^l’}  (<•») 

where  is,  with  reference  to  eq.  (2.25)  of  [1],  the  tt***  element  of  the  matrix 

’  y-l.lTM.lTM  y'l.lTM.lTE  T 
y-l.lTE.lTM  yl.lTE.lTE 

The  above  matrix  is  diagonal.  The  subscript  “rec”  in  (4.11)  stands  for 
“rectangular  waveguide".  Equatimi  (4.11)  is  recast  as 

{Kl  Y  2) 

s;(->«?nr)  ^  (4.12) 

where  “Im”  denotes  the  imaginary  part.  Numerically,  of  (4.12)  should 
be  equal  to  P/**  of  (4.7). 
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An  altemntive  expremon  ior  the  time>avence  power  radieted  by  the 
magnetic  current  in  the  right-hand  rectangular  waveguide  oi  Fig.  2  is 
given  by 

fS'  -  -R*  (u!”)'  ■  mil'”)  i»}  («is) 

where  is  pven  by  eq.  (2.12)  of  [1]  whidi  is  expreesed  in  abloeviated 
Ibnn  aa 


E  ViM4 

(aKl4'l 


where,  aa  in  the  main  program. 


K2b2*K1. 


(4.14) 


(4.15) 


Substitutioa  of  (4.14)  into  (4.13)  giva 


E  >r 


(4.16) 


where  is,  with  reference  to  eq.  (2.27)  of  [1],  the  (i— Kl,  i— Kl)*^  element 
of  the  matrix 


yrjTMJTM  ytJTMJTE  ’ 
ySJTEJTM  y/^JTEjTE 


The  above  matrix  is  diagonal.  Equatioa  (4.16)  is  recast  as 


-u.j_g(-i,yr)i^r). 


(4.17) 


Numerically,  (4-16)  should  be  1  to  P^*^  of  (4.8). 
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4.4  Reduction  of  the  C’s  of  (4.3)~(4.6)  to 
Those  of  [2]  When  RJf*  =  1  and  All 
the  Other  R’s  Are  Zero 

When  »  1  and  all  the  other  S’s  are  zero,  the  excitation  emplojped 

in  the  present  report  is,  as  stated  in  Section  3.5,  that  in  [2]  nmltiplied  by 
In  this  case,  the  V’s  in  (4.3H4.6)  •«  those  of  (21  mul¬ 
tiplied  by  In  view  of  (3.27),  this  multiplication  factor  is 

.  It  Mows  that  the  (Tn  of  (4.3)-K4.6)  ate  equal  to  those 
of  eqs.  (5.32),  (5.33),  (5.36),  and  (5.37)  of  [2]  when  =  1  and  all  the 
other  B’s  are  equal  to  zero. 
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Chapter  5 

The  Propagating  Modes  in  the 
Circular  Waveguide 

If  several  modes  propagate  in  the  drcular  waveguide,  then  eq.  (6.73)  of  [2] 
has  to  be  changed  because  this  equation  was  written  assuming  that  only  the 
even  TMoi  and  even  and  odd  TEn  modes  propagate  and  that  the  excitation 
is  the  even  TMoi  wave  whose  d^ric  and  magnetic  fields  are,  respectivdy, 
and  given  by  eqs.  (5.1)  and  (5.2)  of  [!].♦ 

5.1  Derivation  of  an  Expression  for  the 
Electric  Field  of  the  Propagating 
Modes 

From  eqs.  (B.26),  (B.37),  (B.30).  and  (B.33)  of  [1]  for  odd  TM  modes  and 
from  eqs.  (B.55),  (B.56),  (B.59),  and  (B.62)  of  [1]  for  odd  TE  modes,  we  have 

5s1,2,.-.  (5.1) 

Actually,  (5.1)  defines  the  above  because  these  £m’s  were  undefined 

before  (5.1)  was  written.  From  eqs.  (3.57)  and  (3.58)  of  [2], 

X,,  <  ka  (5.2) 

left-haiKl  side  of  «q  (51)  of  (1]  slweld  b*  In  u  early  version  of  [1],  it 

was  nustakenly  writua  as  ■ 
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where 


(5^) 


0^a  *  yj{hay  -  **,  (5.4) 

(5.5) 

In  view  of  eqs.  (5.1)-(5.3),  substitution  of  eq.  (6.55)  of  [2]  and  (3.7)  into  eq. 
(6.1)  of  (2)  gives 
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(5.9) 


(5.10) 


CTfi»  _  rTB« 

fd - 


+  2^ 


n(“)(? 


The  first  doable  summation  in  (5.6)  is  over  all  nonnegative  integers  r  and  all 
positive  integers  s  such  that 

x„  <  fco.  (5.11) 


The  second  double  summation  m  (5.6)  is  over  all  positive  integers  r  and 
3  such  that  (5.11)  holds.  The  third  double  summation  in  (5.6)  is  over  all 
nonnegative  integers  r  and  ail  positive  integers  s  such  that 

x',  <  ka.  (5.12) 

The  fourth  double  summation  in  (5.6)  is  over  all  positive  integers  r  and  s 
such  that  (5.12)  holds.  The  evanescent  waves  in  (5.6)  are  tho^  of  all  the 
evanescent  modes  in  the  circular  waveguide.  These  modes  consist  of  all  the 
even  and  odd  TM„  modes  for  whidi  x„  >  ka  and  all  the  even  and  odd  TE,, 

modes  for  which  x^,  >  ka.  _ 

The  transverse  part  of  the  term  in  (5.6)  reduces  to  y 

when  =  1  and  z  -  Lz.  This  reduced  term  carries  unit  time-average 
power.  The  multiplier  of  in  (5.6)  was  chosen  such  that  the  part 

of  C^*  attributable  to  direct  reflection  of  the  incident  even  TMr, 
about  r  =  I3  is  -1  when  Bjl**  =  1.  Similar  remarks  apply  to  the 

and  terms  in  (5.6).  _ 

The  term  in  (5.6)  reduces  to  when  *  1  and 

z  =  I3.  This  reduced  term  carries  unit  time-average  power.  The  multiplier 
of  io  (5.6)  was  chosen  such  that  the  part  of  attributable  to 

direct  reflation  of  the  incident  even  TEr.  wave*  about  r  =  X3  is  -1  when 
B'^  s  1.  Similar  remarks  apply  to  the  JSj^  and  the  terms  in  (5.6). 


5.2  The  Coefficients  and 

In  this  section,  C™*  and  of  (5.7)  and  (5.8)  are  expressed  in  a  form 

tThe  electric  field  of  this  incident  wave  is  the  term  in  (5.6). 

*The  electric  field  of  this  incident  wave  is  the  **nn  in  (5.6). 
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suitable  for  calculation.  In  (5.7)  and  (5.8),  and  5^7*"  given  by  eqs. 
(6.56)  and  (6.57)  of  [2]; 

3  /.TM3\ 

5?^  =  -v5E(-irE{V)  E  ‘hCi4'“  (S.13) 

Tr«l  imO  \  ®  /  pmO 

P+*#0 


s?^  =  -v?E(-i)‘'*'’’Ef^)  E  (S.14) 

•yml  fmO  V  ^ 

S+«»0 


where,  because  of  the  additional  factors  y/^  and  (— l)**^  in  (5.13)  and  (5.14), 
Cl  is  given  not  by  (6.61)  of  [2]  but  by 


Cl  *  9 


(5.15) 


The  factor  of  y/ri  explicit  in  (5.13)  and  (5.14)  will  cancel  out  when  and 
5^"  are  substituted  into  (5.7)  and  (5.8).  Equations  (6.28)  and  (6.29)  of  [2] 
are 


<  =  #l'’c«(^) (5-16) 


Comparing  eqs.  (6.15)  and  (3.76)  of  [2],  we  see  that,  since  Zr*  <  ka, 


^  =  1G;,  5  =  TM,TE.  (5.18) 

Substitution  of  (5.18)  into  (5.13)  and  (5.14)  gives 

5™'  =  -f  EC-irEf^r  E  (5-19) 

-psl  *=0  psO 


5^  =  -f  E(-1)''*'’’E<?™  E  (5-20) 

nnl  9*0  psO 

P+^ 
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Equation  (5.19)  is  witniW  to  eq.  (6.86)  [2].^ 

Substituting  (3.27)  and  (5.19)  into  (5.7),  we  obtain,  upon  uae  ci  (5.15), 

(5.21) 

where,  in  view  of  eq.  (6.60)  of  (2l  and  (3.110),  and  are 

given  by  the  right-hand  sides  of  (3.57),  (3.58),  (3.63),  and  (3.64)  with  (m,n,t) 
replaced  by  (p,q,i).  Substitution  of  (3.27)  and  (5.20)  into  (5.8)  fpves,  upon 
use  of  (5.15), 

(5.22) 

where,  in  view  of  (3.111),  and  are  given  by  the  right-hand 

sides  of  (3.65),  (3.66),  (3.69),  and  (3.70)  with  (m,n,t)  replaced  by  (p,9>i)* 

5.3  The  Coefficients  CJJ*  and 

In  this  section,  and  of  (5.9)  and  (5.10)  are  expressed  in  a  form 
suitable  for  calculation.  In  (5.9)  and  (5.10),  5^  and  5^  are  given  by  eqs. 
(6.58)  and  (6.59)  of  [2]  with  replaced  by 


“  **  pmO 

H-f#0 


Hv) 


^There  is  an  error  in  eq.  (6.86)  of  (2].  The  righuband  side  of  this  equation  should  be 
nniltiplied  by  —1/2. 
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E(-irE  E 

■7«1  fmC 

rH#0 


where,  because  of  the  additional  factors  y/rj  and  (-*1)'^  in  (5.23)  and  (5.24), 
Cl  and  C%  are  given  not  by  eqs.  (6.61)  and  (6.62)  of  [2]  but  by  (5.15)  and 

Equations  (6.30)  and  (6.31)  of  [2]  are 

(5^) 

CO.  (i)+,>{»«n(^). 

Comparing  eqs.  (6.16)  and  (3.108)  of  [2j,  we  see  that 


i— = 

c  2  * 


Substitution  of  (5.18)  and  (5.28)  into  (5.23)  and  (5.24)  gives 


51^  = 


p+«to 


(5.25) 

(5.26) 

(5.27) 

(5.28) 


.{r6r«?c.+  (2^)(^5^^)c,)  (5.29) 

2wX^^  ***1  f*0  fmO 


j-r6«<c.  +  (2^)  c«)  •  (5 :») 


65 


Equations  (5.29)  and  (5.30)  are  generalizatimis  oi  eqs.  (6.93)  and  (6.94)  c£ 
(2l.t 

Substituting  (3.30)  and  (5.29)  into  (5.9),  ere  obtain,  upon  use  oi  (5.15), 
(5.25),  and  eq.  (6.60)  of  [2], 


(5.31) 

where  is  given  by  the  right-hand  sides  of  (3.71)  and  (3.79)  with 

(m,n,t)  replaced  by  (p,9,i)  and  is  given  by  the  right-hand  sides 

of  (3.72)«d  (3.80)  with  (m,n,i)  replaced  by  (p,g,i).  On  th^  rig^t-hand 
sides,  and  are,  before  replacement  of  (m,n,i)  by  (p,9,i),  given 
by  (3.112)  and  (3.il3). 

Substitution  of  (3.30)  and  (5.30)  into  (5.10)  ipves,  upon  use  of  (5.15)  and 
(5.25), 


r+f»*o 

(5.32) 


In  (5.32),  is  given  by  the  right-hand  sides  of  (3.81)  and  (3.86)  with 

(m,n,*)  replaced  by  (p,9,i)  and  is  given  by  the  right-hand  sides 

of  (3.82Wnd  (3.87)  with  (m,n,»)  replaced  by  (p,g,j).  On  these  ri^t-hand 
sides,  and  7^*  are,  before  replacement  of  (m,n,i)  by  (p,g,i),  given 
by  (3.114)  and  (3.il5). 


^Equation  (6.94)  of  (2]  is  not  correct.  Please  correct  eq.  (6.94)  of  [2]  by  replacing 
therein  by 
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5.4  The  Time- Average  Power  of  the 

Propagating  Modes  Produced  by  the 
Magnetic  Current  in  the  Circular 
Waveguide 


The  electric  field  in  the  circular  waveguide  of  Fig.  2  is  produced  by  the  sources 
And  therein.  The  partial  electric  field  produced  by  the 

magnetic  currents  and  — is  pven  by 


+EE 


•fevanescent  waves. 


(5.33) 


Equation  (5.33)  was  obtained  by  removing  from  of  (5.6)  the  electric  field 
produced  by  The  latter  electric  field  consists  of  thetiouble  summations 
on  the  right-hand  side  of  (5.6)  with  C™*,  C™®,  and  replaced 
by  —5™*,  and  respectively.^  As  in  (5.6),  the  first 

double  summation  in  (5.33)  is  over  all  nonnegative  integers  r  and  all  positive 
integers  s  such  that  (5.11)  bolds.  The  second  double  summation  in  (5.33)  is 
over  all  positive  integers  r  and  s  such  that  (5.11)  holds.  The  third  double 
summation  in  (5.33)  is  over  all  nonnegative  integers  r  and  all  positive  integers 
3  such  that  (5.12)  holds.  The  fourth  double  summation  in  (5.33)  is  over  all 
positive  integers  r  and  s  such  that  (5.12)  holds. 

»If  and  were  zero,  the  S7.^\  sj^,  and  Sj®*  terms  in  (5.7)- 

(5.10)  would  be  zero  so  that  the  C»  would  reduce  to  the  negatives  of  the  B’». 


67 


The  time-Average  power  of  <rf  (5.33)  it  Prm  giw 

by 

ft-  - 

raO*Bl  rmXaml 

+ZE|cS**+ftir'r+EE|c7*+ft?'r  •  (*•«) 

fvO«Bl  ml«Bl 

The  quantity  Prm  of  (5.34)  it  the  time-average  power  of  the  propagating 
modet  produced  by  the  magnetic  currentt  and  in  the  circular 

waveguide.  In  brief,  Prm  of  (5.34)  it  the  time-average  power  radiated  hy  the 
magnetic  current  in  the  circular  waveguide.^ 


5.5  The  Time- Average  Power  Required  to 
Maintain  the  Magnetic  Current  in  the 
Circular  Waveguide 

The  time-average  power  required  to  maintain  the  magnetic  current  in  the 
circular  waveguide  is  Prma  by 

ft-.  =  -Re{/jfil-fl(S,JO<i»}  (5.35) 

where  M.  represents  the  combination  of  the  magnetic  currents  and 

— in  the  circular  waveguide,  K{Q.^M.)  is  the  magnetic  field  due  to  this 
combination,  and  5  represents  the  surfaces  on  which  M.  exists.  Now,  Prma  can 
be  viewed  as  an  alternative  expression  for  the  time-average  power  radiated 
by  the  magnetic  current  in  the  circular  waveguide.^  With  and 
given  by  (4.10)  and  (4.14),  (5.35)  becomes 

P,«.  =  Re{V-r3v}  (5.36) 

where  V  is  the  column  vector  of  the  V's  in  (4.10)  and  (4.14)  and  is  the 
square  matrix  given  by  eq.  (2.29)  of  [1]. 

^The  subscript  “rm"  oo  Prm  stands  for  "ladiated  by  the  magnetic  current”. 

^The  subscript  ”a"  on  Prma  stands  for  i^temative. 
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(2.22)  of  [1]  is  written  ccmdsriy  ss 

(i'*  +  y’+K*iv  =  /. 


FVom  (5.37), 


(5.37) 


y>P-/*-[y'*  +  y^(>.  (5.38) 

Piemultiplying  (5.38)  by  and  taking  the  leal  part,  we  haw 

Re{v-l'»<>}-Ile{«~/}-Re{<^(y>  +  K«]i>}.  (SJ9) 

Now, 

-  «  N1  N1 

V(y'+y«]V-j:inV5P  +  El^.i«nlK+m|’  (5.40) 

•al  tal 

where  and  the  same  as  in  (4.U)  and  (4.15).  Taking  the 

real  part  of  (5.40)  and  using  (4.11)  and  (4.15),  we  obtain 

Re{('(r‘  +  K>l<^}  =  j4'l  +  i4»>.  (5.41) 

Substituting  (5.41)  into  (5.39)  and  using  (5.36),  we  obtain 

=  {^/}  -  Pi’>  -  pS^  (5.42) 

which  is  recast  as 

=  -Im I [-jy/iiJ]  J  -  ^  -  pS^  .  (5.43) 

Numerically,  Prm*  of  (5.43)  should  be  equal  to  P^m  of  (5.34). 


5.6  The  Time- Average  Power  of  the 
Propagating  Modes  in  the  Circular 
Waveguide 


The  time-average  power  of  the  r-traveling  (incident)  waves  in  expression  (5.6) 
for  the  electric  field  in  the  drcular  waveguide  is  P^  given  by 

- -  rml  tml 


rmOtml 


+EE|B^r+EE|BS^r- 


rsOfail 


rml  tm\ 


(5.44) 


The  subscript  <»  stands  for  *mddeiit*.  The  time-average  povrer 
of  the  — s-travding  waves  in  expression  (5.6)  is  P,  given  by 

Pr  -  EE|c~r+EEtc~r 

+EE|cJ^r+EE|c"*r .  (S.4S) 

r^«al  ralMl 

The  subscript  "r”  on  P,  stands  tor  ^reflected*.  Numerically,  the  time-average 
power  of  the  x-traveling  waves  in  the  circular  wav^ui^  must  equal  the 
sum  of  the  time-averaM  power  P,  reBected  in  the  circular  waveguide  and  the 
time-average  power  i|  ^  transmitted  into  the  rectangular  waveguides: 

P*  *  Pr  +  Pf  **  +  P/*’ .  (5.46) 

According  to  (5.46),  part  of  P^c  is  reflected  in  the  circular  waveguide,  and 
the  rest  of  P^e  is  transmitted  into  the  rectangular  waveguides. 

5.7  Comparison  of  the  Coefficients 
Cif*,  and  C7f®  with  Those  of  [2] 

The  excitation  in  (2)  is  the  r*traveling  wave  whose  electric  field  is 
given  by  eq.  (5.1)  of  (!].♦  This  excitation  gives  rise  to  the  electric  field 
in  the  circular  waveguide.  Assuming,  as  in  [2],  that  only  the  TMoi  and  TEn 
modes  can  propagate  in  the  circular  waveguide,  eq.  (6.75)  of  [2]  gives 


4>e\'ancscent  waves.  (5.47) 

The  right-hand  side  of  (5.47)  is  the  electric  field  that  would  exist  in  the 
circular  waveguide  if  the  exciutimi  in  (2]  were  the  r-traveling  wave  whose 
electric  field  is  . 

^The  excitation  in  m  tbe  aaiiw  m  that  in  (I).  The  left-hand  side  of  eq.  (5.1)  of  [1] 
should  be  la  an  early  vmion  of  (1).  it  was  mistakenly  written  as 
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s  1  it  the  <mly  nonzero  B  in  (3.1),  then  the  exdtntkHi  would  be 
the  z-traveling  wave  whooe  electric  field  iz  ^ 

as  in  the  previous  sentence.  If  only  the  TMn  and  TEn  mo^  can  propagate 
in  the  circular  waveguide,  then  this  exdtatioo  would  give  rise  to  the  electric 
field  of  (5.6)  in  the  circular  waveguide: 


+  evanescent  waves.  (5.48) 

The  right-hand  sides  of  (5.47)  and  (5.48)  are  equal  to  each  other  because 
th^  are  both  the  dectric  fidd  produced  by  the  same  exdtation — namehr,  the 
z-traveling  wave  whose  dectric  fidd  is  Therefore, 

the  C’s  in  (5.48)  are  equal  to  those  in  (5.47).  This  means  that  if  »  1 
is  the  only  nonzero  B  in  (3.1)  and  if  only  the  TMoi  and  TEn  mo^  can 
propagate  in  the  drcular  waveguide,  then  the  codfidents  C^**,  C^,  and 
in  (5.6)  are,  respectivdy,  equal  to  those  in  eq.  (6.75)  of  [2].  In  Appendix 
B,  we  show  that  expressions  (5.21),  (5.31),  and  (5.32)  for  and 

reduce  to  the  right-hand  sides  of  eqs.  (6.88),  (6.96),  and  (6.99)  of  [2] 
when  B^**  s  1  is  the  only  nonzero  B  in  (3.1)  and  when  only  the  TMqi  and 
TEii  modes  can  propagate  in  the  drcular  waveguide. 
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Chapter  6 

The  Tangential  Electric  Field 
in  the  Apertures 


In  eqs.  (3.328)«  (3.332),  (3.336),  and  (3.339)  of  [3]  for  the  nonnelised  ^  and 
z-components  of  the  dectric  field  in  tlw  ^>ertuxes,  is  the  electric  fidd 

of  the  Z'travding  wave  that  would  exist  in  the  circular  waveguide  if  the  aper* 
tures  were  dosed.  The  subscript  *nns*  on  in  eqfi.  (3.328),  (3.3^), 

(3.336),  and  (3.339)  of  [3]  denotes  the  root  mean  square  value  of  the  trans* 
verse  part  over  the  waveguide  cross  seetkm  at  z  s  0.  In  the  present  report, 
the  electric  field  of  the  z>travding  waves  that  would  exist  if  the  apertures 
were  closed  is  not  £#?***  but  £***^(^^,11)  given  by  (3.1): 


(6.1) 


Hence,  _ should  be  replaced  by  j£****(»r^,ll)j^  in  eqs.  (3.328), 

(3.332),  (3.336)7Iod  (3.339)  of  [3]. 

By  definition, 

|£(3+)(^  jj)|^  - 
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where  the  superscript  denotes  the  cmuplex  conjugate  and  the  subscript 
'V  denotes  the  transverse  part.  On  the  right-hand  side  of  (6.2),  stands 
for  Substituting  (6.1)  into  (6.2)  and  using  eqs.  (B.l),  (B.26), 

(B.35),  (B.55),  and  (B.64)  [1],  we  obtain 

(6.3) 


Substitution  of  (3.27)  and  (3.30)  into  (6.3)  gives 


From  eq.  (7.11)  of  [2]  and  (6.4),  we  have 


Multipiying  both  sides  of  eqs.  (3.328),  (3.332),  (3.336),  and  (3.339)  of  [3]  by 
(6.5),  we  have^ 


|£(«)(_^  J)|^ 


tTbe  left^hmad  akks  of  tei  (3^).  (3433).  (3436),  ud  (3.339)  of  [3]  are  incorrect. 
The  correct  left-hand  sidai  are  tJw  aegatieei  of  those  written  in  these  equations  (the  left- 
hand  sides  of  eqs.  (7.13)-<7  )4)  of  (3)  are  also  lacorrsct.  The  correct  left-hand  sides  are 
the  negatives  of  thiMe  smtten  w  thw  eqvatioas).  The  corrected  eqs.  (3.328),  (3.332), 
(3.336),  and  (3.339)  of  [3]  were  mohiplird  by  (6.3)  to  obtain  (6.6)-(6.9). 
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|£'«"(^.aT  *  5  5  vTivj 

r+«#o 

•{(f)(=f)-(?)(^))-(¥)“«f'« 

*■  S  S  V  4  ii„6j 


where  5«  is  given  not  by  eq.  (3.329)  of  [3]  but  by 

-•"\P 

5.. 


EE 


+E 


/X  «? 

/ 


r 


(6.10) 


The  right-hand  sides  of  (6.6)-(6.9)  are  those  of  eqs.  (3.328),  (3.332),  (3.336), 
and  (3.339)  of  [3]  with  Sa.  of  eq.  (3.329)  of  [3]  replac^  by  5«  of  (6.10)  and 


VlTM 

-22 _ _ _ 

ylTE 

V 

V 

VtTM^0Sf‘La 
_ES _ 

y^TM 

yTTE 

*? 

V 

(6.11) 

(6.12) 
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where  means  ‘‘replaced  by*. 

Throughout  this  paragraph,  it  is  assumed  that  »  1  and  all  the 
other  B*s  in  (6.10)  are  aero.  With  such  ^’s,  (6.10)  reduces  to 


5.* 


(6.13) 


The  5.  of  (6.13)  is  that  of  eq.  (3.329)  <d  [3]  multiplied  by  ~-^0n^/k 
FVirthennore,  according  to  (6.1),  the  excitation  in  the  present  report  reduces 
to  the  product  of  with  the  excitation  in  [3].^  Hence,  in  view 

of  (3.27),  the  V’s  in  (6.6)-(6.9)  reduce  to  the  products  yjkj 
with  those  in  eqs.  (3.328),  (3.332),  (3.336),  and  (3.339)  of  [3].  As  a  result, 
the  numerical  values  of  the  right-hand  sides  of  (6.6)-(6.9)  reduce  to  those  ci 
eqs.  (3.328),  (3.332),  (3.336),  and  (3.339)  of  [3]  multiplied  by 


^The  excitation  in  [3]  is  the  same  as  that  in  [1] — namely,  Eg 
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Chapter  7 

The  Main  Program 


The  main  program  of  the  present  report  is  a  modification  of  the  main  pro¬ 
gram  described  in  Sections  3.1  to  3.10  o'  [3]  and  listed  in  Section  3.11  of  [3]. 
Henceforth  in  Chapter  7,  the  main  program  of  the  present  report  will,  for  the 
sake  of  brevity,  be  called  the  main  pn^am.  In  Sections  7.1  to  7.7,  only  those 
main  program  statements  that  are  different  from  statements  in  the  main  pro¬ 
gram  of  Chapter  3  of  (3]  are  described.  These  main  program  statements  are 
described  by  defining  the  computer  program  variables  in  them  in  terms  of 
variables  in  Chapten  1  to  6.  In  Section  7.8,  the  entire  main  program  is  listed. 
The  FORTRAN  statements  described  in  Sections  7.1  to  7.7  are  included  in 
this  listing. 


7.1  Specification  Statements 

In  the  common  block  labeled  PHI.  PH1(150),  PH2(150),  PH3(150),  and 
PH4(150)  appear  in  the  main  program  as  opposed  to  PHl(lOO),  PH2(100), 
PH3(100),  and  PH4(100)  in  the  main  program  listed  in  Station  3.11  of  [3]. 
Since  this  common  block  also  appears  in  the  subroutine  PHI,  it  was  necessary 
to  replace  the  allocations  PHl(lOO).  PH2(100),  PH3(100),  and  PH4(100) 
in  the  subroutine  PHI  by  PHl(I50).  PH2(150),  PH3(150),  and  PH4(150). 
Except  for  this  small  change,  all  the  subprograms  called  by  the  main  pro¬ 
gram  are  the  same  as  those  in  |3].  The  variable  VI  and  the  arrays  CVTME, 
CVTEE,  and  C\^EO  appear  in  the  second  and  third  COMPLEX*8  speci¬ 
fication  statements  in  the  main  program  of  [3]  but  not  in  the  main  program. 
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In  the  third  R£AL*4  specification  statement,  ClOUTS  appears  instead 
of  ClOUTS,  and  PTOTAL  is  absent.  It  was  a  mistake  to  write  ClOUTS 
in  the  third  REAL*4  specification  statement  in  the  main  program  of  [3]; 
ClOUTS  should  be  there.  Howe^,  this  mistake  is  not  serious  because  it 
does  not  affect  any  numerical  results. 

The  fifth  R£AL*4  specification  statement  does  not  contain  the  variables 
PTMS,  PRECMS,  PTMR,  and  BKAG  present  in  the  fifth  R£AL«4  speci¬ 
fication  statement  in  the  main  program  of  [3].  It  was  a  mistake  to  specify 
the  variables  PTMS,  PRECMS,  and  PTMR  as  R£AL*4  in  [3]  because  these 
variables  are  not  used  therein. 

As  defined  in  the  statement  after  statement  155  in  the  main  program  of 
[3],  SA  is  complex.  Therefore,  the  ensuing  arrays  SINP  and  SINQ  and  vari¬ 
ables  SINQQ,  SINPP,  BMNJPP,  BMNJQQ,  BMNJPQ,  BMNJQP  should 
have  been  specified  complex  instead  of  real  in  the  main  program  of  [3]. 
This  erroneous  specification  will  generally  lead  to  incorrect  tangential  electric 
fields  in  the  apertures.  However,  the  tangential  electric  fields  listed  in  Section 
2.2.2  of  [3]  were  not  affected  by  this  erroneous  specification  because  the  value 
of  L3  was  chosen  such  that  ARC  s  x  in  the  statement  after  statement  155  in 
the  main  program  of  [3].  The  REAL*4  specification  of  the  arrays  SINP  and 
SINQ  and  variables  SINQQ,  SINPP,  BMNJPP,  BMNJQQ,  BMNJPQ,  and 
BMNJQP  is  correct  in  the  main  program  because  the  value  of  SA  calculated 
by  statement  155  in  the  main  program  is  purely  real. 

The  last  eight  specification  statements  in  the  main  program  are  new. 


7.2  The  Allowable  Range  of  ka 

As  for  the  control  statements  involving  BKA  early  in  DO  loop  48,  statement 
96  of  the  main  program  is  executed  only  if 

BKA  >  1.84118378  (7.1) 

but  statement  96  of  the  matin  program  of  [3]  is  executed  only  if 

2.40482556  <  BKA  <  3.05423693.  (7.2) 

Because 

i',,  »  1.84118378,  (7.3) 
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(7.1)  states  that  ka  must  be  lai:ge  enough  so  that  at  least  the  TEn  modes 
propagate  in  the  circular  waveguide.^  Here,  is  the  first  root  of  the  derivar 
tive  J{  of  the  Bessel  function  Jii  *u/«  is  the  cutoff  wavenumber  of  the  TEii 
circular  waveguide  modes.  Because 

XQt «  2.40482556  (7.4) 

^  3.05423693,  (7.5) 

(7.2)  states  that  ka  must  be  such  that  only  the  TEii  and  TMqi  modes  propa¬ 
gate  in  the  circular  waveguide.  Here,  xoi  is  the  first  root  oi  Jo,  and  is  the 
first  root  of  J^;  xoi/a  and  x^j/a  are,  respectively,  the  cutoff  wavenumbers  of 
the  TMoi  and  TE^i  circular  waveguide  modes. 


7.3  Calculation  of  the  Ps  Defined  in 
Section  3.4.6 

The  block  of  statements  beginning  with  statement  46  and  ending  with  state¬ 
ment  77  in  the  main  program  calculates  the  /*s  that  are  defined  by  (3.110)- 
(3.115).  This  block  of  statements  replaces  the  block  of  statements  beginning 
with  statement  46  and  ending  with  statement  80  in  the  main  program  of  [3]. 
The  latter  block  of  statements  calculates  the  excitation  vector  described  in 
Chapter  4  of  [2],  the  normalized  amplitude  of  the  wave  of  the  -z- traveling 
even  TMqi  circular  waveguide  mode  due  to  each  magnetic  current  expansion 
function  in  the  apertures,  and  the  normalized  amplitudes  of  the  waves  of 
the  — z-traveling  even  and  odd  TEn  circular  waveguide  modes  due  to  each 
magnetic  current  expansion  function  in  the  apertures. 

7.3.1  Preliminary  Statements 

The  block  of  11  statements  after  statement  12,  DO  loop  114,  and  the  three 
statements  before  DO  loop  114  define  variables  that  will  be  used  in  the  block 
of  statements  beginning  with  statement  46  and  ending  with  statement  77. 
The  ninth  and  tenth  statements  after  statement  12  set 

PMAX2  =  2(m„,„-|-l)  (7.6) 

^There  are  four  TEu  inodes:  the  ±2-traveliag  even  and  odd  ones. 
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where  rnmmK  is  the  maximuTn  value  of  m  (»  the  right>hand  sides  of  (3.110)- 
(3.115).  The  three  statements  before  DO  loop  114  set 


ARG  ~ 

(7.8) 

CSl  as  cos(r^«) 

(7.9) 

SNl  SB  sm(r^a) . 

(7.10) 

Here,  r,  which  appears  in  (3.110)-(3.115),  »  related  to  the  index  R  of  DO 
loop  19  by 

RsbT+I. 

(7.11) 

In  DO  loop  114, 

PHl(MJ)  =  #£'co.(r*.)  - 

(7.12) 

PH2(MJ)  =  *S'eo6(r*.)  +  *£'riii{'’*.) 

(7.13) 

PH3(MJ)  =  *S’cos(r*,)  -  *£> 

(7.14) 

PH4(MJ)  =  *£>co>(r*,)  +  #2’nnW.) 

(7.15) 

where 

MJ  =  M  +  PMAX2 

(7.16) 

and  m  is  related  to  M  by 

M  =  m  +  1. 

(7.17) 

Equations  (3.90),  (3.91),  and  (3.93)-(3.98)  reduce  (7.12)-(7.15)  to 

PHl(MJ)  = 

(7.18) 

PH2(MJ)  =  *£ 

(7.19) 

PH3(MJ)  =  ^2 

(7.20) 

PH4(MJ)  = 

(7.21) 
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7.3.2  Expressions  (3.110)  and  (3.111)  Associated 
with  the  TM  Propagating  Modes  in  the 
Circular  Waveguide 

The  block  of  stetements  extending  from  the  statement  after  statement  46  to 
statement  29  calculates  7^  of  (3.110)  and  of  (3.111).  This  block  of 
statements  will,  as  directed  by  control  statement  46,  be  executed  only  if  the 
TMr«  modes  propagate  in  the  circular  waveguide.  Here,  the  TMr«  modes  are 
the  ±2-traveling  even  and  odd  TMr«  nuxles^  for  the  fixed  value  r  related  to 
the  index  R  of  DO  loop  19  by  (7.11)  and  the  fixed  value  s  related  to  S  by 

S  «  s  (7.22) 


where  S  is  the  index  of  DO  loop  20.  In  the  previously  mentioned  block  of 
statements. 


o 

II 

(7.23) 

(7.24) 

OTM 

BTM(S,R)  =  ^ 
k 

(7.25) 

RTM  =  r+1 

(7.26) 

STM(R)  =  s. 

(7.27) 

In  (7.24), 

II 

1 

IV  II 
o 

• 

(7.28) 

The  value  of  TMS  in  (7.24)  is  its  value  inunediately  before  BTM(S,R)  is 
defined.  The  quantity  in  (7.25)  is  stored  in  BTM(S,R)  for  later  use 

in  calculating  the  curly  bracketed  quantity  in  expression  (6.10)  for  Sa-  The 
variables  RTM  and  STM(R)  will  be  used  only  after  exit  horn  DO  loop  19. 

^Even  TM,,  inodes  exist  for  r  =  0, 1, 2, .  ■ However,  odd  TM,, 

modes  exist  only  for 

*■=1,2,.... 
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The  vilues  of  RTM  and  STM(R)  after  exit  from  DO  loc^  19  are  given 
as  follows.  If  no  TM  mode  propagates  in  the  circular  waveguide,  then  the 
values  RTM  » 1  and  STM(l)  s  0  set  by  the  two  statements  afta  statement 
12  prevail.  However,  if  at  least  oae  TM  mode  prt^>agates,  then  nested  DO 
loops  19  and  20  set 


ItTM  =  I  ^ 
STM(R)  -  / 


the  maxiwniTn  value  of  R  for  whidi 
TMr.1,1  circular  waveguide  mode  pn>pagates 

the  maximum  value  of  S  for  which 
TMr_i^  circular  waveguide  mode  propagates 


I  (7.29) 
}.  (7.30) 


Equation  (7.30)  holds  for  R  =  1,2,  *  *  •  ,RTM. 
Upon  entry  into  DO  loop  52, 


where  n  is  related  to  the  index  N  of  outer  DO  loop  29  by 

N  =  n  + 1. 


(7.31) 


(7.32) 


In  (7.31),  n  appears  on  the  right-hand  sides  of  (3.110)  and  (3.111).  The  index 
M  of  DO  loop  52  is  given  by  (7.17)  in  which  m  appears  on  the  right-hand 
sides  of  (3.110)  and  (3.111).  After  being  increased  in  DO  loop  52, 


MTE  =  m  -I- 


n  =5  0 


^MM(0,  n  >  1 


(7.33) 


where  m  and  n  are  related  to  the  DO  loop  indices  M  and  N  by  (7.17)  and 
(7.32).  With  MTE  given  by  (7.33),  eq.  (3.9)  of  [3]  applies: 


BMN(MTE)  =  k^nb- 


(7.34) 


Inunediately  before  JTME  is  increased  inside  DO  loop  52,  TTM  is  given 


TTM  * 


(7.35) 
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where  7^  is  given  by  (3.60).  After  being  increased  in  DO  loop  52, 

+  «>1 

k  Iml 

where  ju/t  is  wbntever  JTME  was  upon  entry  into  DO  loc^  29.  The  ri^t- 
hand  side  of  (7.36)  is  jta  plus  the  right-hand  side  of  eq.  (3.147)  of  [3]  with 
M  and  N  given  by  (7.17)  and  (7.32).  In  DO  loop  52, 

ITME(JTME)s;J^  (7.37) 

where  is  pven  by  (3.110).  If  R  1,  JTMO  b,  after  bong  increased  in 

DO  loop  52,  given  by 

0,  n  »  0 

JTMO  =  i.,  +  m  +  .  (7.38) 

.  iml 

where  juo  b  whatever  JTMO  was  upon  entry  into  DO  loop  29.  If  R  1  in 
DO  loop  52, 

ITMO(  JTMO)  s  (7.39) 

where  “  given  by  (3.111). 

7.3.3  Expressions  (3.112)-(3.115)  Associated  with 
the  T£  Propagating  Modes  in  the*  Circular 
Waveguide 

The  block  of  statements  extending  from  the  statement  after  statement  65  to 
statement  77  calcubtes  of  (3.112),  of  (3.113),  of  (3.114), 
and  of  (3.115).  This  block  of  staten^ts  will,  as  directed  by  contrd 
statement  65,  be  executed  only  if  the  TE^*  modes  propagate  in  the  drcular 
waveguide.  Here,  the  TEr«  modes  are  the  ±z-traveling  even  and  odd  T£r« 
modes^  for  the  fixed  values  of  r  and  s  related  to  the  DO  loop  indices  R  and 

^Ev«n  TEp,  modes  exist  for  r  «  0, 1,2, •  **.  However,  odd  TEr*  modes  exist  only  for 
»■= 
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S  by  (7.11)  and  (7.22).  In  the  previously  mentioned  block  of  statements, 


GE- 

(7.40) 

(7.41) 

BTE{S,R)  =  ^ 

(7.42) 

(7.43) 

H  (*’‘)  ^))  ( * )  (“'*)  (e  J 

RTE*r+l 

(7.45) 

STE(R)  «  s. 

(7.46) 

The  value  of  TES  in  (7.41)  is  its  value  immediately  before  BTE(S,R)  is 
defined.  The  quantity  in  (7.42)  is  stored  in  BTE(S,R)  for  later  use 

in  calculating  the  curly  bracketed  quantity  in  expression  (6.10)  for  5..  The 
variables  RTE  and  STE(R)  will  be  used  only  aftCT  exit  from  DO  loop  19. 

After  exit  from  DO  loop  19, 

• 

_  r  the  maximum  value  of  R  for  which 

~  1  a  TEr.ij  circular  waveguide  mode  propagates 

1  (7.47) 

»Tp/o\  .  /  maximum  value  of  S  for  whidi  ] 

^  ^  a  TEr.i^  circular  waveguide  mode  propagates  j 

..  (7.4«) 

Equation  (7.48)  is  valid  for  R  w  2, 3,"*, RTE.  Because  of  (7.1),  the  TEu 

circular  waveguide  inodes  propagate  so  that  RTE  >  2.  If  the  TEn  circular 
waveguide  modes  propagate,  then  STE(l)  is  given  by  (7.48).  However,  if  the 
TEoi  circular  waveguide  modes  do  not  propagate,  then  the  value  STE(l)  s  0 
set  by  the  fourth  statement  after  statement  12  prevails  after  exit  from  DO 
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loop  19.  The  third  rtatenaent  after  statement  12  aaiires  that  RTE  ^  1  in  the 
unlikely  case  that  the  TEii  circular  wav^uide  modes  are  so  close  to  being 
cutoff  that,  due  to  floating  point  imuidoff  error,  ITM  s  2  when  the  indices 
R  and  S  of  nested  DO  loops  19  and  20  are  respectivdy  equal  to  2  and  1  so 
that  these  DO  loops  fail  to  define  RTE.^ 

Upon  entry  into  DO  loop  78, 


TTEPN 


(T.49) 


TTEZN 


<5i*'  (r.so) 


where  n  is  related  to  the  index  N  of  DO  loop  77  by  (7.32).  In  (7.49)  and 
(7.50),  n  appears  on  the  right-hand  ndes  of  (3.112)-(3.115).  The  index  M  of 
DO  loop  78  is  given  by  (7.17)  in  which  m  appears  on  the  right-hand  sides 
of  (3.112)-(3.115).  After  being  increased  in  DO  loop  78,  MTE  is  given  by 
(7.33).  With  MTE  given  by  (7  J3),  BMN(MTE)  is  given  by  (7.34). 

Upon  arrival  at  statement  71, 


TEPM 

TE2M 


where  is  given  by  (3.75)  After  statement  71  is 

0, 


JTEE  *  Jm  m  < 


n 


i;MM(0.  n 


Gi‘> 

executed, 
=  0 
>  1 


tel 


(7.51) 

(7.52) 


(7.53) 


^The  validity  of  (7.1)  asMtsi  Umi  tkc  TCu  areolar  woregnide  modes  i^opagate.  As 
calculated  in  the  sobrootiae  OGN  described  and  listed  in  Giapter  9  ot  (3],  ITEsl  when 
Rs2  and  Ssl  if  the  TEn  areolar  waveguide  modes  propagate  and  ITEsS  when  Rs2 
and  Ss  1  if  the  TEn  nreotar  waveguide  modes  do  not  propagate.  Therefore,  the  validity 
of  (7.1)  dioald  aasore  that  ITE  *  1  when  Rs2  and  S  w  1.  However,  if  (7.1)  were  valid  but 
ba^yso,  then  rooodoff  error  eoold  result  w  ITE *3  when  Rs2andSsl. 
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where  jn  is  whatever  JTEE  wu  upon  entry  into  DO  loop  77.  The  right- 
hand  ride  of  (7.53)  is  jn  plus  the  right-hand  ride  of  eq.  (3.147)  [3]  with 

M  and  N  given  by  (7.17)  and  (7.32).  In  DO  kx^  78, 

ITEEP(JTEE)  *  (7.54) 

ITEEZ(JTEE)  =  (7.55) 


where  and  are  given  by  (3.112)  and  (3.113).  If  R  1,  JTEO  is, 
after  being  increased  in  DO  loop  78,  given  by 


JTEO  *  jto  +  n»  + 


0, 

tml 


n  s  0 
n>  1 


(7.56) 


where  jto  is  whatever  JTEO  was  upon  entry  into  DO  loop  77.  If  R  /  1  in 
DO  loop  78, 

ITEOP(JTEO)  *  (7.57) 

ITEOZ(JTEO)  *  1^*  (7.58) 


where  and  7^^  are  given  by  (3.114)  and  (3.115). 


7.4  Calculation  of  the  Normalized 

Elements  of  the  Excitation  Vector 

The  normalized  elements  of  the  excitation  vector  are  given  by 

(3.13). 

7.4.1  The  Br«'s  in  (3.13) 

The  coefficients  and  Bj^  in  (3.13)  will  be  put  in  the 

computer  program  variables  BTME(s.r-H),  BTMO(s,r-|-l),  BTEE(s,r-H), 
and  BTEO(s,  r  -f- 1 ).  respectively.  Passing  through  the  block  of  statements 
beginning  with  statement  68  and  ending  with  statement  43,^  one  arrives  at 
nested  DO  loops  191  and  192.  These  nested  DO  loops  set  to  zero  all  the 

^This  block  of  ttatemeau  is  the  same  as  in  the  main  program  of  [3]. 
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BTM£*8  and  BTMO’s  asodated  with  the  TM  pn^Mgating  modes  in  the  dr- 
ciilar  waveguide.  Nested  DO  ktops  193  and  194  set  to  zero  ail  the  BTEIE’s 
and  BTEO’s  assodated  with  the  TE  propagating  modes  in  the  circular  wave¬ 
guide.  Actually,  the  BTM0(S,1)  and  the  BTEO(S,l)  that  could  be  set  to 
zero  for  some  values  of  S  are  not  coeffidents  in  (3.13).  We  did  not  bother  to 
avmd  executing  the  second  statement  in  DO  192  when  Rs  1  nor  did  we 
bother  to  avdd  executing  the  second  statement  in  DO  loop  194  when  Rs  1. 

The  BTME’s  and  BTMO’s  appear  in  (2.25).  Nonzero  values  of  them  are 
read  in  inside  DO  loop  156.  The  BTEE’s  and  BTEO’s  i^pear  in  (2.26). 
Nonzero  values  of  them  are  read  in  uiside  DO  loop  167. 

7.4.2  The  Contribution  of  the  Even  and  Odd  TM 

Propagating  Modes  in  the  Circular  Waveguide 

Nested  DO  loops  72  and  73  accumulate  in  XITM,  XlTE,  X2TM,  and  X2TE 
the  terms  of  the  summation  in  (3.13)  due  to  the  even  and  odd  TM  proper  ' 
gating  modes  in  the  drcular  waveguide.  After  exit  from  nested  DO  loops  72 
and  73, 


xiTM  = 


►TM€  f  lTM.TMe  ,  dTMo  rXTM.TMo 


+ 


0 


XlTE  =  23  (B^T**/**^™*  + 
X2TM  *  Y. 

X2TE  *  Y 


(7.59) 

(7.60) 

(7.61) 

(7.62) 


where  the  subscript  t  on  the  /,’s  denotes  dependence  on  the  integers  m  and 
n  in  (3.14).  The  eight  /,’s  on  the  right-hand  sides  of  (7.59)-(7.62)  are  given 
by  (3.57),  (3.58),  (3.63),  (3.64).  (3.65),  (3.66),  (3.69),  and  (3.70)  where  m 
and  n  are  related  to  the  indices  M  and  N  of  outer  DO  loops  74  and  16  by 
(7.17)  and  (7.32).  The  quantities  and  in  (3.57)  and  (3.65)  were 
previously  calculated.  See  (7.37)  and  (7.39).  Expressions  (7.59)-(7.62)  are 
the  contributions  of  all  the  even  and  odd  TM  circular  waveguide  propagating 
modes  to  the  summation  on  the  right-hand  side  of  (3.13)  with  afi  s  ITM, 
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ITE,  2TM,  and  2TE,  respectively.  The  calculation  of  these  expressions  is 
described  in  the  next  three  paragraphs. 

Both  MTM  and  MTE  are  set  to  sero  before  entry  into  DO  loop  16.  The 
integer  MTM  will  be  treated  in  Section  7.4.4.  After  being  increased  inside 
DO  loop  74,  MTE  is  given  by  (7.33).  The  indices  R  and  S  of  DO  loops  72  and 
73  are  related  to  the  summation  indices  r  and  s  in  (7.59)-(7.62)  by  (7.11) 
and  (7.22).  Because  JTME  and  JTMO  are  increased  by  KTE  of  eq.  (4.13) 
of  [3]  immediately  after  use  in  DO  loop  73, 

ITME(JTME)  *  (7.63) 

in  the  first  statement  in  DO  loop  73  and 

ITMO(JTMO)  =  (7.64) 

in  the  tenth  statement  in  DO  loop  73.  Consequently,  these  statements  set 

Z1  =  (7.65) 

and 


(7.66) 

respectively. 

The  third  and  fourth  statements  in  DO  loop  73  set 

W1  =  (-l)'+*n  (7.67) 

W2  *  (-1)'**  (7.68) 

so  that,  in  the  fifth  through  eighth  statements  in  DO  loop  73, 

-Wl.Zl  = 

-W2*Z1  =  ^TM.^lTE.TMe 
FNl.Zl  =  ^^.jrrM.TMe 
FMl.Zl  = 


(7.69) 

(7.70) 

(7.71) 

(7.72) 


where  the  /,’s  are  given  by  (3.57),  (3.58),  (3.63),  and  (3.64). 


The  lest  six  statements  in  DO  loop  73  are  not  executed  for  R  1  because 
there  are  no  odd  TMo«  modes  in  the  circular  waveguide.  In  the  lest  four 
statements  in  DO  loop  73, 


wivzi = 

(7.73) 

W2*Z1  * 

(7.74) 

# 

FN1*Z1  = 

(7.75) 

FM1*Z1  s 

(7.76) 

where  the  //s  are  pven  by  (3.65),  (3.66),  (3.69),  and  (3.70).  It  is  evident 
from  (7.69)-(7.76)  that  XITM,  XITE,  X2TM,  and  X2TE  will  be  given  by 
(7.59)-(7.62)  upon  exit  from  nested  DO  loops  72  and  73. 


7.4.3  The  Contribution  of  the  Even  and  Odd  TE 

Propagating  Modes  in  the  Circular  Waveguide 


Nested  DO  loops  168  and  169  add  to  what  was  already  stored  in  XITM, 
X1T£,  X2TM,  and  X2TE  the  terms  of  the  siimmation  in  (3.13)  due  to  the 
even  and  odd  TE  propagating  modes  in  the  circular  waveguide.  After  exit 
from  nested  DO  loops  168  and  169, 


XITM  * 

+  L 

+ E  + B,?*/?"-”*) 

X2TM  ss  5; 


(7.77) 


(7.78) 


(7.79) 


# 
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(7.80) 


+  E 

r,» 


where  the  subscript  i  on  the  /,’s  denotes  dependence  on  the  integers  m  and  n 
in  (3.14).  There  are  two  summations  on  the  right-hand  side  of  eadi  of  (7.77)- 
(7.80).  The  first  summation  is  what  was  already  stored  (see  (7.59)-(7.62)). 
The  calculation  of  the  second  summation  is  discussed  in  this  section.  The 
eight  /,*s  in  the  second  stunmations  on  the  right-hand  sides  of  (7.77)-(7.80) 
are  given  by  (3.71),  (3.72),  (3.79),  (3.80),  (3.81),  (3.82),  (3.86),  and  (3.87) 
where  m  and  n  are  related  to  the  indices  M  and  N  of  outer  DO  loops  74 
tnd  16  by  (7.17)  ud  (7.32).  The  qeeetities  ud  jTSf* 

in  (3.71)  and  (3.81)  were  previously  computed.  See  (7.54),  (7.55),  (7.57), 
and  (7.58).  The  second  summations  in  (7.77)-(7.80)  are  the  contributions 
of  all  the  even  and  odd  TE  circular  waveguide  propagating  modes  to  the 
summation  on  the  right-hand  side  of  (3.13)  with  a/J  s  ITM,  ITE,  2TM,  and 
2TE,  respectively.  The  calculation  of  these  second  sununations  is  described 
in  the  next  two  paragraphs. 

The  indices  R  and  S  of  DO  loops  168  and  169  are  related  to  the  summation 
indices  r  and  s  in  (7.77)-(7.80)  by  (7.11)  and  (7.22).  Because  JTEE  is 
increased  by  KTE  of  eq.  (4.13)  of  [3]  immediately  after  use  in  DO  loop  169, 
the  first  and  second  statements  in  DO  loop  169  set 


W5  = 

(7.81) 

W6  =  . 

(7.82) 

tements  in  DO  loop  169  set 

W3  =  /f™ 

(7.83) 

W4  = 

(7.84) 

W1  =  //TM-TEe 

(7.85) 

VV2  = 

(7.86) 

Z1  = 

(7.87) 
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where  the  /,’s  are  given  by  (3.71),  (3.72),  (3.79),  and  (3.80).  In  the  ninth 
through  twelfth  statements  in  DO  loop  169, 


W1*Z1  = 

(7.88) 

W2*Z1  * 

(7.89) 

W3*Z1  = 

(7.90) 

W4*Z1  = 

(7.91) 

The  last  twelve  statements  in  DO  loop  169  are 

not  executed  if  R  =  1 

because  there  are  no  odd  TEo,  modes  in  the  drcular 

waveguide.  Because  the 

third  of  these  statements  increases  JTEO  by  KTE  of  eq.  (4.13)  of  [3],  the 

first  two  of  these  statements  set 

ws-/?jr» 

(7.92) 

W6- 

(7.93) 

The  fourth  through  eighth  of  these  statements  set 

W3  s  _^rTM,TB» 

(7.94) 

W4  * 

(7.95) 

W1  * 

(7.96) 

^  ^IT£.TE* 

(7.97) 

Z1  *  B^ 

(7.98) 

where  the  A’s  arc  given  by  (3.81),  (3.82),  (3.86),  and  (3.87).  In  the  last  four 

statements  in  DO  loop  169, 

WNZl  « 

(7.99) 

W2.Z1  « 

(7.100) 

-W3.ZI  « 

(7.101) 

-W4.Z1  « 

(7.102) 

It  is  evident  from  (7.59)-(7.62).  (7.88)-(7.91),  and  (7.99)-(7.102)  that  XlTM, 
XITE,  X2TM,  and  X2TE  will  be  given  by  (7.77)-(7.80)  upon  exit  from 

nested  DO  loops  168  and  169. 
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7.4.4  Storage  of  the  Normalized  Elements  of  the 
Excitation  Vector 

The  eighth  horn  the  last  through  the  fifth  from  the  last  statements  in  DO 
loop  74  are  not  executed  if  either  N  s  1  or  M  s  1  because  there  is  neither 
a  TMon  nor  a  TM^o  rectangular  waveguide  mode  expansion  function.  The 
seventh  from  the  last  and  the  fifth  from  the  last  statements  in  DO  loop  74 
use  XITM  of  (7.77),  X2TM  of  (7.79),  and 


•»-» 

II 

(7.103) 

TI(MTM)  « 

(7.104) 

TI(Jl)  =  -jy/riir^. 

(7.105) 

The  right-hand  sides  of  (7.104)  and  (7.105)  are  given  by  (3.13)  with  afi  = 
ITM  and  2TM,  respectivdy.  These  are  the  normalized  elements  of  the  excita¬ 
tion  vector  corresponding  to  the  TM  rectangular  waveguide  mode  expansion 
functions.  In  (7.104),  MTM  is  given  in  terms  of  the  indices  N  and  M  of  DO 
loops  16  and  74  by  eq.  (3.146)  of  (3).  In  (7.105), 

J1  =  MTM-|-K1  (7.106) 

where  Kl  is  the  total  number  of  T£  and  TM  rectangular  waveguide  mode 
expansion  functions  given  by  eq.  (3.10)  of  [3]  where  KTE  and  KTM  are  given 
by  eqs.  (4.13)  and  (4.14)  of  (3]. 

The  third  from  the  last  and  the  last  statements  in  DO  loop  74  use  XlTE 
of  (7.78),  X2TE  of  (7.80),  and  (7.103)  to  set 

TIlJl)  e  (7.107) 

THJ2»  * (7.108) 

The  right-hand  sides  of  (7.107)  and  (7.108)  are  given  by  (3.13)  with  afi  — 
ITE  and  2TE,  respectively.  These  are  the  normalized  elements  of  the  excita¬ 
tion  vector  corresponding  to  the  TE  rectangular  wav^ide  mode  expansion 
functions.  In  (7.107). 

JlsMTE-J-KTM  (7.109) 
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where  MTE  is  given  in  terms  of  n  and  m  by  (7.33).  Alternatively,  MTE  is 
given  in  terms  of  the  indices  N  and  M  of  DO  loops  16  and  74  by  eq.  (3.147) 
of  [3].  In  (7.108), 

J2*MTE  +  KTM  +  K1.  (7.110) 

7.5  The  Normalized  Coefficients  of  the 
TEio  Modes  in  the  Rectangular 
Waveguides 

The  normalized  coefficients  of  the  waves  of  the  TEio  modes  in  the  rectangular 
waveguides  are  given  by  (4.3)->(4.6)  as  opposed  to  those  of  [3],  which  are  given 
by  eqs.  (5.32),  (5.33),  (5.36),  and  (5.37)  of  [2].^  Equations  (4.3)  and  (4.4) 
differ  firom  eqs.  (5.32)  and  (5.33)  of  [2]  only  in  that  the  right-hand  sides  of 
(4.3)  and  (4.4)  have  the  multiplier  yjkwftt  while  the  right-hand 

sides  of  (5.32)  and  (5.33)  of  (21  have  the  mulUplier 

Equations  (4.5)  and  (4.6)  differ  from  eqs.  (5.36)  and  (5.37)  <>f  [2)  only  in  that 

the  right-hand  sides  of  (4.5)  and  (4.6)  have  the  multiplier 

while  the  right-hand  tides  of  cqt.  (5.36)  and  (5.37)  of  [2]  have  the  multiplier 

In  the  main  program  of  (3}.  the  nght-hand  sides  of  eqs.  (5.32),  (5.33), 
(5.36),  and  (5.37)  were  calculated  with  fri  and 

stored  in  the  array  W  In  the  main  program  (of  the  present  report),  the 
matrix  equation  (2.22)  of  (1],  cooaseiy  written  as 


rf  .  r 

(7.111) 

Y\  m  TI 

(7.112) 

where  Y,  V,  and  TI  are  arrays  i«  the  main  program.  The  array  Y  rantains 
the  elements  of  The  array  TI  contains  the  elements  of  From 

^There  is  a  mispnal  «  (5  sf  !7]  Tlw  Itfk-baad  side  of  (5.37)  of  [2]  should  be 

/-•JTE— 

0x0  • 
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(7.111),  (7.112),  and  the  definitions  of  Y  and  T1  in  (7.112),  one  sees  that  the 
array  V  contains  the  dements  of  Vly/rj.  Two  of  these  elements  are 
«nd  V^/v^.  Hence,  and  in  (4.3)-(4.6)  are  treated  in 

the  main  program  exactly  the  same  as  and  in 

eqs.  (5.32),  (5.33),  (5.36),  and  (5.37)  of  [2]  were  treat^  in  the  main  program 
of  [3]. 

The  block  of  statements  beginning  with  the  statement  after  statement  16 
and  ending  with  statement  70  writes  Tl,  sol^  the  matrix  equation  (7.112), 
writes  the  solution  V,  calculates  the  normalised  coefficients  (4.3)-(4.6)  ci  the 
waves  of  the  TEio  modes  in  the  rectangular  waveguides,  and  writes  these 
normalized  coefficients.  In  this  block  of  statements,  only  the  statement  that 
defines  SI  is  different  from  that  in  the  main  program  of  [3].  In  the  main 
program, 

SI  =  (T.U3) 

as  opposed  to 

SI  =  (7.IU) 

in  the  main  program  of  (3].  The  Si  of  (7.113)  is  the  product  of  1/2  with 
the  factor  in  (4.3)”(4.6)  while  the  Si  of  (7.114)  is  the  product  of  1/2 

with  the  factor  ^Ao/dlo?*  (5.32),  (5.33),  (5.36),  and  (5.37)  of  [2].  In 
the  main  program.  ClOUT,  ClIN,  C20UT,  and  C2IN  are  Cjo®”  of  (4.3), 
Cio^  of  (4.4),  of  (4,5),  and  of  (4.6),  respectively.  These  are 

the  normalized  coefficients  of  the  waves  of  the  TEio  modes  in  the  rectangular 
waveguides. 

The  seven  statements  following  statement  70  are  the  same  as  the  seven 
statements  following  statement  81  in  the  main  program  of  [3].  These  seven 
statements  calculate  and  wiite  the  squares  of  the  magnitudes  of  the  normal¬ 
ized  coefficients  of  the  waves  of  the  TEio  inodes  in  the  rectangular  waveguides 
and  the  time-average  power  transmitted  into  the  rectangular  wave¬ 

guides.  Here,  P/*’  and  are  given  in  terms  of  these  squares  by  (4.7)  and 
(4.8),  respectivdy.  The  fifth  statement  after  statement  70  sets 

PT*s/^**  +  P/’*. 


(7.115) 


7.6  The  Normalized  Coefficients  of  the 
Propagating  Modes  in  the  Circular 
Waveguide 


The  block  of  statements  beginning  with  the  statement  after  statement  82 
and  ending  with  statement  190  calculates  and  writes  Cf^  of  (5.21), 
of  (5.22),  of  (5.31),  Cf^  of  (5.32),  and  the  time-average  powers  Prm 


of  (5.34),  ci  (5.44),  and  Pr  of  (5.45).  This  blodc  of  statements  also 
calculates  the  computer  program  variable  W7  defined  by 


wt-L  E 


+EE 


(7.116) 


The  above  expression  spears  in  (6.10).  Now,  Cj?**,  and 

are  the  normalised  coefficients  of  the  waves  Of  the  -s-traveling  even  TMr« , 
odd  TMr« ,  even  T£,« ,  and  odd  TE.«  circular  waveguide  modes,  respectivdy. 
The  quantities  Prmi  Pme%  *nd  Pr  depend  on  these  normalized  co^dents. 
Here,  Prm  i*  the  time-average  power  that  the  combination  of  magnetic  cur¬ 
rents  — and  — in  Fig.  2  would  radiate  in  the  absence  of  P*e 
is  the  time-average  power  of  the  r-traveling  waves  in  the  circular  waveguide, 
and  Pr  is  the  time-average  power  of  the  — z-traveling  waves  in  the  drcular 
waveguide.  The  right-hand  side  of  (7.116)  is  the  term  in  brackets  in  the 
denominator  of  (6.10).  In  this  term,  the  first  double  summation  is  over  all 
the  even  and  odd  TM..  propagating  modes  in  the  drcular  waveguide  and 
the  second  double  summation  is  over  all  the  even  and  odd  TEr.  propagating 
modes  in  the  drcular  waveguide. 


7.6.1  The  Normalized  Coefficients  of  the  TM 

Propagating  Modes  in  the  Circular  Waveguide 

The  normalized  coeffidents  and  of  the  waves  of  the  even  and  odd 

TM,«  propagating  modes  in  the  drcular  waveguide  are  calculated  in  nested 
DO  loops  172  and  173  whose  indices  R  and  S  obtain  r  and  s  according  to 
(7.11)  ’  (7.22).  These  DO  loops  also  add  to  W7  the  TM;,  and  TM;.  terms 

in  (7.U(  >  'pedfically,  the  fifth  statement  in  DO  loop  173  adds  (iB^^P  -i- 
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{\B'^^)0^lk  to  W7.  Inside  nested  DO  loops  172  and  173,  nested  DO 
loops  174  and  176  accumulate  the  triple  summations  (5.21)  and  (5.22)  in 
XTME  and  XTMO,  respectively.  After  exit  from  nested  DO  loops  174  and 


176, 

XTME  *5:5; 

-T-ir-o 

L 

1  y7TM.TM«  1 

1  ^/n  J 

|+77”*™'| 

1  (7.117) 

XTMO  *  i;  z 

L 

pmO 

[  y/n  ] 

1  >aI 

J .  (7.118) 

The  calculation  of  the  right-hand  sides  of  (7.117)  and  (7.118)  is  described  in 
the  next  six  paragr^hs. 

Inside  nested  DO  loops  174  and  176,  the  pf*^  terms  in  (7.117)  and  (7.118) 
are  added  to  XTME  and  XTMO,  respectively,  where  q  and  p  are  related  to 
the  DO  loop  indices  N  and  M  by 


N«f  +  1 
M  «  p+  1. 

There  are  possibly  four  pf*^  terms  in  (7.117).  They  are 


(7.119) 

(7.120) 


(7.121) 


Preliminary  to  the  calculatioo  of  these  terms,  the  third  and  fifth  statements 
in  DO  loop  176  set 


FMl  «  ^  (7.122) 

W5  *  (7.123) 

where  fJJJ*  is  TTJI*  of  (3.110)  with  (i.m.n)  replaced  by  C;,p,9).  The  quantity 
7^^*  was  previouslv  calculated  and.  according  to  (7.37),  stored  in  ITME 
(JTME). 
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The  levath  throu|^  thirteenth  statements  in  DO  lot^  176  account  Iot 
the  first  and  third  terms  in  (7.121).  These  statements  are  not  executed  if 
either  p  »  0  or  f  as  0  because  these  terms  do  not  exist  then.  Quantities 

^  ^  because  there  are  neither  any 
TMof  nor  any  TM^  rectangular  waveguide  modes.  K  executed,  the  ninth 
through  12*^  statements  in  DO  loop  176  set 

mxM 

Z1  .  (T.124) 


r^tTM 

23  *  -BL_ 

(7.125) 

—  jSTM,TI4« 

(7.126) 

VVl  ““ 

(7.127) 

where  and  are  obtsined  by  replacing  (t,m,n)  by  C7,p,9) 

in  (3.57)  and  (3.63),  respectively.  In  the  13*^  statement  in  DO  loop  176, 

(7.128) 

W3.Z3  -  /f"™-  . 

(7.129) 

Thus,  the  13^  statement  in  DO  loop  176  adds  to  XTME  the  first  and  third 
terms  in  (7.121). 

Always  executed,  the  14*^  through  21**  statements  in  DO  loop  176  account 
for  the  second  and  fourth  terms  in  (7.121).  The  17*^  through  20*^  statonents 
in  DO  loop  176  set 

y^iTE 

(7.130) 

ym 

(7.131) 

W4  s 

(7.132) 

W2  « 

(7.133) 
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where  /JTE.'ni*  ^  ©bteined  by  repUdng  (i,m,n)  by  (J^p^q) 

in  (3.58)  and  (3.64),  respectivdy.  In  the  21**  statement  in  DO  Imp  176, 

-W2*Z2  =  (7.134) 

W4*24  *  .  (7.135) 

Thus,  the  21**  statement  in  DO  loop  176  adds  to  XTME  the  seomd  and 
fourth  terms  in  (7.121). 

There  are  possibly  four  pg*^  terms  in  (7.118).  They  are 

/f«™*  j ,  „d  /nt™.  p  13g) 

The  last  nine  statements  in  DO  loop  176  account  for  these  terms.  These 
statements  are  not  executed  when  r  s=  0  because  then  XTMO  is  not  needed 
because  there  is  no  Preliminary  to  the  calculation  of  the  terms  in 

(7.136),  the  eighth  from  the  last  statement  in  DO  loop  176  sets 

ws  =  Ij^  (7.137) 

where  /^JJ*  is  Z,™*  of  (3.111)  with  (t.m,n)  replaced  by  C7,p,g).  The  quantity 
77^  was  previously  calculated  and,  according  to  (7.39),  stored  in  FTMO 
(JTMO). 

The  sixth  from  the  last  through  fourth  from  the  last  statements  in  DO 
loop  176  account  for  the  first  and  third  terms  in  (7.136).  These  statements 
are  not  executed  if  p  s  0  or  g  «  0  because  these  terms  do  not  exist  then. 
Quantities  and  Vj™  do  not  exist  because  there  are 

neither  any  TMof  nor  any  TMpo  rectangular  waveguide  modes.  If  executed, 
the  sixth  ^m  the  last  and  fifth  from  the  last  statements  in  DO  loop  176  set 

\V3  s 
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(7.138) 

(7.139) 


where  ^  ©bUined  by  repUcmg  (t,m,n)  by  (j,P»9) 

in  (3.65)  end  (3.69),  respectivdy.  In  the  fourth  from  the  last  statement  in 
DO  lo<9  176, 

W1*Z1  »  7?'^™*  (7.140) 

W3*23  *  7*™-™* 

Thus,  the  fourth  from  the  last  statement  in  DO  loop  176  adds  to  XTMO  the 
first  and  third  terms  in  (7.136). 

The  last  three  statements  in  DO  loop  176  account  for  the  second  and 
fourth  terms  in  (7.136).  The  third  from  the  last  and  second  from  the  last 
statements  in  DO  loop  176  set 

W4  =  7/™-™*  (7.142) 

W2  a  7>TrE.TM.  ^7  J43J 

where  and  are  obtained  by  replacing  (t,m,n)  by  (jtPti) 

in  (3.66)  and  (3.70),  respectivdy.  In  the  last  statement  in  DO  loop  176, 

W2«Z2  a  7**^™*  ('^) 

W4.Z4  *  7f^™*  (^)  ■ 

Thus,  the  last  statement  in  DO  loop  176  adds  to  XTMO  the  second  and 
fourth  terms  in  (7.136).  It  is  evident  from  the  discussion  in  this  paragraph 
and  the  preceding  five  paragraphs  that  XTME  and  XTMO  are  given  by 
(7.117)  and  (7.118). 

The  first  nine  statements  after  statement  174  use  XTME  to  calculate 
C™*  of  (5.21).  TheM  statemenu  also  add  of  (5.44)  to  PIN,  |CJ^  + 

of  (5.34)  to  PRM,  and  of  (5.45)  to  PR.  The  quantities 

of  (5.44),  Prm  of  (5.34),  and  P,  of  (5.45)  will  be  accumulated  in  PIN,  PRM, 
and  PR,  respectively.  The  first  two  statements  after  statement  174  set 

XTME  «  C™*  +  B™*  (7.146) 

BT  *  B™*  (7.147) 
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where  Cf^**  is  given  by  (5.21).  The  third  statement  after  statement  174  sets 

PINS  *  (7.148) 

so  that  the  next  statement  adds  to  PIN.  The  fifth  through  seventh 

statements  after  statement  174  set 


CTME  =  (7.149) 

XTMES  *  jC™*  +  (7.150) 

CTMES  *  |C^|’  (7.151) 


so  that  the  next  two  statements  add  to  PRM  and 

to  PR. 


The  second  through  10*^  statements  after  statement  162  use  XTMO  to 
calculate  of  (5.22).  These  statements  also  add  of  (5.44)  to 

PIN.  |C^  +  of  (5.34)  to  PRM,  and  of  (5.45)  to  PR.  The 

above-mentioned  statements  are  not  executed  if  rsQ  because  there  is  neither 
nor  Bj^*.  The  second  and  third  statements  after  statement  162  set 


XTMO  =  C™'  -I-  Bji^  (7.152) 

BT  =  (7.153) 

where  is  given  by  (5.32).  The  fourth  statement  after  statement  162 
sets 


PINS  «  (7.154) 

so  that  the  next  statement  adds  to  PIN.  The  sixth  through  eighth 

statements  after  statement  163  set 


CTMO  =  C™*  (7.155) 

XTMOS  «  (7.156) 

CTMOS  =  (7.157) 

so  that  the  next  two  statements  add  4-  to  PRM  and 

to  PR, 
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7.6.2  The  Normalized  Coefficients  of  the  TE 

Propagating  Modes  in  the  Circular  Waveguide 

The  nonnalized  coefficients  and  <A  the  waves  of  the  even  and 
odd  TEra  propagating  modes  in  the  circular  waveguide  are  calculated  in 
nested  DO  loops  180  and  181  whose  indices  R  and  S  obtain  r  and  s  according 
to  (7.11)  and  (7.22).  These  DO  loops  also  add  to  W7  the  TE;.  and  TE^,  terms 
in  (7.116).  Specifically,  the  fifth  statement  in  DO  loop  181  adds  (|J9^P  + 
(I  to  W7.  Inside  nested  DO  loops  172  and  173,  nested  DO  loops 

Inside  nested  DO  loops  180  and  181,  nested  DO  loops  182  and  75  accumulate 
the  triple  summations  of  (5.31)  and  (5.32)  in  XTEE  and  XTEO,  respectively. 
After  exit  from  nested  DO  loops  182  and  75, 

s+f#o 

”“-SS  S 


The  calculation  of  the  right-hand  sides  of  (7.158)  and  (7.159)  is  described  in 
the  next  six  paragraphs. 

Inside  nested  DO  loops  182  and  75,  the  terms  in  (7.158)  and  (7.159) 
are  added  to  XTEE  and  XTEO,  respectively,  where  q  and  p  are  related  to 
the  DO  loop  indices  N  and  M  by  (7.119)  and  (7.120).  There  are  possibly 
four  terms  in  (7.158).  They  are 


(7.160) 


Preliminary  to  the  calculation  of  these  terms,  the  third  statement  in  DO 
loop  75  sets  FMl  equal  to  the  right-hand  side  of  (7.122).  The  fifth  and  sixth 
statements  in  DO  loop  75  set 


W5-/J 


W6  = 


rTEM 


100 


(7.161) 

(7.162) 


where  ai  (3.112)  with  (t,m,  n)  replaced  by  (;,p,  q)  and  is 

of  (3.113)  with  (t,m,n)  replaced  by  (i,p,9).  The  quantities  and 
were  previously  calculated  and,  according  to  (7.54)  and  (7.55),  stored 
in’  nE£P(JTEE)  and  ITE£Z(JTEE),  respectively. 

The  eighth  through  14*^  statements  in  DO  loop  75  acccnmt  few  the  first 
and  third  terms  in  (7.160).  These  statements  are  not  executed  if  p »  0  <» 
q  s  0  because  these  terms  do  not  exist  then.  If  executed,  the  10^  throu^ 
13*^  statements  in  DO  loop  75  set 


Z1  = 


ylTM 


ytnA 

W3  * 


W1  = 


(7.163) 

(7.164) 

(7.165) 

(7.166) 


where  jnM.VEt  ^  obtained  by  replacing  (i,m,n)  by  (i,p,q) 

in  (3.71)  and  (3.79),  respectivdy.  In  the  14^  statement  in  DO  loop  75, 


W3.Z3  = 


(7.167) 


(7.168) 


so  that  this  statement  adds  to  XTEE  the  first  and  third  terms  in  (7.160). 

Always  executed,  statement  175  and  the  seven  statements  after  it  ac¬ 
count  for  the  second  and  fourth  terms  in  (7.160).  The  third  through  sixth 
statements  after  statement  175  set 


ylTB 

(7.169) 

i/CTE 

Z4«-^ 

(7.170) 
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W4  s 
W2  » 


jlTB,TB« 


(7.171) 

(7.172) 


where  and  jtrt.'n*  ^  obtained  by  replacing  (t,m,n)  by  (j,Pt9)  “ 

(3.72)  and  (3.80),  respectivdy.  In  the  seventh  statement  aft«  statement  175, 


W2*Z2  *  (^) 

W4*Z4  * 

Thus,  the  seventh  statement  after  statement  175  adds  to  XTEE  the  second 
and  fourth  terms  in  (7.160). 

There  ate  possibly  four  pq*^  terms  in  (7.159).  They  are 


jtTM.TE» 


jtTH.VU 


(7.175) 


The  last  10  statements  in  DO  loop  75  account  for  these  terms.  These  state¬ 
ments  are  not  executed  when  rsO  because  then  XTEO  is  not  needed  because 
there  is  no  C^.  Preliminary  to  the  calculation  of  the  terms  in  (7.175),  the 
ninth  and  eighth  from  the  last  statements  in  DO  loop  75  set 

W5  -  (7.176) 

W6  «  /™r  (7.177) 


where  is  of  (3.114)  with  (t,m,n)  replaced  by  C7,p,7)  and  is 
of  (3.115)  with  (i,m,n)  replaced  by  (i.p,?).  The  quantities  and 
were  previously  calculated  and,  according  to  (7.57)  and  (7.58),  stored 
in  ITEOP(JTEO)  and  ITEOZ(JTEO),  respectively. 

The  sixth  from  the  last  through  fourth  from  the  last  statements  in  DO 
loop  75  account  for  the  first  and  third  terms  in  (7.175).  These  statements 
are  not  executed  if  psO  or  f  «0  because  these  terms  do  not  exist  then.  If 
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executed,  the  sixth  from  the  last  and  fifth  from  the  last  statements  in  DO 
loop  75  set 


W3  == 


r2TM,TE« 

*3 


(7.178) 


W1 


_  rl7M,TEo 
—  *3 


(7.179) 


where  and  /JTM.tEo  ^  obtained  by  replacing  (t,m,n)  by 

in  (3.81)  and  (3.86),  respectivdy.  In  the  fourth  from  the  last  statement  in 
DO  loop  75, 


W1*Z1  =  1 

3  1 

1  / 

-W3*Z3  =  1 

1  / 

(7.180) 

(7.181) 


so  that  this  statement  adds  to  XTEO  the  first  and  third  terms  in  (7.175). 

The  last  three  statements  in  DO  loop  75  account  for  the  second  and  fourth 
terms  in  (7.175).  The  third  from  the  last  and  second  from  the  last  statements 
in  DO  loop  75  set 


W4  = 
W2  = 


(7.182) 

(7.183) 


where  and  are  obtained  by  replacing  (»,m,n)  by  {j,p,q)  in 

(3.82)  and  (3.87),  respectivdy.  In  the  last  statement  in  DO  loop  75, 


W2«Z2  «  1 

/y^lTEx 

3  1 

1  y/V  J 

-W4.Z4  =  /”E.TE»  1 

3  1 

1  yA  J 

(7.184) 

(7.185) 


so  that  this  statement  adds  to  XTEO  the  second  and  fourth  terms  in  (7.175). 
It  is  evident  from  the  discussion  in  this  paragraph  and  the  preceding  five 
paragraphs  that  XTEE  and  XTEO  are  pven  by  (7.158)  and  (7.159). 

The  first  nine  statements  after  statement  182  use  XTEE  to  calculate 
of  (5.31).  These  statements  also  add  of  (5.44)  to  PIN, 
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of  (5.34)  to  PRM,  and  of  (5.45)  to  PR.  The  quantities  of  (5.44), 

Prm  of  (5.34),  and  P,  of  (5.45)  are  being  accumulated  in  PIN,  PRM,  and  PR, 
lespectivdy.  The  first  ti*o  statements  afta  statement  182  set 


XTEE  *  (7.186) 

BT  s  (7.187) 

where  is  given  by  (5.31).  The  third  statement  after  statement  182  sets 

PINS  *  (7.188) 

so  that  the  next  statement  adds  \Bj^^  to  PIN.  The  fifth  through  seventh 
statements  after  statement  182  set 

CTEE  =  (7.189) 

XTEES  -  |C™*  +  (7.190) 

CTEES  -  (7.191) 


so  that  the  next  two  statements  add  \C^  +  ^  ^ 

PR. 

The  second  through  10*^  sutcments  after  statement  165  use  XTEO  to 
calculate  of  (5.33).  These  sutcments  also  add  of  (5.44)  to  PIN, 

|C^  +  of  (5.34)  u  PRM.  and  |CCS®*|’  of  (5.45)  to  PR.  The  abo^ 
mentioned  statemenu  are  not  executed  if  rvO  because  there  is  neither 
nor  Bj^.  The  second  and  third  statements  after  statement  165  set 


XTEO  mCjy^Bj^  (7.192) 

BT  -  S”*  (7.193) 

where  is  given  by  (5.33)  The  fourth  statement  after  statement  165  sets 

PIN'S  -  \Bjyf  (7.194) 

so  that  the  next  statement  adds  to  PIN.  The  sixth  through  eighth 

statements  after  statement  165  set 

CTEO  -  C?*  (7.195) 

XTEOS  -  |Cj**  -►  (7.196) 

CTEOS  -  |CJV|*  (7.197) 
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so  that  the  next  two  statements  add  |C^  +  to  PRM  and  to 

PR. 

After  exit  from  nested  DO  loops  180  and  181, 


PIN  «  /W  (7.108) 

PRM  *  Pr^  (7.199) 

PR  *  Pr  (7.200) 


where  Prmt  *n<l  Pr  ue  pven  by  (5.44),  (5.34),  and  (5.45),  respectivdy. 


7.7  All  Statements  Beyond  Statement  190 
in  the  Main  Program 

The  group  of  statements  beginning  with  the  statement  after  statement  190 
and  ending  with  the  last  statement  in  the  main  program  are,  except  for 
minor  differences  noted  in  the  following  paragraphs,  the  same  as  in  the  main 
program  of  [3]. 

The  second  statement  after  statement  190  is  new.  This  statement  sets 

PINC(KA) «  .  (7.201) 

No  value  of  BKAG  is  calculated  after  statement  91  because  BKAG  no 
longer  appears  in  the  two  statements  that  define  PTA  and  PRMA.^  These 
two  statements,  immediately  following  statement  91,  set 

PTA  «  Pu  (7.202) 

PRMA  «  Prr^  .  (7.203) 

where  Pu  and  are,  in  omtrast  to  eqs.  (3.316)  and  (3.323)  of  [3],  given 
by 

-  -ta [^[j  (7-204) 

iU.  -  -R.  -  Im  (-jV5«  (^)  }  ■  (7-205) 

^Tbe  left-baod  side  of  eq.  (3-317)  of  (3]  is  iocorrect.  This  left-hand  side  should  be 
BKAG  instead  of  BKAB. 
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In  (7.204),  yj*  i«  the  same  as  [K*  +  y*]«  in  eq.  (3.316)  <rf  [3]: 

ir=(i''+j^L-  (t'Sot) 

To  gain  physical  interpratations  oi  the  quantities  Pu  and  Prm»i  ^  first  recast 
(7.204)  as 


Pu-^pS^+pS^ 


(7.207) 


where  and  P^^  are  given  by  (4.12)  and  (4.17),  respectively.  As  stated 
in  Section  4.3,  is  an  alternative  expression  the  tinie>average  power 
radiated  by  the  magnetic  current  in  the  left-hand  rectangular  wavq^de 
and  pS^  is  an  alternative  expression  for  the  time-average  power  radiated  by 
the  magnetic  current  in  the  ri^t-hand  rectangular  waveguide.  Hence, 
Pta  is  an  alternative  expression  the  time-average  power  radiated  by  the 
magnetic  currents  in  the  rectangular  waveguides.  Using  (7.207),  we  see  that 
the  right-hand  side  of  (7.205)  is  the  same  as  that  of  (5.43).  Hence,  P„^  of 
(7.05)  is  an  alternative  expression  for  the  time-average  power  radiated  by  the 
magnetic  currents  in  the  drculax  waveguide. 

Upon  exit  from  DO  loop  101,  we  want  to  have 


E3A1P(J)  * 

E3A2P(J)  = 

E3A1Z(J)  = 

E3A2Z(J)  = 


£<^»)(^M»),0: 


|£(3+)(^2) 

1 

£!'»>(0,xS''>) 

l£(*+)(^,Q)| 


(7.208) 

(7.209) 

(7.210) 

(7.211) 


where  the  right-hand  sides  of  (7.208)-(7.211)  are  given  by  (6.6)-(6.9),  re¬ 
spectively.  According  to  eqs.  (3.369)-(3.372)  of  [3],  the  main  program  of  [3] 
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Hurfji 

(7.212) 

(7.213) 

EJAizw)  - 

(7.214) 

E3A2Z(J)  . 

(7.215) 

where  the  right-haiid  tide*  of  (7.212H7*215)  ere  given  by  eqs.  (3.328), 
(3.332),  (3.336),  and  (3.339)  of  (3),  respectively. 

The  main  program  treat*  each  parcnthesiaed  quantity  proportional  to  one 
of  the  V'»  in  (6.6)**(6.9)  exactly  the  tame  a*  the  main  program  of  [3]  treat*  the 
corresponding  parenthenaed  quantity  in  eq*.  (3.328),  (3.332),  (3.336),  and 
(3.339)  of  [3].  For  example,  the  main  program  store*  (6.6)  and 

(6,8)  in  the  same  variable  a*  the  main  program  of  [3]  stores 
of  eqs.  (3.328)  and  (3.336)  of  [3|.  Hence,  a*  far  as  both  main  program* 
are  concerned,  the  right-hand  sides  of  (6.6)-(6.9)  are  those  of  eq*.  (3.328), 
(3.332),  (3.336).  and  (3.339)  of  (3j  with  5.  of  eq.  (3.329)  of  [3]  replaced  by  5. 
of  (6.10).  The  difference  in  5«  is  accounted  for  in  statement  155.  Statement 
155  in  the  main  program  of  |3|  and  the  statement  after  it  have  been  replaced 

^The  ri^(*baad  wOw  of  «ui  (3  9ttM3173)«f{3]ai«uieocr*ct.  The  collect  right-band 
•ides  are  the  ocgativas  of  Uwar  vniua  la  thew  >qaatioos.  Eqnatioos  (7,212)-(7,215)  aie 
the  collected  eqs.  (3  )69).  (3  371)  (3  370).  and  (3373)  of  [3],  lespeetivdy.  The  sign 
enois  in  eqs.  (3.369H3373)  of  (3)  veie  dee  le  incorrect  signs  on  tte  left-hand  sides  of 
e<H.  (7.12>-(7.14)  of  p;  la  fact,  every  r/* '(•'/' ’.0).  £l"’(d5^’\0),  £{-***(»,  and 

£<aa)(o 

in  Scctwa  3  10  of  (3)  has  the  wioog  sign  attached  to  it.  Hoeever,  nested 
DO  loops  101  and  104  of  the  maia  peegram  of  p]  still  accumnlate  the  ri^ii-baad  sides  of 
eqs.  (3328),  (3.332),  (3  336).  and  (3  3)9)  of  (3)  in  E3A1P(J).  E3A3P(3),  E3A1Z(3).  and 
E3A2Z(J),  rcspectivd|r 
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by  fUtement  155.  The  Utter  etetrment  155  aets 

SA  «  5.  (7.216) 

where  5«  is  pven  by  (6.10)  instead  of  eq.  (3.329)  of  {3}.  In  this  statement, 
W7  is  given  by  (7.116). 

In  the  main  program  of  [3] ,  the  value  of  the  time-average  power  incident  in 
the  circular  waveguide  is  not  written  out  because  the  time-average  transmit¬ 
ted  and  reflected  powers  calculated  therein  are  thoae  per  unit  time-average 
incident  power.  In  the  main  program,  however,  the  value  of  the  time-average 
incident  power  is  written  out  because  the  time-average  transmitted  and  re¬ 
flected  powers  calculated  therein  are  absolute  time-average  powers  rather 
than  ratios  of  them  to  the  time-average  incident  power.  The  time-average 
incident  power  in  the  main  program  depends  on  the  anq>litudes  ci  the  waves 
of  the  propagating  inodes  incident  in  tJte  circular  waveguide.  The  statement 
after  statement  149  writes  out  (PINC(I),  I »  1,  KAM)  where 

PINC(I)  *  (7.217) 

where  is  the  time-average  incident  power  given  by  (5.44)  when  the  index 
KA  of  DO  loop  48  is  I.  When  KA  x  I,  Ire  is  given  by  the  right-hand  side  o( 
(2.15)  with  KA  replaced  by  I. 

7.8  Listing  of  the  Main  Program 

niPuczT  iuL«e  (i-i.o-z) 

coRNOi  /iiaocs/K.nia.KTii.m.iDi(so).Biii(ioo).Biin(ioo) 

OMIIOI  /fl/fl 

commi  /nuz/iiuz 
coiin)i  /us/in.sHU 

comoi  /pix/u.iis.nui.t.scft.pii(iso) .raadso) .rasdso) . 

IPMdSO) 

comoi  /Dei/S.'^n.U.C.CI.nB.M(SO).^4(M).P6C 
COMPLIX^ie  ZL%  .a3.0.UD.S4.TTS.TTI.K1l(S0)  .0Tm.0TI(M}  .OTB 
CORPLEKie  03.03I.Zl.Z2.a.Ze.S.Sl.S3.S4.SS 
C0IIPLS*l  miP(100).t3AlZ' 100). 13321(100) .I3A3Z(100) 

C0mLn*l  Tm.TBI.Tm,m,Ti34330).n(lM).T(lM) 
comLn«8  cna.cTK.cno.ciooT.cizi.czooT 
coiinjz*8  c2zi,yuc(iM) 

IE3L*S  Ll.L2.L3.ZJ(300).ZJP(300).eT1l(B0).6n(S0) 
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UtLf  1IP(M).1HI(S0).TIP(S0).TBi(S0).0qn(S0).0QTI(B0) 

iiu.*4  raxi(ioo).raxa(ioo).z(ioo) 

iEiL«4  mAi(ioo).nirL(ioo).Bun.T(ioo) 

iiAL«4  ciaoTs.eizis.canns.cazis.FT.eTaB.CTBs.cTios.Hi 

iiu.*4  szip(ioo).sziQ(ioo).isAin(ioo) .nzandoo) .mizscioo) 

BIAL«4  B3A2ZS(100).8Ziqq.SZirP.UIJP 

UAL«4  inuq.BnjK.iajQQ.Bajrq.aaijqp 

zmooi  i.fti.s.auz.p.nuz.pi.P2.P3.Q.qi.pni.FTi.Qni.QR.xrs(iM) 

zmeB  qra.nsdoi) 

zraen  tiM.m.mzq.sndO.sndo) 

zmen  txin.iTni.smdO.smdo) 

IBALM  Zm(3000).ZTRO(3000).ZTBI?(3000).ZTIB(aOOO) 

iiAL«4  ZTBOP(20oo).znca(aooo).?zi.m.raMi.Fxicdoo) 
iiAL*4  mdo.io).mdo,io) 

coiiPLa«8  Bmdo.io).Bnodo.io).inEdodo).iRodo.io) 

COIIPLn«C  BT.ZTMB.ZTKO.mE.mO.CnO 
C(MPLn«8  Z1TM.UTE.Z2T1I.Z3T8 
GPBI(aiIT*20  .FILE-  *  Z1  .OAT*  .STATDS*  *010' ) 

0PEI(0IIT-21  .FZI£>*(KJT.0AT*  .STAT1»«*0L0*) 

1810(30.10)  B.C.U.L2.U.BEM.ZN.ZL1.ZU 

10  FQMUT(4014.7) 

V1ITE(21.11) 

11  FQllUT(*B.C.U.U.U.BEM.ZM.ZLl.aL2') 

VtZTECOl.lO)  8.C.U.U.U.BSM.ZN.ZL1.ZL2 
1810(20.144)  EAN.BE10.0B81.I£3N.IPn.IZ 

144  r011UT(Z4. 2014. 7.3X4) 

«irrE(21.14S)  KlN.KlO.OiXl.KESN.lPn.R 
146  F01M1T(*E1N«'.X4,*.  U10>'.014.7. * .  0B81>' .014.7/ 

1'8E3N>*.I4.'.  fP8I>*.Z4.*.  1Z*'.Z4) 

1810(30 . 14«) (883(1) .I«l .SE3M) 

146  F01IUT(1SI4) 

VirrE(21 . 147) (883(1) . I-l .8E3N) 

147  F01IUT('883*/(1SI4)) 

PI«3 . 141S93663688700«0 
BC-l/C 

PC>PI*BC 
B81I2>B8N«B8N 
CALL  MOOES 
VirTE(21.103)  BMIZ 
103  F01M1T('BM1Z>*.I4) 

V1ITE(21 . lOi) (MM(I) .Z«l .BMIZ) 

108  F01N1T(*MM*/(1SI4)) 

B1*BTII*BT8 

V1ZTE(21.183)  8TM.8TE.81 
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ISS  muUTC Z4.*.  RI-.Z4,*.  I1-M4) 
VdHAZ.OT.O)  00  TO  Its 

mzn(2i.tio) 

110  rauuTC'in  n  too  auuo 

S1W 

lis  CUL  WtMU 
flTmfll.lhC 

pzBc«n*ic 

PZS>PZ*.SD«0 

Iil«l 

00  40  KA-l.IAM 

KA«Ki04(IA-l)«0IXi 

MU2-KA*MCA 

ZP(KA.eT.1.04110370D40)  00  TO  00 
OtZTKSl.OS) 

OS  POaXiT(>IXA  ZS  TOO  SHALL*) 

STOP 

04  ZF(C.tT.B)  00  TO  00 

Hkznai.oo) 

00  POUUT(*C  ZS  tor  USS  TIAI  1*) 

STOP 

00  ID«BU*1 

Htzncsi.iso)  to 

ISO  P(nRAT(*BO«*.tl4.7) 

ZP(BO.GT.PZ)  00  TO  110 
vuTKai.iii) 

111  POI1UT(*BO  ZS  TOO  SHALL*) 

STOP 

110  ZPCtO.LT.PZa.AB.tO.LT.PZiC)  00  TO  113 
HtZTSCai.llS) 

113  FOSRATC'tO  ZS  TOO  LAlOt*) 

STOP 

113  B03>KB«ta 
10*1. 0^/10 
0>(0.IH0.1.»*0) 

BO*-KE*0 

BS«t«.S»*0 

BZS«0ASZt(BS) 

BZ«3.IK0«tSS 

ZB-l.lHO/U 

Zl«Ll/i-ZB 

Z3«L3/t*U 

JTE-0 

JTH«0 


no 


00  13  Q-I.nux 

F3-1 

xr(Q.«.i)  OM 

n-«(Q) 

00  14  p^.oa 

JTI*JII^l 

JTIl«JTS«Rli 

JIB-JTU^Il 

aM3>Mn(JTi)-BBa 

zr(o.a.a.oft.q.B.i)  «o  n  is 

oiT»oaq>T(-«i«a) 

Ai«osr*xi 

eA^)COS(Al} 

SA«0tZl(Al)*0 

si-oiT»an 

niC(JTIl)*(CA<»ZLl*SA)/(SA4ZLl*CA)««l 

Aa>air*x2 

CA«0e0B(A3) 

SA«osa(Aa)*o 

TtlCCJTI3)«(eA4ZL2«SA)/(SA4n.2*CA)«Sl 

00  TO  17 

IS  OAM-oaotTCOAiia) 

Trt»-OAX*SKI 

nic(jrni)«rTi 

7IIC(JTI2)«m 

17  ZF(P.Ki.l.OS.Q.B).l)  00  TO  14 
TTR-OD/OAR 

Tni«JTll<»l 

JT1I2«JT1I>K1 
TUC(JTM)*TT1I 
TUC(JTII2)-mi 
14  COITXlOB 
13  COITZIDB 
I3>K1«3 

VXZTI(31 . 100) (TIIC( J) . J>1 .13} 

100  F01fUT('TUCV(«14.7)) 

CS«C«.SIHO 

2SS>1.0^/CS 

niAZ«iii(i} 

XZ-S/OX 

Tzni-s. •sxs«osQEr(sc/ox} 

Tzn-Tzm/sc 

TA«FX2/0XA 

302*1 .0«0/0SQET(2.0^) 
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Tci«osqKr(n«K) 

TONTCl/SQa 

1C1«TC1*IKA 

ait^ 

CStafMQKrd 

KMQ4C1 
M  13  ZT»l.n 
Y<XT)«0. 

13  COWUW 
tHNl 

•nd)^ 

in«i 

md)-o 

JTMI-O 


nw>o 

nUX3«fMAX*3 

T«i»qKr(t.iKO*rxie) 

S03*-1.0«0 
00  10  Bffl.SOO 
3Ga«-S63 

ki>t-i 

tS«kl«tl 

CIXI.  BIS(k.ZJ.XJP) 

Zr(SlUX.QT.300)  STW  M 

xr(SMZ.iq.o)  eo  to  3S 
cui.  rii 

lM«kl«BU 

csiaocosdtc) 

SBl»DSZBdIC) 

DO  114  N«1.PIUZ 
RJ«B«nUZ3 

nii(iu}-rai(R)«csi-m(i)«sii 
ra3(iu)>ra3(H)«csi«rai(R)«sii 
PI3(IU)«PIS(N)«CSl-ra4(R)«Sll 
n4(RJ)>ra4(R)«CSl^raS(R)«SRl 
114  COBTZm 

00  30  S«1.SRAZ 

CALL  06Id.ZJ,XXTR.mi.eARTR.Tir.TRR,DTR,6TR.0qni.aCSTR.aC31R. 
1ZBITR.ZZTR .  ZOITR.ZOOTR) 

CALL  0Cl(3.XJR.Zm.ZTI.eARTE,TV.TER,inS,eTB,0qTl.aCSn.0C3TB. 

iziin.zzTB.zom.xooTi) 

u-xm-is 
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«a»>ia/xaNum 

W6«n«xni 

VMM/U 

«S«Et««K 

M*zsa/xxR*ta 

nmmua/QMKn 

HS-VS/OMTI 

xrcii.n.o)  00  TO  46 

«1«V1*.CIH0 

VMIS«.60«0 

m^w.ny^o 

40  ir(nM.IQ.2)  60  TO  os 
OINOAMXV/IU 
T1S-T/DSQ>T(6II) 
XF(ftl.lQ.O}  TMS«TRS*SQ2 
ITM(8.ft)*0> 

STM-t 

sntwms 

KTE*0 

00  39  ■•i.nux 

II2>1 

Tni>TIB*6TR(I) 
ZF(II.IE.O)  60  TO  76 
II3«3 

Tiii«Tia*sq3 

76  Ms-mid) 

00  63  RalQ.IO 

MTE-m^l 

Rl-N-1 

TTR>TRI/BRI(RTE) 
Zr(Rl.BQ.O)  TTR«TTR«SQ3 
JTRE>JTRE«1 
RJaR^mX3 

ZTRE(JTRB)>TTR«n3(RJ) 
ZP(Il.IQ.O)  60  TO  S3 
JTR0«JTR0«1 

ZTR0(  JTR0)-TTR«m  (RJ) 
63  COITZIDB 
39  COITXIDB 

66  ZF(ZTE.IQ.3)  60  TO  68 
6B>8U/6ARTI 
T8S«T/DSQRT(CE*ZR) 


Xr(Il.|Q.O)  m«T15*302 
■R(s.i)«ai 

TTIZ»1IS*lZ«ZXn/0AITI 

in«ft 

sn(»)>s 

MTIN) 

00  TT  ■•l.nu 

Ma«i 

nVI«llV*6R(I) 

TTm«TTIZ«04(I) 

zrcfi.n.o)  oo  to  er 

■3^ 

TTIFIaTTlfl«SQ2 
TTixi«rnm*sQ3 
67  Ks-nid) 

00  70  ■•■3.MS 

onni«Mi(Ri) 

Townw/omni 

Tiz»iTnzi/OMni 

ZF(Ml.ll.O)  60  TO  71 

Tira«Tini«sQ3 

Tizii«txzR«sqa 

71  JT1I«JTII«1 
RJ«H4nUZ2 

XTII0(  JTII)>Tini*rai  (El) 
ZTiiz(jni)«Tizii*ra6(u) 
Xr(Ot.lQ.O)  60  TO  76 
jnoajTVKi 

XTI00(JTI0)*Tini*?l3(U) 

ZTIOZ(JTIO)>TIZII«m(IJ) 

76  COITZIDI 

77  COmiOE 
66  IT1I*0 

00  31  ■•i.nuz 

>3>1 

ll«l*l 

XF(II.IQ.O)  ■3«3 
ii3>ni(i) 

070*0 
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0 


9 


Qn*o 

ntiB>ii*ic 

IIO«l-Hll-3«(ll/2) 

00  ID  (18.S1).  zn 
It  Taai«iiai(i) 
nra>iv(i) 
imoMira(i) 

80  ID  03 

81  oom-oQTad) 

TKi»8cnM»Dqnai 

na-oesn^ooni 

83  80  TO  (84.88).  ZTl 

84  Tiai>in(i) 
nra>iip(i) 

OTa«on(8) 

D38*03(l) 

80  TO  88 

88  DOTM«oqTl(l) 

Tll«8CSTI*0QTa 

TI3«8C3TB«0QT1I 

ni*8i 
ra6«rii«88c 
88  M  22  Q«l.nUZ 
P2-1 
Ql-0-1 

Zr(Ql.lQ.O)  r3«3 
p3«iai(Q) 

V8«Ql*ni8 

IQBO«lEO^qi-3«(01/3) 

IIQQ«2 

ZFCll.Sq.Ql.ilD.Ql.IE.O}  IEOQ-1 
QP1«2 

ZrCQl.IQ.O.UD.ll.EQ.O)  QPB«1 
60  TO  (31.32).  ZTR 

31  ZlI-11-01 
Zr«I14Ql 
Tllliq>110l(Q) 

TIIPQ«T1ir(Q) 

FTM-rardP  .TW8  .TIIIIQ)-m(X«.TIIIII  .-THHQ) 
i-m(zp  .Tma  .TiirQ)«nT(  ZR.Txra  .-THPO) 
Zl-tfl*(FTN<»6TII(Q)*DTlll) 

60  TO  47 

32  60  TO  (48.107.117).  8QC0 
48  60  TO  (118.119).  QPl 
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4 


118  ZIIpZZTX 
60  TO  120 

119  ZIBpZBTM 
00  TO  120 

107  ziMzasni 

GO  TO  120 

117  Zll*ZaOTK 

120  zii*zi8*T«i*oqni(q) 

GO  TO  (121.122).  RQQ 

121  Zll>Zlt-TK2 

122  ZlB»lll«Zlt 

47  60  TO  (36.37)  ITZ 

36  6TV)«6TI(Q) 

040^(0) 

ZM>I1-01 

ZP>I14qi 

TnQ«Tza(Q) 

TIPQ«TO(Q) 

F1«PZT(Z1.TIIII  .'TIMQ) 
r2*FZT(zr.Tm.TnQ) 
F6«rZT(XF.TIia.TDQ) 
r7«rzT(zi.Tira.*TOQ) 

Pi-n-f'i 

F0*F2*ri 

Fl-r6Hr7 

F2-F64F7 

m-F8-Fl 

F3*F9«F2 

F4«F8'^1 

F8-F9-F2 

Z2«ir2«(m*6TlBiOIU) 

Z3«V3«(F3*GT»«On) 

Z4— W(F4i6<0«OTU) 

ZS*«f  (F8»G60»0») 

XFdl.ZQ.Ql.m  81.88.6)  ZM«6«ZI 
60  TO  49 

37  60  TO  (139.124.138).  8QB 

123  60  TO  (126.137).  OFO 

136  ZI«ZZT1 
60  TO  138 

137  Zt«ZZZTf 
60  TO  138 

134  Z8>Z0m 
60  TO  138 
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126  ZMZOOR 

126  zi*a«Tii*DQn(q) 

Z2I«ZE 

00  TO  (129.130),  nOQ 

129  Z2B«Z2ft<-TI2 

130  231»PIQ*Z31 
Fqi«qi 
pqo>rqi*Kc 
zap9qo*z23 
ZSft-Fie*PQG*ZK 

00  TO  (131.132).  IBOQ 

131  Z8t>Z8a«TB2 

132  Z2I*>H>Z21 
Z3a>ir3*Z3ft 
Z41«W3*Z4i 
Z63»U6*ZS& 

00  TO  (133.49)  IBOO 

133  ZS3«V9«ZSI 

49  ZP((nM4>ZTB).re.3)  00  TO  97 
00  TO  (63.64). no 

63  Z3«Z23 
Z3«Z3t 
Z4«Z4& 

Z6«Z6E 
00  TO  07 

64  Zl-ZIE 
97  MTH-fTK 

MTB-nS 
DO  23  II>II2.M3 
EMI>i 
M1>N-1 

ZF(ll.BQ.O.OE.Ol.iq.O)  00  TO  20 
EMI<i2 
OTM-inO^l 
26  OTE-mO^l 
V9«M1«Q1 
TK1E>K1/9C 
TB>TE/BIII(ini) 

PTflaQTO 
PTB-QTB 
DO  24  F>F2.73 
KPQ«1 
P1»P-1 

ZF(Ql.IQ.O.OE.Pl.EQ.O)  00  TO  27 
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IM)«2 
rai«raKi 
ar  nB«ni«i 

n>n/in(rTB) 

XF(HI.IQ.O)  T1«T1*8Q3 

zp(ii.v.o)  n>n*SQ2 
xr(Pi.»).o)  Ti«Ti*sq2 

XP(Ql.lQ.O)  T1*T1*SQ2 

raip«pii(p) 

nDP«pf2(p} 

piap*pis(p) 

PI4P«PI4(P) 

o*i«aaii^Kpq)/4 

IX-(PTII-l)*Ka 

BBlaim^RlB 

aR«ni«KiM 

nM«KiHni 

Kia*(KTM^PTB-l)HC3 

KIIB«BB2<»imi 

in»ES4>ni 

K«0 

HJaN 

DO  38  J«l.a 
Plliu-PIKMJ) 

PI3IU>PI3(RJ) 

PB3IU«PI3(IU) 

PI4IU«P14(NJ) 

PAC1«PI2P«PI3MJ-PI1P*PB41U 

PA63«PB2P«PB2MJ«PBSP«PB1MJ 

Pi63«PB4P*PBlMJ^B3P*PB2HJ 

PAa4«PB4P«PB4MJ^PB3P«PB3IU 

Xr(J.B).l)  60  TO  30 

PA6la^i61 

PA04>^A64 

30  2F((nE4ZTR).Eq.4)  60  TO  68 
50  S1*P162« (21-23) 

S3>PA61*23 
S4«Pi63*24 
S6«Pi64*26 
64  60  TO  (33.38),  KB 
38  Xiai«KlBl4K 

Tiai«Tl*  (  V8*Sl4ira*S3-V10*S4-Vll*S5) 
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r(2W)*T(2n)4TM 

33  00  TO  (34.30)  .DQ 

30  zn-ni^ 

*  TI|iin*(1l0^1-V0*S3-«ll*S4^t0*S5) 

T(Zn)«T(ZIR}«TBi 

34  00  TO  (30.40) .HOI 
40 

Tn«Tl*(«10«Sl^lfll«S3«ll0«S4«V0*S5) 
T(ZB)>T(Xn)«TK 
30  m-KB^K 

TII>T1*(V11*S1-V10«S3^II0«S4>W*S0) 

T(III)-T(XII)«TB 

00  TO  07 

00  S10«000a« (310-230) 

S30-P001«330 
S40«P003«340 
SU«»Ae4«Z00 
00  TO  (00.00).  n 
00  xnianoi^i 

Tm-Tl«(«0«S10««0*S30-«10«S40-«ll*SOO) 
T(Z)n)«T(zni)«THi 
00  00  TO  (00,01).  tfQ 
01  XIM>KIM«I 

Tn«Tl«(VO«OlO-«l«S30-«ll««40««10*tOO) 
T(Zn)«T(X3R)«tai 
00  00  TO  (03,03).  DU 
03  XMB-DB^I 

TllE«Tl«(W10*Slt««ll*S30*WO»S41««0*SOt) 

T(XllB)*T(Xia)«TK 
03  XEE>Ka*l 

TIE>T1* («1 1 •Slt-«10*S3t*«0*S40-«0*S00) 
T(Xn)«T(XII)«TB 
07  Ml 
iu>N*nux 

30  corrxiol 

'  34  COITXIOE 

33  COITXlOB 

Qm-rni 

QTB-ni 
33  COlTXm 
ITM-RIII 
ITB-IRI 

31  COmiOK 
30  COITXfDI 
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IB  canuiii 

II1XTI(31.44) 

44  FQMUT(*&  IS  TOO  U18B>) 
STW 

2S  Xr(S.B.l)  00  TO  41 
1iBXTS(ai.42) 

42  rOUUT  (*ZS  IS  TOO  SMALL*) 

sm 

41  14<l«sa 
il«0 

00  so  >1.11 

00  SS  X«1.I1 

J1«J141 

J2«J14S1 

JS^Jl-HCS 

J4«JSHC1 

T(JS)«T(Ja) 

T(J4)»T(J1) 

88  COMTIMOS 
J1»314M1 
SB  COITXBUB 
SB-S2«1 
XT»1 

00  43  X«1.B2 

T(XT)*T(XY)4TSSC(X) 

XT-XTHB 

43  COirXIOB 

00  IBl  E«1.STM 
SIUX«STN(1} 

Xr(SMU.BQ.O)  CO  TO  IBl 
00  1S2  S-l.SNAZ 
BTIIE(S.S)«0. 

BTNOCS.D-O. 

102  COlTXlOB 

101  cosmos 

00  103  S*1.STI 
SMAX-srS(S) 

Xr(SIUZ.BQ.O)  GO  TO  103 
00  104  S«1,SIUX 
BTBE(S.S)«0. 

BTE0(S.S)*0. 

104  cosmos 

103  COSTXSOS 
lSA0(a0,100)  STNS.STES 
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iM  raauT(ioz4) 

HUTKai.M)  nn.im 

66  faUUTClXn»M4.>,  tTli*M4) 

tiABCao .  106)  (sm(i)  .I*!  .tin) 
lltZTI(21.rf)(ST«(E)  .1^1. moi) 

79  rOUUT(*SnE*/(10Z4)) 

•UB(».  109)  (smcft) .  ft*i  .im) 

91XTI(21  .M)  (sm(t)  .1*1  .HIE) 

M  7aEEiT(*STlE*/(10Z4)) 

00  166  E«1.ETME 

sfliz«snat(E) 

Zr(SIUZ.IQ.O)  M  TO  166 

E1«E-1 

ElAO(aO.  164)  (BTH(S.E)  .S-l.SEAZ) 

164  F0EMAT(4014.7) 

1IEZTS(31.90)  E.(inB(S.E).  S«l.SEi2) 

80  F0EIUT(>BTEBC.Z4.*.S)V(4D14.7)) 
ZF(E1.1Q.0)  GO  TO  166 

EBA0(a0. 164) (BTMO(S.E) .S«1,SIUZ) 
VEZTB(31.16S)  E.(ETE0(S.E),  S«1.SEAZ) 
163  FOEIUT('BTMO(«.Z4.*.6)V(4D14.7)) 

166  COlTZlOB 

00  167  Eal.EIEE 
SMU«STBE(E) 

ZFCSEU.BQ.O)  60  TO  167 
E1«E>1 

EEA0(30 . 164)  (BTEECS  .6)  .S«1.S1UZ) 
VEZTE(21.184)  E.(BTIS(S.E).  S«1.SNU) 
184  F0EIUT('BTSE('.Z4,'.S)*/(4014.7)) 
ZF(El.Eq.O)  60  TO  167 
IEA0(30 . 164) (BTG0(S .E) .S-l.SlUZ) 
VEZTE(21.61)  E.CrmXS.E).  S't.SlUZ) 

81  F0EIUT('BTB)C.Z4.*.S)'/(4D14.7)) 

167  comiuE 
«EZTE(31.186)  ETK.ETE 

186  F0B1UT('ETE>'.Z4.*.  ETE-'Z4) 

VEZTE(21 . 166) (STN(E) .I*! .EXE) 

186  F0EIUT('STM*/(10Z4)) 

VEZTBCEl , 187) (STE(E) .E«l .ITE) 

187  F0EIUT('STl*/<i,014)) 

IITM*0 

IRS*0 

00  16  i-i.nuz 

II3«1 


121 


ni«ii 

XP(ll.Bi.O)  12*2 

■MHCI) 

00  74  MB.n 

rai-m/M 

zin-o. 

ZlT|i«. 

Z2ni«o. 

Z2T|i«. 

SOMl. 

MTI«RI«1 

Jni*MTI 
jno«Hn 
00  72 

set««9aft 

s>Az*sni(ft) 

Xr(8IUZ.IQ.O)  60  TO  72 
00  73  SM.SlUZ 

zi«nn(miB)*tiiii(s.t) 

vi*s6i*ni 

V2>S6l*r« 

ZlTM«Hft*Zl«ZlTI 

ZlTI«-«2«Zl«ZtTI 

Z2TN«ntl«Zl*Z3TW 

Z2TI-f1ll*Zl*Z3TI 

Zr(El.lQ.O)  60  TO  73 

zi«zTio(jTin>«tnn(s.t) 
JTiio«JTinHm 
ZiTN«Vl«Zl«ZtTl 
Z1TB>H2«Z1*Z1TI 
Z3TR*ni«Zl*Z2T« 
Z3TI«l1ll*Zt«Z3TE 
73  COITXlOE 
72  COITXIDE 
361*1 . 

JnSaRTI 

xnoaicn 

00  IM  1*1.171 
S61*-S61 
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aux*8n(i) 

XF(Sia.lQ.O)  M  TO  IM 

00  IM  8-l.SMU 

iS«ZTBr(JTB) 

M"XTBX(JTIB) 

M«nu«vs>ni«M 

«i>softMia 

Zl«tTB(S.a) 

ZlT»ill*Zl<»Z11V 

Z1TIMB«Z14Z11I 

Z3n>ts«zi<»z3n 

ZaTI««4*Zl«ZZTE 

mot .10.0)  00  TO  too 

«S«nBP(JTIO) 

VO>XTBZ(JTIO) 

JTlOaJTBKKTI 

«3«ni«vs>nii««t 

«1*S0B»«3 

«3>set*«4 

zi«im(s.t) 

ZlTN««l«Zt«ZlTR 

ZlTl«ll3«ZS«ZtTl 

ZZTNa-VS^Zt^ZZn 

Z2TI*>«4«Zt«ZZTI 

too  comioi 
too  cnruoi 

mit.n.o.oo.ot.Q.o)  00  TO  m 
iiTii«in)i«t 
TZ(nil)«-0«XtTll 
JtaRTIHKt 
TZ(Jt)>-0«ZaTR 
t7t  Jt«irri«KTR 
TZ(Jt)>-0«ZtTI 
Ja«Jt«Kt 
TZCJ2)«-0«ZZTI 
74  COOTIWl 

to  comm 

HOZTICZt .a04) (TZ(Z) .Z»t  .IZ) 
a04  rOBIUT('TXV(4tt4.7)) 
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CALL  BKaV(l3.XP8.T) 

CALL  SOLTlQD.ZPS.T.R.f) 

MtZTKai  .2M)  (V(Z)  .Z>1  .la) 
aoe  nttAT('T*/(4U4.r)) 

■ir»MqKr(Kaa-Hia(i))/a 

Aia-aD«(u-zz) 

cs-Dcai(Aaa) 

SMtzi(Aaa) 

iv*inKi 

si> .  iD*o*Mqtr  (HT/ncA) 

SA«Sl«f(IV)/(CS4U«0*8l) 

ciaoT«(ZLi<»i.D«o)«oc»u(es.sa)*SA 

ClZI*(ZLl-l.D«0)*OCIPU(CS.*Sa)*SA 

Hizncai.m)  cinrr.cin 
M  FaMUT(*cioor«*.aii4.7.'.  eizi*\ati4.7) 
Aio>iir«(La-zz) 
es>Mii(AM) 
si«oszi(Aia) 

Kf«CT4Cl 

SA«si«f(i:f)/(cs«ZLam«si) 
caaoT«(ZLa«i.D«o)*ocaru(a.sa)«SA 
cazi*(ZLa-i  .i>*4)«oc»u(cs.>si)«SA 
iiizTi(ai.7o)  caotrr.can 
70  nnuuT('caooT«*,aii4.7.*.  eazf«*.aEi4.7) 
ci(m»<iaoT*cau6(ciaoT) 
cizis>cizi«coije(cizi) 
caooT3»caooT*coiJc(caooT) 
caiis«caxa«coiJ6(cazi) 

PT«ciooTS-cizis<aooTS*cazis 
«UTi(ai.«a)  ciooTS.cizis.caooTs.cazis.FT 
sa  F0IUUT('Cl(RrrS«Mi4.7.«.  CtZIS«M14.7.'.  CaOOTS«'. 
lC14.7/*caZIS«Ml4.7.'.  n«Ml4.7) 

W7«0. 

SGMl. 

jm«o 

JTR0*0 

nM>. 

FXf>0. 

DO  ITa  l>l.tTR 
SCM- aOE 
sMAi-snid) 

XP(SaAZ.IQ.O)  CO  TO  173 
Al>t-1 
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00  170  s>i.auz 

«l«Xl<QiJ0(Zt) 

za«tnD(s.i) 

«M3<aiji(za) 
•T««r«(«i««2)«rai(s  .0) 
znMo. 


m«o 

00  174  lal.HUZ 

■a-i 

m>xi 

zrui.^.o)  M2m7 
n>ni(i) 

00  174  o-n.n 

■1*N-1 

7II1«I1 

ni«7xi/oc 

jniajni^i 

0S-XTK(J7n) 

Zr(Rl.lQ.O.Oft.ll.lQ.O)  CO  70  174 

RTRaim^l 

Jl-RTIHKl 

zt-7(ifni) 

Z3>T(J1) 

V3«rii*vs 

V1«SGR*II3 

ZTMI>-«l«Zl«H3«Z3«niB 
174  RTB-inX^l 
ja-RTi^nii 
J3>J3^1 
Z3>f(J3) 

Z4«7(J3) 

V4«FN1*HS 

V3«SeO*ll4 

rm—n*z2*V4»u*mii 

zr(Ki.iq.o)  CO  to  174 

JTNOaJTIRH} 

V4>XTR0(JTII0) 

XF(Rl.iq.O.Ot.ll.lQ.O)  CO  TO  174 

V3*FI1*«S 

V1-SGMV3 
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mo^*zi>«3*29»zno 
179  ««*nii*vs 

zno>ii3«z3m*Z4»znn 
17«  CWTIllll 
174  COITIIW 

xnB«.f*ini 

RaraKS.I) 

ni»4r<nja(iT) 

pa"nM*nn 

em^-n*mE 

zms«nnKaija(zm) 

cnis«cnB«cnjo(cna) 

ra»inM4znaB 

ra«ni4cms 

mnicai.iM)  t.s 

160  FaMUT(*l*M4.*,  S«*.Z4) 

nni(ai.i6i)  zna.xnn 

161  FaUUT(>ZTIII>*.ail4.7.'.  ZniS*M14.7) 

mzTi(ai.i6a)  cm.cms 

163  raUUTC'Cm**. 3614.7.*.  CXIIS>* .114.7) 

zr(fti.n.o)  00  TO  173 

zno«.i*rnn 

IT>tnB(S,l) 

FXIS»IT*C0IJ6(fT} 

PXI*PXI4FXIS 

Cm»iT4ZTN0 

ZT«05>ZTll0«C0fJ6(ZTM) 

CTM0S<TN0«C(MU6(CTN0) 

PM*PW4ZTII0S 

pft-pt4cnRS 

vazn(31.1i7)  ZTNO.XTMH 
167  POUUTC'ZTW-'. 3114.7. *.  ZTinS*' .114.7) 
wun(31.l60)  CTW.CTMOS 
lit  P0I1UT(*CTN0«*.S14.7.*.  CT1B6«'.I14.7) 
173  COWXIW 
173  COITXIOC 
S66-1. 

JTII>0 

JTtOmO 

DO  180  lal.in 
DOMm-DBH 

siuz«sri(i) 

XP(SIUZ.IQ.O)  60  TO  160 
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DO  Itl  S-l.SHU 

»-ITB(S.I) 

viaZKOijeczi) 

Z>iTB(S.I) 

i3«Z2<aiJ0(Z2) 

ZTBIoO. 

ZTB0>0. 

MnNO 

DO  113  ■■l.niiZ 

M3-1 

ll«f-l 

rai«ii 

ir(fl.BQ.O)  IQ«3 
■3«nf(i) 

DO  rs  N-R3.10 

FMl-Kl 

FXl«ni/BC 

v$«nzEr(jnE) 

V6«mEZ(JTiX) 

ZF(Nl.EQ.O.Ot.ll.EO.O)  CO  TO  ITS 

Jl-NTIHKl 

Z1>T(KHI) 

Z3>T(J1) 

V3«nii*v6«nii«v6 

vt-ses*H3 

ZTEB>Vl«Zl^V3«Z3«Zm 

ITS 

J3-KTB^KTR 

JS-JZ-HCl 

Z2>T(J3) 

Z4«?(J3) 

V4«nii«vs«ni«ira 

V2«SeB*V4 

ZTZB-«3*Z3«W4*Z4«ZTEZ 
ZF(Bl.Eq.O)  CO  TO  TS 
JrE0«JTE0^1 
VS-ZTBOrCJTtO) 

«6>ZTE0Z(JTI0) 
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XrCMi.lQ.O.aK.ll.lQ.O)  00  TO  its 

ii3«rai*wt>Fiii*iie 

«i-so>*fia 

m0>Hl«Zi*ll3«ZS*ZTI0 

Its  v^nt*m*ni*w 
«2>sa*v4 

ZTHMI2«Z3HM«Z4«ZT(0 
7S  COWUUl 
Its  COITXIDI 
ZTB>.S«ZTII 
BT-RB(S.E) 

PXlS>BT«COIiQ(BT) 

PZI>PIl4PZtS 

cmB«BT«ms 

ZTIBS«ZTIB*C01J6(ZTEI) 

CTiis>cm*coije(cm) 
pitapM^rrEts 
Pt«Pft4CTIBS 
ttZTKSi.sto)  t.s 
vtZTKSi.iM)  zm.niis 
ISO  POMUTCzm-.siM.r.*.  nns-Mio.T) 
vKmoi.stt)  cm.  ems 
its  pollute *€111- *.3114. 7.*.  ems«Mi4.7) 
XP(ll.lQ.O)  cn  TD  til 
mo>.i«mo 

■T>ITID(S.E) 

pzis-BT«caije(iT) 

PZI>PZI«PXIS 

cTEOa-iT^nn 

zTios>zTioKaijc(im) 

CTE0S«CT10*CaiJ«(CTn) 

PIII«PIH«t7IOB 

PE*PE«CTIOS 

vEznoi.iM)  mo.smi 

tee  POElUT('nO>*.Sl4.7,'.  IT»«*.tl4.7) 
vEZTKsi.to)  em.enot 
to  P0EIUTCCn0«*.Xt4.7.*.  Cm*Ml4.7) 
iti  COITZlOE 

ISO  comioE 

VEZTKSt.liO)  PXI.PM.Pt 

190  P0EIUT(*PXt*Ml4.7.*  Pt»»Ml4.T.*  PE>Mt4.7) 
nUPLT(U)*IIS 
PZIC(KA)>PXI 
PTEA1(E1)«PT 
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ra»L(Ki)«m 
rTOTAL*mni 
«tni(2i.is4)  rroTAL 

tS4  FaUUT(*PT0TAL-M14.7) 

TBM 

TOM 

00  91  Z«l.Ka 

Tii>ccija(v(z)) 

Tn-Tiii4Tue(z)«f(z)«m 

Tn»tMi*TZ(Z)*TB 
91  C09TZ101 

m— AZMiGCTn) 
nMi«-PTA-AZRAG(Tn) 
mzn(2i.93)  rn.puu 
93  F0MUT(*m«*.E14.7.*.  PUU*S314.7) 
ZP(KA.n.IX3(IlI))  CO  n  49 
IAI«UI«1 
DtL»n/(Mni~i) 

00  109  J-l.lFlZ 
I3A1P(J)«0. 

■3A3P(J)>0. 

PJ1«J>1 

PJ10«FJ1«0BL 

pni(j)«pz>rjio 
FlZ2(J}*rJlO 
106  COfTZlOB 
DIL«PZ/(n«l) 

00  loe  >i.R 

l3AlZ(J)-0. 

l3A2Z(J)-0. 

FJl-J-l 

Z(J)«FJ1«0EL 

106  comm 

ZFdAl.n.Z)  60  TO  lU 
mzn(21 . 162)  (FB2(Z)  .Z-1  .IFIZ) 

162  F0llUT('Fn2V(ai4.7)) 
«UTI(21.161)(Z(Z}.Z«t.R) 

161  F0IIUT('Z*/(6B14.7}) 

166  SA>-PZ2«0SgtT(FZ«9C/ir7) 

00  136  J-l.FMAZ 
FJl-J-1 

SZIF(J)>0SZI(FJ1*FZ6)«6A 

SZIQ(J)«SZIF(J)«R 

136  comm 
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JTINO 

m^o 


00  iol  q-i.mz 

roiaqi 

pqic>iqi/c 

M-l 

xr(qi.iq.o)  n-o 

n*ni(q) 

sziqq-sxiq(q) 

DO  104  r>n.F3 

JTi»jn+i 

JRlajn^KTI 

MU-l./BnKJTI) 

P1«P-1 

FPl-Fl 

ZP(qi.Eq.o)  Biii>Bn>sq3 
ZF(Pl.iq.O)  BlllJ«BKIJ*sq3 
SXIPP-SIIP(P) 

Biatjp>Biiij*rpi/B 

Biajq>Biiij*pqic 

injpp>iiaijp«szipp 

Bnijqq>BiiiJQ«sziqq 

Tiai«BiaijqQ*T(jTEi) 

TEll>BiaiJPP*V(JTBl) 

THE— BIIIjqQ*V(JTO) 
TCS-BIIBJPP«V(m2) 
ZF(Pi.Bq.o.oB.qi.Bq.o}  co  to  134 

jm-jTii+i 
JT1I3>JT1I4K1 
BiaiJPQ-BIIBJP«SZBQq 
Biaijqp>BMBjq«szBPP 
TIIM>-BIIBJPq*T(J11f)^TIOI 
TBII»BIIBJQP*f<JTH)>TEII 
TlIBaBniJPO*T(JTIQ}^TIfE 
TEB-BiaiJqP«T(JTie)<»TEE 
134  PT-FPl 


DO  136  J«1.BP1Z 

l361P(J}aE3AlP(J)<»TI0l*C0S(PT«PBZl(J}} 
E362P(J)«E3i2P(J)-»TIIE«C0S(PT*PBZ2(J)) 
136  COBTZBDE 
PT-Fqi 

DO  148  J«1.BZ 
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P1*C0S(PT*Z(J)) 

isiiz(j)«miz(j)’»TBM«ra 

l3A2Z(J)-B3A3Z(J)-^TBB*Pt 
i  148  C08TI1UB 

104  COITZIDB 
101  COITZIDS 
.  00  137  Jal.IPIZ 

I3A1?S( J)*CABS(E3A1P( J) ) 

I3A3P8( J)«CABS(ISA2F( J) ) 

137  COITZIDI 

00  138  J-l.IZ 

t3AlZS(J)«CAlS(B3AlZ(J)) 

l3A2ZS(J)-CABS(E3i3Z(J)) 

138  COITZIOB 

HUTB(21 . 140)  (B3A1PS(  J) .  J>1  .IPIZ) 

140  F0BMAT(*E311PS'/(SI14.7)) 
nOTBCai .  141)  (I341ZS(  J) .  J>1.R) 

141  F0UUT(*E3ilZS*/(Stl4.7)) 
inZTB(2l .  142)  (■342n(  J) .  J>1  .IPIZ) 

142  F0MAT(*E3i2PS*/(SI14.7)) 
iniZTB(21.143)(t342ZS(J)  .J«1.R) 

143  F0llUT(*E3A2ZS*/(Stl4.7)) 

48  COITZIOS 

nZTB(21 . 140)(BUn.T(Z)  .2*t.U«) 

140  F0tlUT(*BKAPl.T’/(W14.7)) 

VUTB(21 . 188)  (PZIC(Z) .  X*!  .IA8) 

189  F0BIUT('PZIC'/(IS14.7)) 

VErrE(2l  .R)(r7UI(Z)  .Z«l  .148) 

92  F0BlUT('rnulV(Kl4.7)) 

inaTE(2i.i39)(Pun.(X)  .z«i.Z4a) 

139  F08IUT('PUn.'/(8El4.7)] 

STOP 

ElO 


131 


Appendix  A 

Verification  of  (3.119) 


The  term  in  parentheses  on  the  right-hand  side  o(  (3.119)  is  |pven  by  eqs. 
(4.8)  and  (4.9)  of  [2]:^ 


„  _y 


m  even 
m  odd 

m  even 
m  odd 


(A.1) 


(A.2) 


where  is  Gj**  when  r  =  0  and  s  =  1.  The  left-hand  side  of  (3.119) 

is  given  by  (3.116): 


0  _  (A.3) 

where  it  is  understood  that  is  evaluated  at  (r,«)  =  (0,1).  In  Ap¬ 

pendix  A,  we  verify  (3.119)  by  showing  that  the  right-hand  side  of  (A.3)  is 
equal  to  the  product  of  ^k/0^  with  the  right-hand  sides  of  (A.l)  and  (A.2). 
Substituting  (3.110)  into  (3.57),  (3.58),  (3.63),  and  (3.64)  and  multiplying 

*Tbe  right-hand  tide  of  eq.  (4.8)  of  [2]  it  a  poor  ezprettion  because  it  it  indeterminate 
when  m  s  0.  Bowevet,  the  riglu-band  ride  of  (4.8)  of  p]  should,  at  it  evident  from  its 
derivation  in  [2],  vankh  when  m  «  0. 


a 


A 
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the  resulting  equations  by  —  we  express  the  right-hnnd  nde  of  (A.3)  ns 

-77r"''™*.->nr™*2[(}r*]^,  0-1,2  (A.4) 

«-».*  (AS) 

where  T™  is  gjven  by  (3.60)  with  r  =  0  and  a  s  1.  In  (A.4)  and  (A.5), 
is  given  by  eq.  (6.29)  of  [2]  with  p  replaced  by  m  and  with  r  s  0: 


(A.6) 


Setting  r  =  0  and  p  =  m  in  eq.  (3.41)  of  [2]  and  using  eqs.  (3.44)  and  (3.45) 
of  [2],  we  have^ 


”'ZY 

,  m  odd  J 


(A.7) 


In  view  of  (3.89),  substitution  of  (A.7)  into  (A.6)  gives 

Substitution  of  (3.60)  and  (A.8)  into  (A.4)  and  (A.5)  gives 

^TM.TM« 


aM  _  1*1  /  1 

"*’  »  1  i’  j* 


(A.8) 


Comparing  (A.9)  and  (A.IO)  with  (A.l)  and  (A.2),  we  see  that  expressions 
(A.9)  and  (A.IO)  for  the  right-hand  side  of  (A.3)  are  indeed  equal  to  the 
product  of  with  the  right-band  sides  of  (A.l)  and  (A.2). 


*It  is  evideot  from  eq.  (C.IS)  of  (1]  that  dm*  should  vanish  when  m  s  0. 
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Appendix  B 

Verifi''.ation  of  the  Equality  of 
Modal  Coefficients 


In  Appendix  B,  we  show  that  expressions  (5.21),  (5.31),  and  (5.32)  for  the 
modal  coefficients  C^,  and  i^uce  to  the  right>hand  sides  of 
eqs.  (6.88),  (6.96),  and  (6.99)  of  [2]  when  s  1  » the  only  nonaero  B  in 
(3.1)  and  when  only  the  TMqi  and  TEu  mo^  can  propagate  in  the  drcular 
wavyiide.  The  V'$  in  eqs.  (6.88),  (6.96),  and  (6.99)  of  [2]  are  the  products 

of  with  the  V  ’s  in  the  present  report  because,  as  stated 

in  Section  3.5,  the  excitation  in  the  present  report  is  that  of  [2]  multiplied 
by  Multiplying  each  V  in  eqs.  (6.88),  (6.96),  and  (6.99) 

of  (2)  by  using  (3.27)  to  replace  this  multiplicative  factor 

by  and  using  eqs.  (6.89),  (6.90),  (6.97),  (6.98),  (6.100), 

and  (6.101)  of  (2],^  we  convert  expressions  (6.88),  (6.96),  and  (6.99)  [2]  to 

expressions  in  ternu  of  the  V's  in  the  present  report: 


respectively. 
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Now,  the  objective  is  to  show  that  the  right-hand  sides  of  (5.21),  (5.31), 
and  (5.32)  reduce  to  those  of  (B.l),  (B.2),  and  (B.3)  when  =  1  is 
the  only  nonzero  B  in  (3.1)  and  when  cmly  the  TMoi  and  TEn  inodes  can 
propagate  in  the  circular  waveguide.  Setting  (r,s)  =  (0,1)  and  Bj^  s  1 
in  (5.2n,  substituting  the  modified^  right-hand  sides  of  (3.57)  and  (3.63)  for 
*,  substituting  the  modified  right-hand  sides  of  (3.58)  and  (3.64)  for 
^te.TM€^  and  using 

-  *?’.  (B.4) 

^Heie,  “modified"  mean  with  (m.n,t)  repUced  by  (p,«,  j).  In  (3.57),  (3.58),  (3.63), 
end  (3.64),  is  given  by  (3  110)  in  which  4^  =  "d  is  given  by  (3.60). 
Equation  (3.60)  contains  of  eq  (6.60)  ^  (3). 
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weobuin 


Equation  (B.4)  was  obtained  by  setting  r  =  0  in  (5.17).  Setting  (r^)s(l,l) 
and  s  0  in  (5.31),  substituting  the  modified  right-hand  sides  of  (3.71) 
and  (3.79)  for  substituting  the  noodified  right-hand  sides  of  (3.72) 

and  (3.80)  for  and  using  (3.112),  (3.113),  (3.75),  and  eq.  (6.60)  of 

[2],  we  obtain 


Setting  (r,s)  s  (1,1)  and  b  0  in  substituting  the  modified 

right-hand  sides  of  (3.81)  and  (3.M)  for  substituting  the  modified  ^ 

right-hand  sides  of  (3.82)  and  (3.87)  for  and  using  (3.114),  (3.115), 

(3.75),  and  eq.  (6.60)  of  (2],  we  obt«n 


e+f#o 
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Equati<»u  (B.5),  (B.6),  and  (B.7)  are  identical  to  (B.l),  (B.2),  and  (B.3), 
reapectivdy.  Hence,  our  objective  has  been  achieved. 
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